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I) Introduction 

This program studied the dynamics of the solar photosphere and 
chromosphere. Observations obtained by the Laboratoire de Physique 
Stellaire et Planetaire*s ultra-violet spectrometer onboard the 0S0-8 
spacecraft were analyzed and two papers presenting results were 
published. Three reviews on different topics concerning the theory 
of stellar atmospheric dynamics were presented and published. Dynamical 
models of the chromosphere and the emitted resonance line spectrum were 
calculated, and one article has been accepted for publication and 
another is in preparation. Three contributed papers were delivered 
at scientific conferences. These publications are attached. 



II) Discussion 

Repeated spectral scans of strong, optically thick resonance lines 
formed in the solar chromosphere were analyzed for indications of oscillatory 
velocities and intensities. Figure 1 displays Lyman alpha data obtained 
during one orbit, which has had data losses replaced by interpolation. These 
spectra were taken away from active regions, near the center of the disk, with 
a spatial resolution of 4 x 10 arcseconds. They show only the central 
absorbtion reversal and one of the adjacent mission peaks, and succeeding 
scans are displaced vertically for viewing "clarity" . One can see the geo- 
coronal Lyman alpha absorbtion near grating position 898 and the blue peak 
near 928 (wavelength increases to the left) . A slow displacement of the 
profile to the left can also be seen. Figure 2 shows the average profile 
for this orbit, after the correction for the wavelength drift has been applied. 

Among other indicators of velocity which were studied, the blue peak 
is reasonably well defined, and the position of a parabola fitted by the method 
of least squares was used to define it. Figure 3 shews a series of spectra 
and the fitted parabolae, and Figure 4 shows the resulting variation of the 
peak's position. A sine wave, whose period equals the satellite’s orbital 
period, has been drawn through the measurements, and the positions of data 
losses indicated by arrows. The residuals are shown on Figure 5, and (apart 
fran the measurment at 560 seconds where the data was fit by a parabola with 
a minimum, rather than a maximum) one sees the indication of a variation with 
a period near one thousand seconds. All of the Lyman alpha data had this same 
characteristic oscillatory behaviour. Figure 6 presents the variation of 
the intensity of the blue peak and the curvature of the peak. 
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After the first flush of excitement over the discovery of this "900 
second" oscillation - for that was it.? mean period - had passed and references 
to other observations of chromospheric oscillations with similar periods were 
dredged out of the literature, the real discovery was made that the observed 
periods corresponded precisely to the beat period between the repetition rate 
of the spectra (10.24. seconds) and the rotation rate of the spacecraft (10 
seconds on the average, with slow drifts of several seconds) . The amplitude 
of the intensity fluctuation is too large to be accounted for by displacements 
of the pointing axis, and it appears that it results fran a modulation of the 
transmission with the period of the spacecraft rotation. Extensive additional 
observations and analysis by the P.l. team have demonstrated that the effect 
does not manifest itself in the other channels of the spectraneter, and the 
intermittency of the effect and its lack of correlation with the operational 
mode of the LPSP instrument suggests that it arises fran another canponent 
of the spacecraft. Efforts to identify and correct for the modulation have 
been unsuccessful to date. 

More positive results using statistical methods are contained in 
attachments B, C, and D. 

A major effort of this program was devoted to the dynamical modeling of 
the solar chranosphere with the dual goals of improving our understanding of 
1) the dynamical processes themselves and 2) spectral line formation in the 
dynamic chranosphere. As these lines provide the principal diagnostic tool 
for chromospheric measurements, the calculation of the diagnostics combined 
with the complete knowledge of the physical conditions giving rise to these 
diagnostics gives us a means of evaluating the efficacy (and veracity) of 
the measurement techniques. The combination of hydrodynamics and radiative 
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transfer in this way offers sane exciting new insights into the behaviour of 
the Sun as we see briefly below. The results of this work will be published 
in two papers in collaboration with R.F. Stein of Michigan State University 
and P. Gouttebroze of LPSP. 

A one dimensional, non-linear dynamical model of the solar atmosphere 
fran 1.6 megameters below the visible surface to 3 megameters above was 
excited fran below by pulses or oscillations in a series of calculations 
performed at LPSP. Figure 7 shows the velocity (solid lines) and pressure 
(dotted lines) at a number of different altitudes as a function of time for 
a typical calculation. The important points to note are the standing wave 
(pressure and velocity 90° out of phase) formed in the lower atmosphere (the 
zero of altitude corresponds to the visible surface) and the less regular, 
shorter period, standing wave formed in the chromosphere (1200 to 2000 
kilaneters) . The photospheric oscillation is the well known "300 second" 
oscillation, and the chranospheric oscillation is the - slightly ±ess - well 
known "200 second" oscillation. Variation of the structure of the chromosphere 
and the amplitude of the motions has shown that the chranospheric oscillation 
results fran a trapped, interfering wave - much as the 300 second oscillation 
does - rather than fran a ringing at the cut-off for propagating acoustic 
waves, as had been previously thought. 

Figure 8 shows a series of magnesium resonance line profiles calculated 
for a dynamical model similar to those in Figure 7. A rather low amplitude 
motion was studied so that differences fran the better understood, static line 
profiles would not be too great for an initial study. Displacements and intensity 
fluctuations of the emission peaks as well as of the central absorption reversal 
are evident. Figure 9 shows the variation of intensity at different wavelengths 
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as a function of time. The photosphoric 300 second oscillation dominates 
further than 0.5 Angstroms from line center, while closer to the center the 
200 seoond oscillation dominates. The intensity variations of the emission 
peaks as well as the displacements of the central absorbtion stand out clearly. 
The intensity fluctuation at each wavelength has been expressed as a fractional 
change about the mean at that wavelength, and the fluctuation at each wavelength 
has been normalised to unity. The fractional fluctuation at each wavelength 
is indicated at the right. In the blue peak, the intensity varies by over a 
factor of tvro, even for this model with its relatively low chromospheric 
velocities. 

Figures 10 and 11 show the results of fitting parabolae to the maxima 
and minima of the profiles. The central absorbtion is referred to as K3 
and the emission peaks as K2 (violet and red) . Because of the standing wave 
character of the dynamics, the perturbations remain in phase through a large 
range of heights in the chromosphere, and thus the two enission peaks are 
of equal intensity nearly simultaneously with the maximum cf K3, and the zero 
crossing of the velocity. The very high vertical phase velocities result in 
diagnosed velocities which are remarkably close to those at the heights of 
formation of the spectral features. Perceptible differences in the phase of 
the three diagnosed velocities are nonetheless apparent. 

The contribution to the emergent intensity as a function of altitude 
is shown in Figure 12, as a function of time, for a wavelength justly slightly 
to the blue of line center, within the central absorbtion feature. The variation 
in altitude of the emitting region results primarily frcm the motion of the 
matter up and dcwn during the oscillation. Figure 13 plots the same contribution 
functions against the average height of the matter, so that the abscissa is 
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just the mass depth transformed to read in kilometers. The emission arises 
from essentially the same material independently of the time. One should 
note that only a very, very small range of altitudes contributes. The 
behaviour at the blue emission peak is slightly more surprising (Figure 14) . 
Reasonably large variations of the altitude of the emitting material occur, 
which results from the doppler shifts of the absorption profile with very 
rapid wavelength dependence. Figure 15 shews another projection of the same 
information. The depth above which 0.1, 0.2, 0.3 ... 0.9 of the emission arises 
is plotted for the same wavelength as Figure 12, as a function of time. The 
variations of the thickness as well as the altitude of the emitting region are 

striking. 

Finally, Figure 1C shows the contribution functions on the flanks of the 
emission peak where a dramatic alternation of the emitting material by over 
1000 kilometers occurs. This is but the beginning of the understanding of 
this great wealth of information which premises to not only increase our 
knowledge of the formation of spectral lines and their use as diagnostics, but 
also to enable us to make testable predictions on the bases of various 
hypotheses concerning stellar atcmspheric dynamics. 

This vrork would not have been possible without the generous collaboration 
of the staff of the Iaboratoire de Physique Stellaire et Planetaire. Access 
to the data, even after it had been taken, was possible only as a result of a 
tremendous effort on their part. The non-linear, hydrodynamic calculations 
and the radiative transfer calculations were carried out on the CDC 7600 
computer of the Centre National d* Etudes Spatiale. The data analysis 

and analysis of the results of th dynamical ar.d transfer calculations were 
carried out within the Space Astronomy Group of IMSC and profited greatly 
from their expertise. Finally, I wish to thank NASA for its support of this 
research program. 
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Attainment A 

Published in the Proceedings of the November 1977 OSO-8 Workshop 
Boulder, Colorado; pages 311-339 

Theoretical Aspects of the 300 Second and Related Oscillations 

of the Solar Atmosphere 

John Lcibachcr 
Space Astronomy Group 
Lockheed Palo Alto Research Laboratory 
Palo Alto, California 


I. Introduction 

It seems to me that only extremely rarely does our knowledge of the 
behaviour of the Sun, and astronomical objects in general, advance in a nice, 
clear fashion at all similar to The Scientific Method, whose virtues were 
extolled back in high school physics courses. Only infrequently do real, 
testable predictions emerge from hypotheses formulated to describe known, 
but unexplained phenomena. Even less frequently is the physics sufficiently 
simple that one effect dominates and the predictions are unambiguously confirmed. 
Nonetheless, this has occurred recently in the interpretation of what for a 
long time was known as "wiggly line spectra" - the spatially and temporally 
resolved velocity field of the solar photosphere. It is thus with nostalgic 
pleasure that 1 undertake a description of the denouement of nearly twenty 
years of research into the nature of the "five-minute" oscillation of the 
solar atmosphere. 

After briefly recalling what the problem was, I shall review the model 
of the "five-minute" oscillation which has demonstrated this predictive 
success - first, the atmospheric properties which determine the resonant 
frequencies, the variation of these frequencies with horizontal scale and the 
vertical structure of the resonant motions; and second, the mechanisms which 
have been proposed to drive these motions. While the detailed, numerical 
agreement between theoretical predictions of the resonant frequencies and 
observation has not only confirmed the model but allowed its application as 
a diagnostic of sub-surface conditions ("solar seismology"), our ideas about 
the excitation and damping mechanisms are in a much more preliminary state. 
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Theft, I shall discuss the chromosph- ric concomitants of the observed pholo- 
•pheric notions. (Given the input of the photosphcric oscillations, what do 
we expect to see in the chromosphere as a result ?) Although the mechanism 
for the decrease in energy of the "five minute oscillation" and the increase 
In energy of th t "three minute oscillation" appears to be understood, attempts 
to model the temporal variation of optically thick lines have been lnrgely 
unsuccessful. Although it is possible that our dynamical description is in 
error, it seems clear that the proper combination of hydrodynamic models and 
line formation diagnostics has net yet been used to derive a- riate 

predictions for comparison with observation. Finally, we ar<. -ed to 

conclude that the most naive expectations of shock wave hcati..^ wx the solar 
chromosphere and corona have not been fulfilled, and that while the observational 
quest must continue - guided by more realistic theoretical predictions - 
alternative heating mechanisms should be actively investigated. 


II. Oscillatory Motions of the Solar Photosphere 

Although the oscillatory motions of gravitationally stratified, compressible 
fluids - atmospheres - have been "well understood" since the work of Lamb at the 
beginning of the century, in the absence of sperific problems posed by 
observation, astronomers seem not to have considered the consequences of the 
existence of such hypothetical motions. While Biermann and Schwarzschild had 
conje 'ured that running waves generated in the convection zone could trans- 
port the energy required to account for departures from radiative equilibrium 
in the outer layers of the solar atmosphere, it was assumed that these waves 
would span a broad spectrum of temporal and spatial scales commons urate with 
the "turbulent" convective motions thought to be responsible for generating 
them. When Leighton actually measured the time variation of the velocity at 
a point on the Sun, he found much to everyone's apparent surprise that instead 
of being a more or less random function, the velocity oscillated quite 
regularly with a period near five minutes. What is so special about five 
minutes? How docs the Sun go about selecting just this one narrow range of 
periods? 
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It turns out, embarrassingly enough, that there are a lot of ways to 
do it - none so much as suspected before Leighton's observations directed 
attention to the problem - and within a decade half a dozen distinctly 
different, but "credible" models for the five minute oscillation emerged. 

(The Interested reader is referred to reviews by Michalitsanos(1973) and Stein 
and Lelbaeher (1974 )>. Within the last two years observations have been obtained 
(Deubner, 1975) and corroborated (Rhodes, Ulrich and Simon, 1977) which 
confirm predictions made on the basis of one of these models, and thus in the 
current order of things this model has been elevated to the status of 
"astrophysical truth" - not to be confused with "astrophysical accuracy." 

Let me Qualitatively describe the important aspects of this model. 


The "truth" is that acoustic waves with periods greater than 200 seconds 
will be trapped within a limited region, some tens of megameters thick, just 
below the visible surface of the Sun where they couple to the thermal and 
dynamic state maintained by the outward convective transport of the solar flux. 
These trapped waves drive non-propagating, evanescent waves seen in the visible 
atmosphere. To describe this phenomenon, it behooves us to think a little 
about acoustic waves. 


A change in the density of a gas, brought about by pushing on one end 
of it for example, generally causes a change in the pressure through the 
conservation of energy Since the pressure is just the internal energy densit 
compressing the gas increases the pressure. Since the pressure has not 
increased everywhere, only on the side where we pushed, a pressure gradient 
now exists which will in turn accelerate the gas towards the low pressure. 

This acceleration compresses the low pressure gas and we've gone around the 
loop and start once again: a density increase (via the conservation of energy) 

causes a pressure increase, a pressure gradient (via the conservation of 
momentum) causes an acceleration, an acceleration (via kinematics) causes 
a displacement, a differential uisplaceroent (via conservation of mass) causes 
a density increase. This rompressional disturbance is propagated by the 
internal random, thermal motions of the particles at a speed — /kT/m, which 
la referred to as the speed of sound. 
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la fact, in an infinite , isothermal, homogeneous gas, the sound speed 
represents the only characteristic dimensional quantity, and it is just the 
signal propagation speed for small disturbances. If we look at a plane 
parallel sound wave, the wave crests move at the sound speed along the direction 
<f propagation, and hence they must rave faster than the sound speed along any 
other direction. (It's worth drawing a little sketch.) For example on a 
plane nearly perpendicular to the direction of propagation, the crests will 
move sideways at "nearly infinite" speed. The propagation speed of the crest 
(or trough, or node) is often referred to as the phase velocity. Interesting 
things occur at the interface between two regions of differing sound speed 

(i.e. diff< .ent temperature) . Waves which are propaging towards the region of highc 
sound speed with an arbitrary period and wavelength (just the product of the 
period and the sound speed) will maintain their period and wavelength along 
the interface (that is, the phase velocity along the interface is the same 
on both sides of the interface). But in the region of higher sound speed, 
the wavelength corresponding to a given period must increase and hence the 
wavelength perpendicular to the interface must increase (since the wavelength 
parallel to the interface is fixed by continuity). This is just the same as 
saying that the direction of propagation changes, in the sense chat the wave 
propagates more nearly parallel to the interface. Which is all just 
e.,aivalent to Snell's law in optics - waves propagating into a region of 
decreasing index of refraction (increasing speed of light) arc refracted away 
from the normal. The "interesting things" happens when the temperature 
increases so that the sound speed in the higher temperature region equals, 
or exceeds, the horizontal phase velocity of the incident wave. The wave 
can no longer propagate perpendicular to the interface in the hot region, and 
it is reflected ("total internal reflection"). In general, the wave doesn't 
atop abruptly at the interface, but it's energy falls off exponentially into 
the hot region. This "evanescent wave" is curious in that although at any 
point the velocity varies sinusoidally in time, there is no phase propagation 
away from the interface. The velocity reaches a maximum ( or minimum) every- 
where simultaneously, so it might be more appropriate to think of the phase 
velocity perpendicular to the interface as bding infinite. Furthermore, in 
contrast to the propagating wave 'here the pressure and velocity fluctuate 
in phase * which is just what is required *c ?t a- aunt of work to be done 


4 





•round one cycle - the pressure end velocity are ninety degrees out of phase, 
so that no energy is transmitted. 


A wave propagating at some angle to the temperature stratification in the 
direction of increasing temperature will eventually get to a region (temperature) 
where the sound speed equals the horizontal phase velocity and the wave will be 
turned around. It is worth noting that if we look at waves of the same 
frequency but varying. angle of propagation with respect to the temperature 
stratification, or equivalently varying horizontal wavelength if we think of 
• "vertical" temperature strat-ification, the more nearly horizontally 
propagating waves will be reflected first (at the lowest temperature) and the 
sound spped to which the waves can propagate is iust proportional to the 
horizontal wavelength; or inversely proportional to the sine of the angle to 
the vertical. Thus waves propagating "straight in", that is waves with infinite 
horizontal wavelength, will never be reflected. The existence of such re- 
flections is of interest in that it can exclude waves from the high temperature 
region (the deep interior. of the sun) and, more significantly in the present 
context, in that two reflections will form :a cavity, or duct (an organ pipe 
or cello string) in which interference will occur and certain frequencies will 
be selected from a broad band excitation. For example, one of the first 
models for the "five minute oscillation' 1 - proposed by F.D. Kahn but later 
shown to be unacceptable on several points - was based upon the trapping of 
sound waves near the temperature minimum by the higher temperatures above and 
below. 


Other "interesting things" happen to sound waves when they propagate 
in a region of changing density - such as occurs in the presence of a gravi- 
tational field, but which would also occur if the mean molecular weight were 
stratified for example. If this "atmosphere" were displaced upwards or down- 
wards "slowly", then sound waves would run back and forth, and after they damp out 
and the atmosphere would return to its stratified equilibrium. That is all 
well end good, but what is "slowly"? Whereas the isothermal. gas we considered 
above had only a characteristic velocit associated with it, the density 
stratification introduces a characteristic length - the scale height. If we 
move the atmosphere periodically up and down with a period sufficiently long 
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to that a small density fluctuation would travel many ("many" turns out to 

equal in) scale heights in one period , then the atmosphere more or less stays 

together and moves up and down in phase. However* if we jiggle the atmosphere 

up and down very rapidly - more rapidly than the atmosphere can respond - 

then waves sip on up at the sound speed. Since the sound speed remains constant, 

these waves propagate without reflection. However* since the atmospheric 

density varies* the velocity amplitude must vary inversely to maintain the 

2 

constant flux of energy (pv c, where c is the velocity of sound). Thus waves 
propagating upwards* into less dense ga& will increase in amplitude as a 
As the period increases to that at which the atmosphere can respond in phase - 
the cut-off period for propagation - the energy propagation velocity decreases 
to zero and the wavelength increases to infinity. For periods slightly longer 
than the cut-off, the atmospheric fluctuations vary exponentially rather than 
Sinusoidally with height* and these waves arc again referred to as "evanescent." 

If the gas is stratified gravitationally, the scale height is inversely 
proportional to the temperature, so the time for a density disturbance to 
propagate a scale height is just proportional to T//T or /l. Once again, 
the propagation characteristics change with a change in temperature so 
that when a wave propagates into a cooler part of the atmosphere - with 
a smaller scale height - waves of sufficiently long period will no longer 
propagate, that is they will be reflected. Thus a sound wave propagating 
into a region of continuingly decreasing temperature will eventually reach 
a region where its period times the sound speed equals i*~. times the density 
scale height, at which point energy will be propagated no further. If sound 
waves covering a range of periods propagate into a region of decreasing 
temperature, the longest period waves will be reflected first and the short 
period high frequency waves will propagate down to low temperatures (small 
scale height). Combining this behaviour with the reflection that occurs 
at high temperature (when the sound speed equals the horizontal phase speed) 
we see (I hope) that in an atmosphere there exists a range of temperatures 
within which a given acoustic mode (specified pair of period and horizontal 
wavelength ) may propagate. Whereas we have already seen how a cavity capable 
of trapping acoustic waves is formed at the temperature minimum, we now have 
discovered that there exist additional cavities in the rising temperature 
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regions above and below the temperature minimum. In fact, the second model 
for Che "five minute oscillations" - proposed by Bahng and Schwarzschild but 
alas, also shown to be unacceptable on several points • suggested that waves 
trapped in the chroaosheric temperature rise could account for the observations. 

There followed many other descriptions; some using running acoustic 
waves, others buoyancy waves trapped in the temperature minimum; before the 
’mportance of acoustic waves trapped in the temperature rise towards the 
solar interior was pointed out by R. Ulrich and independently by Bob Stein 
and myself. As simple as this description now seems, it took an embarrassingly 
long time to arrive at it. In my own case, 1 was considering a very different 
aspect of the problem and was extremely annoyed that the model developed 
beautiful five minute oscillations long before it had a chance to develop wnat 
I wanted to study. After a lot of hard work trying to make the oscillations 
go away, I conceded defeat and accepted what the model was trying to tell us. 
Ulrich recounts a similar serendipitous experience, having originally thought 
that the "overstability", or growing amplitude of the oscillations was a 
numerical difficulty of his computer program. Sut I an getting ahead of the 
story. Let us consider in a little more detail the properties of the sub- 
photospheric cavity, before discussing how these trapped waves are excited. 

For any given wave mode (period and horizontal wavelength), in a propagating 
region, there exists two waves corresponding to upward and downward propagation. 
If a reflection exists somewhere, then the wave propagating away from the 
reflection is just the wave propagating towards it, but turned around. Of 
course, the wave experiences a phase chance at the reflection. This would 
be the situation in a semi- infinite uniform gas with a wall at one end, for 
example. If we introduce a second reflection the situation becomes much more 
intersting. Consider an excitation operating within the cavity; a diaphragm 
moving up and down for example. For an arbitrary periodic excitation, the 
waves reflected back to the excitation will arrive senewhae out of phase 
with it and will impede the excitor. Thus the gas may be moving upward at 
the diaphragm while the diaphragm is trying to push the gas back down. Not 
a very efficient way to transfer energy, I hope we all can agree. For certain 
periods of excitation, however, the reflected wave will arrive back at the 
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driving exactly in phase with it and the excitor can put more energy into 
the gas and the wave will grow in amplitude . That is, the gas would be caving 
with the diaphragm. A slightly more concrete analogy, would be trying to push 
a child on a swing. If you don't push in phase with his swinging, it doesn't 
work very well. These resonant modes, often referred to as "normal modes", 
arc of interest to us here just because of their resonant response to broad 
band excitation. That is, while waves spanning a broad range of periods ... J 
horizontal wavenumbers will be trapfwrd in the sub-photosphcric cavity, these 
resonant modes should be preferentially excited to large amplitudes. 

In general, a cavity formed by two reflections will have many resonant 
modes, that is many wavelengths that can satisfy the rcquirenet.u that after 
two reflections they arrive back at the "starting point" in phase. For example, 
the cello string has as resonant modes all of those sine waves which go to 
zero at both ends. Thus the longest wavelength, resonant mode consists of 

a half cycle of a sine wave, but 1, l£, 2, 2%, cycles will also resonate. 

Thus each successively higher order resonant mode has an additional zero 
crossing along its length. Since the reflections at the end are indpendent 
of the wavelength (ideally), this cavity possesses an infinite number of 
resonant modes. In addition, for this simple system, the wavelength and the 
period of a wave motion are linearly related so that these resonances have 
a simple harmonic relationship between the periods. Such systems, for which 
waves of all periods propagate at the same speed arc said to be "dispcrsionlcss" 
a packet of energy consisting of a range of wavelengths docs not disperse as 
it propagates. As we have seen above, the propagation speed and reflection of 
acoustic waves in the solar atmosphere depend both upon the period and the 
horizontal wavelength. Hence they are dispersive. Thus while vertically 
propagating waves with periods of 300 and 400 seconds will be reflected down- 
ward just below the visible surface because the scale height gets too short, 

200 second waves will just propagate into the visible atmosphere and 100 
second waves will propagate right past the temperature minimum and on up 
into the chromosphere. So in fact, the cavity does not exist for these short 
period waves. Purely vertically propagating waves will go all the way in to 
the center of the sun without being reflected, so the cavity will be extremely 
thick for them, while waves propagating very obliquely will be reflected 
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back upwards before they have penetrated very deeply. 



At a fixed horizontal wavelength, the longest vertical wavelength 
which resonates has the longest period, and higher harmonics have shorter 
periods. Since the downward reflection no longer exists for sufficiently 
short period waves, only a finite number of resonant modes exist below the 
solar surface - in contrast to the situation for the cello string. In addition, 
while for the nice, uniform cello string all of the higher order resonant inodes 
had the wave energy distributed more or less uniformly along the string which, 
of course, has the same length for all of the harmonics, in the subphotosnheric 
cavity the lower order resonant modes propagate more nearly horizontally than 
the higher order ones and thus they propagate in a more limited range of 
depths in the solar convection zone than the higher order (shorter vertical 
wavelength, shorter period) resonant modes. Of course the resonance will 
exist for other horizontal wavelengths, for which - because of the dependence 
of the vertical wavelength on the horizontal wavelength and the period - 
the resonant period will be different. The ensemble of resonant modes with 
the same number of vertical wavelengths trapped in the cavity is often referred 
to e$ a 'Vertical mode" and the number of zero crossing of the pressure 
fluctuation - for example - used as an index. Hie longest vertical wavelength 
mode - which has no zero crossings - is called the "fundamental". For a given 
vertical mode, longer horizontal wavelengths penetrate deeper in the solar 
convection zone to higher temperatures; thus the cavity is mere extended, 
the vertical wavelength longer and hcncc the period longer. The precise 
functional dependence of the period upon the horizontal wavelength results 
from the detailed mathematical formulation, however the qualitative relationship 
can be arrived at quite simply. At any temperature in the atmosphere, the 
characteristic period * the acoustic cut-off period - is just proportional 
to the square root of the temperature, while the characteristic length - the 
scale height - is proportional to the temperature. Thus in an atmosphere 
with a more or less smoothly varying temperature, by keeping the ratio of 
the period to the characteristic period and horizontal wavelength to the 
characteristic length constant, that is since 
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f ~ / 1 


and 


X - T, 


then 


p ~ A, 

the resonant nodes with the sane number of vertical modes all "look the same", 
with the longer period, long horizontal wavelength ones just being shifted to 
higher temperatures - deeper in the convection zone. This square root dependence 
is, in fact, precisely the analytic result for an atmosphere consisting of 
a linear temperature increase, and it obtains as the large horizontal wavelength 
Unit for the aore ccmplicated temperature distribution in the solar sub- 
photospheric cavity. 

Figure 1 displays the low order resonant modes calculated analytically 
for a model solar atmosphere consisting of an isothermal temperature minimum 
of 4000°K, 1.5 Mm thick, between a linear rise of 100°/ km into the "corona" 

and a rise of 10° / km into the interior. (It is generally more convenient 

for theory and analysis to proceed using the Fourier transforms of the real 
variables, so reciprocal periods (frequencies) and reciprocal wavelengths 
(wavenumbers) are most commonly used.) Acoustic waves can propagate above 
the shaded band. The fundamental mode is evanescent everywhere, and it lies 

very close to the compressionless waves for which the divergence of the 

velocity vanishes. "Internal gravity" waves - corresponding to the buoyancy 
restoring force - exist in the lower shaded region, and the first three 
internal gravity wave modes arc shown. Note the square root dependence of 
the frequency on the horizontal wavenumber for small wavenumbers. The 
observed oscillations correspond to horizontal wavenumbers of the order of 
several tenths of a Mn , and the higher horizontal wavenumber behaviour - 
where the modal curves are open upward, rather than downward - which corresponds 
to a constant number of vertical wavelengths in the temperature miniraum- 
reglon-has not been observed. Hie S - shaped bends in the higher harmonics 
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occur as Che node crosses a line corresponding to Che same number of vertical 
wavelengths in Che isothermal temperature minimum as a lower order mode. The 
detailed vertical variation of the energy density in the resonant modes is 
unfortunately not so amenable to a simple qualitative discussion. But it 
la interesting to note that while for horizontal wavelengths greater than 4“ 
times the scale height at the temperature minimum the oscillatory energy 
of the fundamental mode is concentrated in the sub-photosphcric cavity, for 
nailer horizontal wavelengths the energy concentrates in the chromospheric 
cavity - as suggested by Bahng and Schwarzschild. Furthermore, the maximum 
amplitude of the higher order modes alternates from the convection zone to 
the chromosphere from mode to mode. 

Ifcimerical calculations with temperature and density distributions given 
by realistic solar models have been carried out by Ando and Osaki (1975, 

1977 ) and Ulrich and Rhodes (1977). Both groups draw attention to modes 
which are trapped in the chromosphere with a period near 240 seconds. As 
mentioned above, analytic studies predict a whole series of such chromospheric 
modes. Deubner has found some observational evidence suggesting the possible 
existence of power at this frequency. However, quite frankly, it uould be 
somewhat surprising were this mode to exist because of the large amplitude 
of the motion at chromospheric heights. That is, the wave will significantly 
distort the resonating cavity. The least that one might expect is that the 
modes with different horizontal wavelengths, instead of being completely 
independent, would have their phases locked together. As we shall see below, 
the currently discussed excitation mechanisms all operate below the visible 
surface and hence the excitation of the chromospheric mode is not expected 
to be favored. However, it should be pointed out that for comparable amplitudes 
near the visible solar surface, the chromospheric mode should be much more 
visible in the chromosphere than the normal modes with their maximum in the 
sub-photospheric cavity. While the astounding coincidence with observations 
served to establish the sub-photosphcric cavity model, more recently small 
remaining differences between observation and theory have been exploited as 
diagnostics of the thermal structure of the hydrogen convection zone. It 
may be that no plausible choice of a mixing length would be capable of 
precisely matching the observations, thereby vitiating the mixing length 
treatment of convection. The different depth distribution of the energy in 
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the different nodes appears to be capable of providing a diagnostic of the 
variation of the solar rotation rate with depth (Dcubncr, Ulrich and Rhodes , 
in preparation). 

While in what has been said up to now, any horizontal wavelength was 
as good as another; along a model curve, for every wavelength there corresponded 
a period. If , however, the oscillations last for a very "long" time, where 
now "long" is understood to be of the order of t’.e time for the phase to 
propagate around the circumference of the Sun; then even though there arc 
no reflections, the wave can interfere with itself after going all of the 
way around and once again only those inodes with wavelengths which end up in 
phase with themselves will resonate. The necessary observational resolution 
to resolve the horizontal modes is not impossible to imagine, and the existence 
or non-existence of this secondary modal structure should be an exciting 
diagnostic of the lifetime of the modes. In any case, there are several 
hundred thousand modes which contribute to the observed "five-minute" 
oscillation! 

While the fact that the observed structure of the "five minute" oscillations 
coincides so impressively with the normal modes of the sub-photospheric cavity 
gives us confidence in the appropriateness of the model, I haven’t given 
the reader any reason to suppose that there is any means of putting energy into 
these modes; or in the lingo of the trade, to "excite" them. (Just because 
unicorns could exist, one is not compelled to assume that they do.) As often 
seems to be the case, there exist a plethora of mechanisms available and the 
problem seems to be more choosing amongst them. The different mechanisms 
differ not only in their dependence upon the period and horizontal wavelength 
of the mode, but also in their dependence on the structure of the atmosphere. 

Thus for other stellar types the relative importance of the different mechanisms 
will almost certainly change. 

Two fundamentally different ways of putting energy into a wave, or 
oscillatory motion exist: internally or externally, via the energy equation 

or via the momentum equation. By pushing on the gas in a cavity- - mechanically 
through a diaphragm, as discussed above, or via some other motion in the 
atmosphere - at one of its resonant frequencies, the gas will be excited to an 
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oscillation at the amplitude specified by the pushing ("driving" is the more 
conmonly used lingo). The sound wave could exist at any amplitude in this 
formulation, the non-linearities having been neglected. Thus if a resonant 
mode was oscillating away at some amplitude and the driving was turned on at 
a lower amplitude, the driving would reduce the oscillation amplitude and thus 
would really be dissipation. Lighthill has formulated this approach in 
connection with the generation of sound waves by the turbulence in the exhaust 
from jet engines, and it has been used in astrophysics - primarily by Stein - 
to evaluate the flux of sound waves emitted by the "turbulent" flow in the 
hydrogen convection zone. "Turbulent flow" is the buzz word used to describe 
the velocities responsible for carrying the convective flux of energy - as 
specified by the mixing length treatment of convection. Given a characteristic 
velocity, which depends upon an assumed numerical and functional form for the 
mixing length, a turbulent spectrum of temporal and spatial scales is assumed. 
Unfortunately, the emitted flux is very sensitive to all of these assumptions. 
This flux is primarily at very short periods (less than 100 seconds) which 
are not trapped. Although these high frequency waves do not contribute to 
the "five-minute" oscillation, they may well be important in the heating of 
the lower solar chromosphere. Recently Coldreich and Keeley (1977) and 
Keeley (1977) have considered the same mechanism for the longer period, 
resonant modes and they find that although the turbulence acts primarily to 
dissipate the radial modes, it is capable of exciting the non-radial modes 
to amplitudes comparable with observation. This mechanism arises primarily 
in the upper parts of the convection zone where the convective velocities 
and the associated dynamical effects rise rapidly before falling to zero in 
the convective ly stable photosphere. While this pushing acts on the wave 
through sources which appear in the momentum equation, periodic addition of 
internal energy has also been suggested as a mean6 of exciting the resonant 
modes. 

As the temperature and density of the gas vary throughout the oscillation, 
the interaction of the gas with the mean atmosphere will also vary. Depending 
upon the phase between the energy transfer and the temperature fluctuation, the 
effect will increase or decrease the amplitude of the wave. For example, if 
heat is lost when the gas is hot and gained when the gas is cold, the oscillatior 
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will be damped. By reducing the energy density when it is high we reduce the 
restoring force of the comprcsibility and hence the energy of the wave has 
been reduced. If on the other hand, the gas could be made to absorb more energy 
when it was hot and less energy when it was cool, the amplitude of the wave 
would increase. This occurs quite generally for any thermodynamic cycle: 
heat absorbed during the hot part of the cycle adds to the energy of the 
oscillator, heat added during the cool part of the cycle reduces the energy. 

When the energy added is proportional to the amplitude - often called a "linear 
instability" - the amplitude will grow exponentially in time. An oscillator 
whose amplitude increases exponentially in time is said to be "overstablc" 
or "vibrationally unstable". The amplitude will continue to increase until 
some process whose importance depends non-linear ly on the amplitude grows to 
dominate - a parabola, or higher order polynomial, starts rising from zero 
more slowly than any straight line but eventually will exceed it. These non- 
linearities may occur in the excitation or in the dissipation process. If 
the amplitude of the wave in the driving region becomes large, over part of 
the cycle the excitation may change to dissipation. Thus any further increase 
in the amplitude would increase the dissipation more rapidly than it increased 
the excitation and thus the growth of the oscillation would cease. This 
limitation is often referred to as a "saturation" of the excitation. Leakage 
from the cavity may also limit the amplitude. The sound waves penetrate 
slightly through the temperature minimum region and will be dissipated at 
greater heights. This dissipation increases rapidly in its efficiency with 
increasing amplitude; so that were nothing else occur ing, the wave amplitude weui 
increase until the energy dissipated above the temperature minimum just equaled 
the energy put into the wave by the overstability. 

Such overstabilities have been invoked for years to account for the Ccphcid 
pulsation phenomenon. Although several effects are intimately combined, the 
increasing opacity with increasing temperature (at constant density) which 
prevails in the hydrogen and helium ionization zones of stellar envelopes 
takes energy from the outward thermal flux and transforms it to acoustic waves. 
Rather than controlling the deposition of energy into the wave, the opacity 
mechanism regulates the escape of energy through the wave; and, following 
Eddington, this distinction between periodic injection and periodic damming up 
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of a continuous flux is referred to as the "Eddington Value". This fluctuating 
opacity mechanism is often referred to as the "K mechanism" in distinction to 
the contribution made by the fluctuation in the internal energy brought about 
by the fluctuating ratio of specific heats - the "^mechanism". At temperatures 
above the maximum of ionization, further into the interior, the situation 
reverses and increasing temperature leads to a decreasing opacity, and thus 
fluctuations deep in the stellar interior arc damped. Since the driving 
"changes sign" at high temperatures, as the temperature fluctuation in an 
orcillation grows with amplitude, eventually the driving will saturate and the 
amplitude will be limited. Resonant modes with large amplitudes in the damping 
region will not be excited, while those with large amplitudes in the outer 
regions, where hydrogen and helium are ionizing, will be excited. Thus, along 
a modal curve for-say-the first harmonic, since longer horizontal wavelength 
modes will resonate deeper in the interior they will grow less rapidly in 
time than Sorter horizontal wavelength modes; and very long horizontal 
wavelength modes will most probably not be overstable. However, Provost (1973) 
has cautioned against such "intuitive" reasoning. Ando and Osaki (1975, 1977) 
have done extensive calculations of the variation of the growth rate with 
horizontal wavenumber for the different modes. 

/ 

While this overstability arises in the outer region of the solar convection, 
zone, near the depths where the "Lighthill mechanism" is effective, the thermal 
rather than dynamical coupling to the mean atmosphere should uncouple the 
relative importance of these two mechanisms in other stars. Ando (1976) 
has calculated growth rates of the overstability in other stars. 

Another physically distinct overslability, suggested long ago by Cowling, 
revived by Spiegel and discussed recently in detail by Graff (1976), derives 
energy from the super-adiaba-u: temperature gradient maintained by convectin. 
Because the mean atmospheric temperature gradient is greater than 
that of the displaced gas in the sound wave, any dissipative mechanism - such as 
radiative transfer - will add heat when the fluctuation is displaced towards 
the interior; that is f compressed and heated; and remove it during the other 
half cycle. (This type of mechanism has also been suggested to generate MUD 
waves, particularly in connection with the problem of cooling sunspots. 

A recent article by Cowling (1977) is of interest here.) It appears that for 
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the Sun, the K mechanism contributes much more than this "Cowling - Spiegel" 
to the overstability of acoustic modes in the sub*photosphcric cavity. 

Because of the tremendous efficiency of convection in the Sun, the actual 
temperature gradient exceeds the adiabatic gradient by less than a part in a 
thousand throughout almost the entire convection zone. It is only in the outer 
scale height or two that substantially super-adiabatic gradients exist, and 
thus only a relatively small mass is available to participate in the driving. 
However, in other stars where convection may be less efficient and extended 
super-adiabatic zones present, this mechanism may contribute. 


In concluding this section, let me repeat that the theoretical "action'' 
in the near future appears to focus on the generation of these resonant modes 
whose existence has been established. To evaluate the amplitudes attained 
by the overstabilities for comparison with observation and to evaluate the 
energy loss by dissipation, a much more difficult non-linear problem must be 
attacked. In addition, having established a new paradigm, new observations 
must continue to test it. In particular, the interaction of a magnetic field, 
be it in a sunspot, an active region, or in the quiet sun network, with the 
resonant modes is of great interest. The observations repotted on at this 
meeting by Livingston (see also Giovanelli, Livingston and Harvey, 1978) 
provide an excellent stimulus for this undertaking. 
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II. Chromospheric Dynamics 



One of Che principal programs of the pointed instruments onboard 0S0-8 
has been the observation of velocity and intensity fluctuation fields through* 
out the solar chromosphere and th n transition to the solar corona. These 
fluctuating fields are of interest not only as phenomena in their own right 
but also because of their relation to the mechanical energy transport which 
has been speculated tp be responsible for the heating of the chromosphere 
and corona. Because spatial inhomogeneities are known to be much greater 
in the chromosphere Chan in the photosphere and the effect of magnetic fields 
more important, we may anticipate that several processes may simultaneously 
contribute to the chromospheric dynamics and that the isolation and identi- 
fication of the phenomenon may prove more difficult. Thus humbled, it 
behooves us to examine the apparent consequences for the chromosphere of 
the model of resonant modes trapped in the suh-photosphcric cavity. We shall 
show that the ''requencies which cannot propagate near the temperature minimum 
and are trapped in the convection zone, propagate once again in the higher 
temperature chromosphere. In addition, waves at the cut-off period are generated 
near* the temperature minimum and these play an important role in the chromos- 
phere. Finally , even for resonant oscillations with amp’itudcs smaller than 
those observed near the visible surface, the chromospheric amplitudes are 
sufficiently large that non-linear effects become i-'Oortnnt. These effects 
will be illustrated by the results of numerical calculation- of the one dimension 
response of the chromosphere. 


The resonant modes trapped in the sub-photospheric cavity occur because 
waves of the frequencies of these modes cannot propagate at the temperature 
minimum where the scale height becomes too small. An evanescent, exponentially 
varying wave extends into the photosphere and the energy density decreases 
upwards towards the chromosphere. Although the energy density decreases, 
the velocity amplitude actually increases upwards, just not rapidly enough 
to overcome the effect of the rapidly decreasing density. Although the velocity 
becomes large, it is very difficult for these evanescent waves, with pressure 
and velocity out of phase, to dissipate, and they certainly will not form 
shock waves in the normal sense of the term. At large amplitude, they presumably 
couple to higher frequency waves which can propagate and dissipate. If t lie 
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lew temperature region extended to infinty, the reflection would be perfect. 
However, efter approximately seventeen pressure scale heights, the temperature 
rises and the resonant mode frequencies can once again propagate. If the 
amplitudes were very small; a resonant cavity similar to that in the hydrogen 
convection tone would exist - as discussed above - and the two cavities would 
lm coupled by the evansecent waves in between. As we shall see below, the 
waves are sufficiently strong to distort the high temperature reflection 
region of the chromospheric cavity and the chromospheric resonance is 
destroyed. Thus somewhat higher than 2 megaaeters above the visible surf a 
the five minute oscillation will start propagating once again, but the 
amplitude of these modes is now quite large and as soon as they can propagate, 
they start to dissipate their energy. 

Another process, whose existence night not have been anticipated frost the 
modal analysis, seems to be of greater importance however. Lamb - who I 
Indicated earlier discovered practically all of the basic physics of atmospheres 
over fity years ago - showed that when an isothermal atmosphere is subjected 
to a pulse of energy, a wake remains after the pulsc,oscillacing at the 
acoustic cut-off period. This is about 200 seconds for the solar temperature 
minimum region. Like the evansecent waves, this cut-off period wave transmits 
very little energy, and the pressure and velocity fluctuations are nearly 
ninety degrees out of phase. In contrast with the evanescent wave, the energy 
density does not decrease with altitude in the atmosphere. This wave appears 
to be excited by the 300 second evanescent wave in the solar photosphere, with 
which it enjoys a resonance - in that after two 300 second periods, three 
200 second periods have elapsed. This is illustrated in Figure 2 where the 
non-linear response of a plane parallel solar atmosphere to a pulse started 
1.6 megameters below the visible surface has been followed for 1400 seconds. 

A rigid lower boundary was raised about a kilometer in 30 seconds, so that the 
maximum upward velocity of the impulse was 0.05 km/ sec. The resulting 300 
second oscillation is rather weak only 0.1 km/scc at the surface. In contrast 
to the 300 second oscillation -at negative altitudes, above the temperature 
minimum (550 km) a very regular 200 second oscillation has been established. 
Above 1000 kilometers the velocity oscillations become significantly asymmetric 
in time, and the pressure pulses narrow so that by 2400 kilometers the pulse 
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is only about twenty percent of the period. High in the chromosphere, the 
pressure more than doubles in these short pulses, so that the atmosphere is 
blasted upwards and then nearly free falls back down until the descent is 
wrested by the next upward blast. As low down in the chromosphere as 1800 
km, the particle excursions are of the order of 600 km, that is to say several 
scale heights. Thus even for this relatively low amplitude sub-photospheric, 
resonant mode, there remains no semblance of a mean, equilibrium chromosphere. 
Another resonant mode, of infinitesimal amplitude, could no longer interact 
with an equilibrium atmosphere, but rather it would sec drastic time changes. 

Attempts have been made to understand the temporal evolution of strong 
chromospheric lines such as the calcium K resonance line (Keasley (1975), 

Cram (1976)) and the hydrogen Lyman or line (Kneer and Kakagawa, 1976) 
by calculating line profiles from the chrc .osphere as a pulse, like the one 
which initiates the motions in Figure 2, propagates upwards. Durrant, 
Grossman-Doerth and Kneer (1976) have argued that such an isolated pulse 
with pressure and velocity in phase, is not capable of reproducing the 
observed profile changes. In fact, one sees that the situation following the 
initial pulse is distinctly different, in that -.he pressure and velocity 
fluctuations are significantly out of phase - nearly ninety degrees. Thus 
in a real, dynamical atmosphere the 200 second pulses arc always moving into 
an at me. phere f*r removed from equilibrium by the passage of the previous 
pulse. Of course, there may exist other sources of pulses with no special 
phase with respect to the 200 second oscillation - those generated by exploding 
granules for example. But these too will move into an atmosphere whose tine 
averaged properties bear little rescmblcncc to the equilibrium, static 
atmosphere - there is no way that a micro-turbulent pressure is going to model 
even the mean state of an atmosphere with motions greater than the scale height. 
Furthermore, because of the substantial motions of the atmosphere, the velocity 
at constant mass depth - that is, measured in a Lagrangian coordinate - which 
bears a closer relation to what a spectral line will "see" than the velocity 
at a constant geometrical depth, contains striking differences for the upper 
chromosphere. Figure 3 displays the Lagrangian fluctuations for the same 
calculation presented in Figure 2. The phase velocities are reduced high in 
the chromosphere because of the outward impulses. The prises also tend to form 
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pairs, which is seen most dearly in ‘ it pressure fluctuation - the first pulse 
•lews down the existing infall and Liu .econd rc-accelcratcs the material 
outward. The resulting outward motion moves this mass depth away from the 
next outwardly moving pressure pulse, which as a result takes an increased 
length of time to catch up. 

Although the most striking aspect of chromospheric dynamics - the decrease 
of the period of oscillation from 300 to 200 seconds - is easily accounted 
for within our present understanding of sub-photospheric resonant nodes and 
propagation in the photosphere and chromosphere, the detailed evolution of 
optically thick spectral lines, which (presumably) provide us with much 
improved height resolution and a better constrained diagnostic, has not been 
satisfactorily described. While those who expected 0S0-S to observe beautiful 
M wave shock profiles to confirm the shock wave heating hypothesis (and deliver 
them their science on a silver platter) must feel somewht disappointed, it 
certainly appears at this point that the ball has merely been placed back 
in the theoretician's court: it is incumbent on "us" to provide a meaningful 

expectation of the observable consequences of the "astrophysical truth" of 
the photosphere. 


III. Chromospheric and Coronal Keating 

The lack of unambiguous success in identifying propagating and dissipating 
shock waves in the chromosphere has certainly impeded the solution of the problem 
of chromospheric and coronal heating; but in addition to the diagnostic problem 
cited above, other theoretical aspects of a self-contained shock wave heating 
model remain. to be resolved: total flux generated, distribution of the flux 

with period and horizontal wavelength, height distribution of the dissipation 
a gai n as functions of the period and horizontal wavelength, effect of changes 
in the mean atmosphere on the generation and dissipation; to list but a few. 

It has been argued that the Sun is the ideal stimulus for attacking these 
questions not only because of the spatial and temporal resolution which its 
propinquity affords, but also because of the relative weakness of the enromos- 
phere/corona phenomenon. Because only a very small fraction of the total flux 
appears to be required to heat the outer atmosphere above its radiative 
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equilibrium temperature, only very .small modifications to the conditions in 
tiie convection cone - or wherever else the generation of e propagating mechanical 
energy flux from a thermal flux may occur - would seem necessary. That is, the 
Structure in the generation region may be calculated neglecting the mechanical 
flux, which is then just calculated as a perturbation - it exerts no effect 
hack on the local, thermally controlled structure, such as might exist if 
very efficient wave generation refrigerated the atmosphere. However, because 
of the weakness of the pheonotnenon it nay well turn out that no one effect 
dominates and we may advance our understanding by studying the outer atmospheres 
of stars with more significant dynamical effects: ccphcids or stars with more 

Vigorous convection zones. Currently the situation is embarrassingly reversed; 
stellar astronomers seem to be taking unfounded solar hypothesis seriously and 
ere applying them out of context. (See Linsky and Ayres (1977) for a critique). 

These hypothesis have not changed significantly in recent years. Let 
me just outline them for the sake of completeness. Short period (less than 
100 seconds) non-resonant waves generated by turbulence in the convection zone 
dissipate their energy by non-linear processes after propagating a shorter 
distance than the longer period, nearly resonant waves essentially because 
the short period waves have shorter wavelengths. Quite sinply, they don't 
have to go so far in order to break. In addition to the horrible ur.ccrtainty 
in the power generated by this mechanism mentioned earlier, the short wave- 
lengths - comparable to the region over which a spectral line is formed - 
have posed problems to the direct observation of these waves. Dcubner (1977) 
has recently interpreted wavincss in the temporal power spectra of photospheric 
oscillation at these periods as an indication of the existence of these short 
wavelengths. Ulmschncidcr and Kalkofen (1977) argue strongly in favor 
of the importance of this short wavelength heating, while Proderic and Thomas 
(1976) contest the treatrmnt of radiative losses, Jordon (1976) supports 
Ulmschneidcr's position and Cram (1977) questions the whole edifice. These 
short wavelength waves would be primarily of importance in heating the low 
chromosphere, where they dissipate the remaining small fraction of their 
energy not lost by radiative damping as they propagate through the photosphere. 
Essentially none of their energy remains to heat the corona. 
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The Model of heating by dissipation of sub-photosfheric resonant nodes 
in the corona has also suffered the difficulty of fundamental uncertainties 
In the energy generation rate. However, this appears to be in the process 
of being resolved. These nodes have the advantage of being directly measured 
in the photosphere so that an observational normalisation is available. 

The direct calculation of the amplitude achieved by the over stability of non- 
vertically propagating modes is a challenging problem, which unfortunately 
has not been challenged. Since the amplitudes in the driving region are 
very small (less than one thousandth the sound speed), the oscillations will 
find it difficult to saturate the excitation mechanism. Calculation of the 
limitation of the amplitude growth by shock wave dissipation requires the 
coupled treatment of the overstability and of two dimensional shocks. While 
progress may be made using existing one dimensional shock wave computer 
programs, the determination of the importance of two dimensional effects must 
await development of the appropriate programs. According to Bob Stein, 
several groups are attacking this problem via different techniques. 

Several authors have calculated the thermal structure to be expected 
from a balance between shock wave heating, radiative losses and thermal 
conduction using analytical expressions for the dissipation of the mechanical 
flux. Because the wavelengths are long compared to the scale height, approximate 
treatment of the propagation and dissipation of these waves by methods which 
assume the medium to be uniform are not valid. Nonetheless, a great deal 
of work and even more publications are based upon this invalid treatment. 

Stein and Schwartz (1972) have demonstrated the inappropriateness of applications 
of the "weak shock theory", which makes additional assumptions essentially 
equivalent to assuming that the shock doesn't exist, by comparing results with 
the full non-linear treatment. 

In closing, it is appropriate to point out that while there remain many 
theoretical aspects of the hypothesis of shock wave heating of the outer solar 
atmosphere which require development before the hypothesis can be accepted 
or rejected, alternative hypotheses exist which may in time prove to be more 
successful in describing stellar chromospheres and coronae. In particular, 
very strong (greater than 1000 gauss), very small (less than an arc-second) 


22 



MgMtic fields are known to exist in tie super 'granulation network boundaries 
All across the “quiet sun“, and these will support magnetic um /nodes with 
generation, propagation and dissipation prop**rites quite distinct fren the 
purely acoustic aodes described here, lbe work of Defouw ( 1976) , Piddingtor. 
(1973), Dchida end Kaburaki (1974), tfcatxel (1974) 

should be consulted. lr. addition , one can imagine thrt these fields arc 
capable of dissipating their energy ii either eery small f iarc-li kc , 
"catastrophic" events, or more uniformly it*, cioe - to give pla -,-s. physicists 
an easy problem for a change. In any case, I hope that i»; 1.-; clear cl.-tr 
theoreticians have before than a vide range of e.:cit*.ng, r*l»ri-'*iy 1 1 
problems and they may take heart from the recent success in the pre-ilrci-'.i of 
the nodal structure of the five minute oscillation - some problems actciiiy 
do manage to get themselves solved! 

It Is a pleasure to thank Robert Stein and Pierre Gouttebroze for many 
useful discussions and Lise Ouvrard for correcting the rramar. The non- 
linear calculations were performed on the CDC 7600 of the Centre Nations! de 
Etudes Spatiales. This work was supported in part by NASA contracts NASW-3053 
and KAS8-32356, and the Lockheed Independent Research Program. 
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Figure I: Resonant modes for e plane parallel atmosphere with a 1.5 megameter 

thick temperature minimum at 4000°K, a 100°/km temperature gradient above and 
• 10°/km gradient below. Limits for propagating acoustic and internal gravity 
vaves in the temperature minimum region are shaded. The fundamental and 
first three acoustic and .internal gravity resonant modes are shown. 

Figure 2: Velocity (solid lines) and relative pressure (dotted lines) 

fluctuations as a function of time at various altitudes in a model solar 
atmosphere. Altitude (in kilometers) is measured from a zero at '^qqq * 1 in 
the initial atmosphere. The energy equation is adiabatic with a constant 
V Y* 1.4. The velocity and pressure normalisation at each altitude was chosen 
to exhibit the fluctuations after the a tmosphere had recovered from the 
initial pulse, and points out of range have not been plotted. 

Figure 3: Same as Figure 2, except the fluctuations are measured at constant . 

mass depths, corresponding (approximately) to the altitudes displayed in 
Figure 2 for the initial atmosphere. 
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We have observed simultaneously the solar CaK (3933 A), 

CaH (3969 A), Mgk (2796 A), Mgh (2803 A), and H-Lyo 1216 A 
solar lines with spatial, temporaj and spectral resolutions 
of 1" x 10", 10 seconds and 0.02 A, over time sequences ranging 
from 40 to 60 minutes. 


When plotting as a function of time the integrated intensity 
of the Mgk line (over 1.09 A), a clear quas i-pc r i od i c oscillation 
is visible on 48 out of 50 time-sequences ; 16 of the 48 measured 
periods are shorter than 197 seconds, 16 between 197 and 238 
seconds, 16 greater than 238 seconds. 

Within the 10 seconds resolution, the temporal evolutions 
of the intensity of the blue emission peak of the CaK and Mgk 
lines are correlated without phase shift. 

A study of thepower spectrum of the temporal evolution 
of the intensity of the blue emission peak of the CaK and Mgk 
line', for 12 time sequences, indicates that strong, short 
period oscillations (< 200 seconds) occur over quiet areas of 
the sun, such as the center of the supcrgranulation cells, 
whereas weak, long period (>250 seconds) oscillations occur 
over areas with active CaK profiles. 

For 7 time sequences we have computed the average Lyman-a 
profile corresponding to the "blue" and "red" phases of 
the solar oscillations (i.e., for a time sequence of 300 spectra, 
we add the 100 Lyo profiles measured simultaneously with the 
100 Mgk profiles with the strongest blue peak, and repeat one 
other Lyo average corresponding to the 100 Mgk profiles with 
the weakest blue peak). The resulting "blue" and "red" 

Lyman-o profiles have the same total intensity, the same 
intensity at the emission peaks, whereas the corresponding 
CaK and Mgk profiles are strongly different. The central 
reversal of the average "blue" Lyo profile is shifted to the 
red with respect to the average "red" profile, indicating 
that the part oC the solar atmosphere where this central reversal 
comes from oscillates somehow in phase with t he lower ehromos phe rc 
where the' Mgk lines originates. 
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ABSTRACT 

The in-flight performance for the first 18 months of operation of the French, pointed instru- 
ment on board OSO 8 are described. The angular and spectral resolution, the scattered light 
level, and various other instrumental parameters are evaluated from the observed data and 
shown to correspond mostlv to nominal design values. The properties of the instrument are 
discussed, together with their evolution with time. The distribution of the first 8363 orbits between 
various observing programs is given. Preliminary results are also described. They include studies 
of the chromospheric network, sunspots and active regions, pron* nonces, oscillations in the 
chromosphere, chromosphere-corona transition lines, and aeronomy. 

Subject headings: instruments — Sun: chromosph re — Sun: prominences — Sun: spectra — 
Sun: sunspots — ultraviolet: spectra 


!. INTRODUCTION 

One of the major goals of solar physics is to under- 
stand the nature, origin, and evolution of the various 
features present in the solar atmosphere. Some, like 
the granulation, represent dynamical responses to the 
convection zone. Others, like spots or the more dis- 
persed magnetic flux tubes (e.g., network fragments), 
represent symptoms of a magnetic process. Anv 
advance in cur understanding of such interior pro- 
cesses must come from high angular and spectral 
resolution observations of the line profiles of such 
features. 

The NASA orbiting solar satellite OSO 8 , launched 
on 1975 June 21, carried in its pointed section two 
instruments designed for the highest angular and 
spectral resolution achieved by spacecraft to date. One 
of these instruments w*as the responsibility of the 
Laboratory for Atmospheric and Space Physics 
(LASP) of the University of Colorado, the other was 
that of the Laboratoirc de Physique Stellaire et 
Plandtaire (LPSP) of the Centre National de la 
Recherche Scientifique (France). 

In this paper, we describe the performance achieved 
in orbit and outline the main results obtained with tne 
LPSP instrument after 18 months of successful 
operation. Most of these results arc in a preliminary 
state. A complete description of the instrumentation 
has been given in Aruncr et ai (1977), hereafter 
referred to as Paper I. 


• Lockheed Palo Alto Research Laboratory, Palo Alto. CA. 
t Hieh Altitude Observatory, National Center for Atmos- 
pheric Research; and University of Colorado, Boulder The 
NCAR is spoitsoicd by the NSI , 


H. SUMMARY OF THE INSTRUMENT CAPABILITIES 

Because of the limitations imposed by the size of the 
spacecraft (although considerably larger than the 
previous OSO*). the LPSP tdc cope was a Cassegrain* 
ian with a dtametc r - . »u, cm. Consequent I> we limited 
our observation^ to e most interne chromospheric 
lines. Because the chromosphere is of limited depth, 
the spectrometer was designed so as to simultaneously 
observe six lines: 

Ca ii H (396.9 nm) and * 393.4 nm) ; 

Mg n h (280.3 nm) and k (279.6 nm) ; 

H i Lo(l21.6nm)and L0(1O\5 nm' 

A very rapid and versatile spectral scanner made it 
possibles also study the lines of 0 vi (103.2 nm) and 
Si in (120.6 nm) nearly simultaneously with those 
listed above and enabled us to study propagation 
effects and to obtain height resolution from the upper- 
photosphere to the lower corona. Table 1 summarizes 
the main characteristics of the spectrometer, which 
operates with two different spectral resolutions. 

Two metf were available to make specuohelio- 
grams. One was by means of spacecraft rasters: two 
image sizes were available, 44' x 40' and 2!75 x 213 
(nominal). The rcauer is .cferred to Paper 1 (Table 2) 
for more details on these rasters. In addition, the LPSP 
telescope had an articulated secondary mirror which 
was moved by a »wo-axis stepping mechanism. 
Accordingly the so’ ige could be moved step by 
step. Each step w a . „ the solar surface while the 

maximum area covered was 64* x 64". 

A slit wheel mechanism, at th: focus of* the tele- 
scope, was used to select various sli: sizes ranging 
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TABLE I 

Dnc«« OuftArrmspcs or im CSO $ LPS* Stx«OMSNa HiGH>RyiQumow Su ctm m i t ii 


Su ctui toumoN (nominal) 

High Mode SKCTtAi Inciument 


CbimiLwr 

(ami 

Low Mode 

lam) 

(nm) 

Lms-* 

Maximum Stcctrai Range 
( nm) 

HEX CXATINC Sl* P 

(nominal) (nm) 

Can II.. 

30* o 

O.l 

0.0020 

1.51 

305.IM 30K.244 

0.1X1140 

Can K 

vu 

0.02 

0.0020 

1.52 

301 444- 304. 042 

U 00169 

Ma«a.. 


0,1 

00025 

2.67 

277.739-282.264 

0 00210 

Mgn* 

4/9.6 

0.02 

00025 

2.6$ 

277.005-281 S$\ 

000220 

H l La 

121.6 

0.02 

0.0020 

4.03 

120.569-1 22.291 

0 00085 

Hi tg 

1025 

0.1 

0.0060 

17.65 

101 .686-103.229 

0 000^4 


Ken.— Wavelength units are nanometers (ntn). 


from P x I* toT x 40* as well as ~ * corresponding 
to 6* x r. 

m. IN-njCHT PERFORMANCE 

Prior to launch, the instrument was submitted to 
numerous tests and calibrations, the results of which 
are given in Paper I. Here we present the actual per- 
formance measured in flight. 


a) Angular Resolution 

The ground tests of the telescope led us to expect 
an instrumental profHc with a full width at half- 
maximum (FWHM) of 2*. Two methods have been 
used to estimate the angular resolution in orbit : 

i) f mopes of Limb Shape 

Repeated scans of the solar limb, as observed with a 
I* x 1* aperture using the internal raster mode in the 




Fm. I. M Internal raster scans utross the lunar limb in Ca K. The dots show the tcsults one would obtain assuming the tele- 
scope to be diffraction limited. The triangles show the results one obtains assuming the instrumental profile to be that represented 
on (6), the telescope instrumental prolile. 
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far wings of the Ca n line, show a limb darkening with 
a 507. decrease over an angular distance of 2*. 

However, this method does not separate the resolu- 
tion of the telescope from the jitter of the pointing 
system. Rapid, one-dimensional scan c the solar limb 
as well as other methods allowed us to estimate the 
amplitude of the rms jitter as 0T5. 

ii) Partial Solar Eclipses 

On 1976 April 29 and October 23, two partial 
eclipses of the Sun were visible from 0S08. Figure la 
represents spatial scans laterally across the lunar 
limb measured in relative intensity units. The actual 
observations are compared with the results of a com- 
putation assuming the 16 cm Cassegrainian telescope 
to be diffraction limited, and to have an instrumental 
profile as given in Figure 16. The fit using the profile 
of Figure 16 represents a good approximation. We 
therefore conclude that the instrumental profile has a 
FWHM of 275 ± 0T5. 

6) Spectral Resolution 
i) La and LB Channels 

For these we take advantage of the narrow absorp- 
tion line due to geocorona! hydrogen. By applying a 
method developed for the study of interstellar absorp- 
tion lines (Vidal-Madjar et al. 1977), we obtain a 
spectral resolution in flight of 0.002 ± 0.0005 nm at 
La, and 0.006 z 0.001 nm at L3. 

ii) Calcium and Magnesium Channels 

We compare the solar spectra obtained by our 
instrument with ground based (Ca n channels) and 
balloon- or rocket-borne observations ( Mg n channels). 

In Figure 2 we show* the full range Ca it K and 
Mg it k spectra from OSO 8 . The variation with 
wavelength of the sensitivity of the instrument has 
been corrected for, and the spectra are deconvolved 
from the instrumental profile. The comparison with 
the spectrum of the Kitt Peak Preliminary Solar Atlas 
(Brault and Testerman 1972), which .has a spectral 
resolution of 0.0024 nm, shows that our resolution in 
orbit is better than this value. In the case of the Mg n 
channels, comparison with the spectra of Lemaire and 
Skumanich (1973) and Kohl and Parkinson (1976) 
yields a resolution of 0.0025 ± 0.00025 nm. This value 
is equal to the nominal design value (cf. Table 1). 

c) Dispersion and Grating Mechanism 
(Spectral Scanner) Stability 

The dispersion law of the spectrometer was deter- 
mined in orbit by measuring the position, in units of a 
grating step, of II solar absorption lines of known 
wavelengths in the Ca n and Mg ii channels. For the 
LjS channel we used the O i lines at 130.48 and 
130.6 nm and N i at 119.9 nm (which appear in the 
1 1th and 12th orders of diffraction). 

Because of the lack of lines in the La channel, wc 
deduce the dispersion law from that at Lft. The absorp- 
tion line of geocoronal hydrogen provides an absolute 


reference. This proved to be valuable due to the 
appears nee of positioning uncertainties ( ± l grating 
step) in the movable La, L0 exit slit mechanism. The 
dispersion law was measured repeatedly to check for 
long-term variations. Over I year we found that the 
correspondence between absolute wavelength and 
grating step number varied by no more than ± 3 grat- 
ing steps (cf. last column of table 1). 

To check the mechanism stability over one orbit, 
we measured the position of the photospheric line 
391.52 nm in the wings of Ca ii K. Any departure 
from the orbital Doppler effect could be attributed to 
photospheric Doppler shifts and or changes in the 
spectrometer. The result is shown on Figure 3. One 
can easily recognize the 300 s photospheric Doppler 
oscillations after removal of the orbital Doppler shift, 
measured for the first lime from space. The amplitude 
of the residual noise on this curve amounts to 
± 30 m s“ L The stability of the mechanism over a full 
orbit day is better than one grating step and exceeds 
our design expectations. We are able to easily and 
accurately measure Doppler shifts of photospheric and 
chromospheric lines (see § Vc below). 

d) Scattered Light Background and Dark Current 

The level of scattered light in the telescope plus 
spectrometer was determined from partial eclipses of 
the Sun. From Figure 1 we see that at 1 V from the 
lunar limb this level amounts to 2% of the intensity 
of the disk in the calcium channels. For Mg ti it is 
4% . For La and L£ these figures become 10% and 
20%, respectively, which indicates that the scattered 
light level may vary with wavelength, roughly as 1 X 1 . 

No simple and unambiguous method was available 
to measure separately and give an absolute value for 
the amount of scattered light in the spectrometer due 
to wavelengths well away and near the wavelength of 
interest. 

In the case of the Ca it and Mg n channels, we could 
compare the performance of our spectrometer with 
those of other ground-based or rocket-borne instru- 
ments. The result of this comparison appears in Table 
2, where we give the ratio of intensities at Ca n and 
Kj, Mg: ii /ij and k u relative to those of the Ca n and 
Mg h line wings. We notice that our performance is 
excellent for the Ca n channels, for which these ratios 
are smaller than those deduced from the Utrecht 
(Minnaert, Mulders, and Houtgast 1940) and the Air 
Force (Beckers, Bridges, and Gilliam 1976) atlases. 
We also compare our values with those of Linsky 
(1970) and of White and Suemoto (1968) who used 
particularly good optical systems. The Utrecht Atlas 
was used to evaluate the ratio of the H 3 and K 3 
intensities relative to that of the continuum at 400.0 nm. 
The results are: 

Jjh. = 0.053 , « 0.065 . for OSO 8 . 

Mon o • .no o 

while White and Suemoto (I96S) find 0.071 ± 0.0015 
and 0.061 L 0 001, respectively, and Linsky (1970) 



0S08 LPSP 



Mg H k CHANNEL 



CA II K CHANNEL 

Fk; 2 Full range spectral scans in the Ca n K and Mg n A channels using a 1* * V entrance slit The ordinates are counts per 
gate Instrumental sensitivity variations over the wavelength range have been corrected Tor Notice the reversed asymmetries m tne 
Mg it h and k lines It takes approximately .Ml s to go from Mg u k to //. and the asymmetries reflect time variations in the line prohles 
o%v*i 5o s 
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TABLE 2 

Ratio of Intensities at th. "'ositions K u H t , k lm and h t to Intensities Measured 
in tv* Wings of c- a and Mg u Lines in the OSO 8 Channels 


Channel 

Ratio 

OSO s 

Utrecht Athx Air Force A tbs 

(Mtnnacri •** at 1940) (Deckers t\* at l9?eo 

and 

Parkm>on 

Ca it H 

/<J95 IS» 7(H ; > 

0084 

0085 




nm.m 7(H t > 

0 077 

0090 



CanK 

/<J9l.4?)7(Kt) 

0.078 

0088 

0 087 



/(?94.85>7(K t > 

0083 

0096 

0095 


Mg it h 

M»r) 7(2-7 7?) 

0 079 



0048 


/Uir)7<2F2 OH 

0063 



0O45 


/(Aii-V/C??. 73) 

0046 



0 028* 


/<A t r) /<282 01) 

0037 



C 02?+ 

Mg n k 

Hk t r ) /C" *3l 

0 048 



0 04S 


/<A t r) A 277 73» 

0.03 



0 028+ 


N< te. — The results are compared with the same ratios evaluated from Solar Atlases or available published data r and r referred 
to the red and the blue pan of the lines. 

• We are indebted to Dr J Koh! for providing original records of hts spectra 

* These values are probably uncertain due to the difficulty of measuring the solar intensity at k l% . 


0.0409 ± 0.0022 and 0.0434 ± 0.001 1 for the same 
ratios. Our results are intermediate between these two, 
confirming the very good performance of the Ca it 
channels. 

For Mg n we have used the spectrum of Kohl and 
Parkinson for comparison. Excellent agreement is 
obtained in the k channel, our spectra indicating 
nearly exactly the same amount of scattered light as 
in the comparison spectrum. The agreement is. how- 
ever, poor in the case of the h channel. Both the h and k 
lines can be observed in this channel together with the 
reference wavelength at 277.73 nm, allowing direct 
comparison between the two channels. As a result, we 
notice that the h channel has a much higher level of 
scattered light. Tus might be the result of degraded 
spectra! resolution, due to a defective adjustment of 
the common coma and astigmatism corrector used in 
the Mg it channels whose delicate adjustment was 
optimized for the k channel. 

For the La and L£ channels we have made computa- 
tions using the instrumental profiles, determined by 
ray tracing techniques, and the expected properties of 
the baffling inside the instrument ; and we find that the 
level of scattered light is approximately 3% of the 
maximum flux in both channels. 

Dark-current measurements were performed syste- 
matically each orbit during the first year and every 2 
or 3 days in Ihc second year. The dark current was 
found lo he stable* with nearly no change during 
IK months. The values arc respectively 1.3, l.l, and 
0.1 counts s' 1 for Mg »t, La, and L£. respectively. 

The calcium channel dark current was typically 
200 counts s *. This comparatively had performance 
is due to a leak in the enclosure system (skm). 

e) Photometric Standardization 

Photometric sensitivities have been measured 
regularly in orbit relatively to their values at launch. 
Such measurements aic made cveiy day or two mth 


the 1' x 10' entrance slit at disk center quiet Sun, in 
the high spectra! resolution mode for Ca it and Vlg u 
and low resolution modcfo r La and Lf. The number of 
counts at certain standard wavelengths 'see below) 
measures the relative efficiency, whose variation as a 
function of time is shown m Figure 4. 

The interpretation of these curves may be of interest 
to those who plan to utilize similar instruments in 
space. The telescope mirrors* the collimator, the 
grating, and all surfaces in the La, L3. and Me u 
channels were coated at the Goddard Space Flight 
Center, with Al- LiF (Bradford etal. 1969). Elaborate 
precautions were taken in the storing and handling of 
optics throughout the mounting and calibration of the 
instrument. In fact, a special 300 square meter facility 
was built with air cleanliness and with temperature 
carefully controlled and humidity always kept below 
30% (Saivetat 1975). A loss of sensitivity such as the 
one reported here is very unlikely due to a contamina- 
tion in the instrument "before the launch and should 
rather be regarded as caused by the outgassmg of the 
spacecraft and the instrument once placed in the space 
vacuum. In that case, the greater the number of reflec- 
tions, the larger the degradation. The presence of steps 
which appear at nearly the same time on all the curves 
of Figure 4 is probably the signature of sudden out- 
gassing periods. The significant differences noticeable 
between the indi vidua* curves, however, arc indicative 
of causes of degradation proper to each channel, 
a flee ting either the mirrors, the tillers, or the detectors. 

The two Lyman channels show nearly the same 
behavior, with the larger loss at La attributed to the 
larger number of reflections in this channel (>cven at 
La versus live at Lti). At day 540, i.c.. IN months 
after launch, the sensitivity at La and Lr? was 10 J 
and 5.10 ‘ J , respectively, of the value at launch. 
Assuming that each reflection is affected equally by 
the contamination (which is certainly a crude approxi- 
mation), these numbers indicate that each reflection 
has reached 35% .it L*t and *7% at l.**)' ol its value at 
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day zero. Such results are not particularly dramatic 
when compared with the photometric behavior of 
other solar space instruments (Huber et oi 1973). 

However, the combined effect of outgassing and the 
baking of the secondary mirror surface by a flux of 
more than IS “solar constants** is likely to be respon- 
sible in a large propot non lor the \eiiMtivitv loss. A 
simulation made a lew weeks prior to launch at NASA 
on AI + LiF coated samples illuminated by 17 “solar 
constants" and placed in a normally outgassing 
environment showed a decay in efficiency from 63% 
to 58% and from 67% to 35% at La and L 3, respec- 
tively, only 52 hours after pumpdown. 

Obvious solutions, such as closing a shutter in front 
of the telescope during the first orbits when outgassing 
is high, were unfortunately not possible and would 
have delayed the launch several months. 

Assuming, arbitrarily, that the secondary mirror is 
responsible for a loss of a factor 10 ar both La and LS, 
each surface would have reached an efficiency of 46% 
of its value at launch, which is more or less normal. 

Totally unexpected and more striking is the behavior 
of the Mg and Ca channels. Although of yet unknown 
origin, outgassing might also be responsible for the 
degradation observed in these channels, at least until 
day 160 when the sensitivity reaches 1 40 and 1 20. 
respectively, of the value at launch. This corresponds 
to an average loss per reflection of 48% and 55%. After 
day 160 the sensitivity in the *wo Ca channels rises 
again. At the same time a faster decay is observed in 
the Mg channels. This peculiar behavior is attributed 
to interference phenomena, probably complementary, 
in Ca and Mg within thin films of contaminants) 
deposited on any one of the optical surfaces. Deteriora- 
tions of the interference filters which are used in all 
these channels may also contribute. In the case of the 
Mg channel it is also very likely that the detector itself 
is responsible for the loss of sensitivity. This is 
apparently not the case for the two Ca channels since 
their sensitivity follows nearly the same variation with 
time, which more likely reflects a variation in the 
optics used in common. 

The overall loss of sensitivity compromised certain 
aspects of the observing program. However, the 
versatility of the instrument made it possible to obtain 
scientific data of high quality and value throughout the 
mission. 


/) Absolute Calibration 

The calcium channels were calibrated by comparison 
with the data of Linsky (1970) and Livingston and 
White (1978) which represent average quiet Sun 
condition*. The absolute intensity at our standard 
wavelengths II lm and were taken as 0.0751 and 
0.06S7, respectively, in units of the continuum intensity 
at 400 nm. 

For the magnesium channels we attempted to 
improve Bonnet's 1967 results (Bonnet 1968) and 
designed a high spectral resolution instrument cali- 
brated against a blackbody constructed by R. Pevturaux 
at the Institut d'Astrophysique de Pans. This instru- 
ment was launched twice on the LASP rockets number 
21029 on 1975 July 2S,and 2 1030 on 1976 February IS, 
but because of malfunctions in the electronics it did 
not give reliable results. We prefer therefore to rely 
on other recent measurements — e.g., those of Koh! 
and Parkinson (1976). The absolute intensity at the 
standard wavelengths and wcre taken as 8 and 
6 x 10’ 12 ergs cm* 2 $* : sr" 1 cm* 1 , respectively. 

For the Lyman channels we also used the above 
rocket program to carry packages consisting of j m 
Ebert-Fastie spectrometers, measuring the integrated 
solar disk : calibrate the LASP and LPSP instruments 
separately. Only the second flight yielded good calibra- 
tion data. The results are given in Table 3 and are 
compared there with other measurements. The LPSP 
values are somewhat high; however, they are in the 
d section suggested by geophysicists (Levasseur et at . 
!*„ . To use these integrated intensities in La and LB , 
we Uww jiet Sun average profiles computed for the 
whole disk (Fig. 5). We have taken into account the 
center-to-Iimb variation in an approximate way that 
will ultimately be improved upon by means of entire 
Sun raster-generated profiles (i.e., profiles constructed 
from spectroheliogramo), when these become available 
from the data tapes. The absolute flux at the standard 
wavelength (core of the line! used to monitor the La 
relative sensitivity was taken to be 3 x I0 10 photons 
cm’ 2 s“ 1 nm" 1 . 

IV. REAL TIME OPERATION AND PROBLEMS 

The various modes of operation of the instrument 
have been described in Paper I. Here we discuss the 
“real-time" operation mode which allowed one, for 


TABLE 3 

Integrated Solar Flux Measurements Obtained with Calibration Rockets a"* La and L0 


Calibration Rockets 


CaI IIIKATI t> r LUX 

1975 July 2R 
LASI** 

1976 Feb IS 
LI’SI* 

197ft Feb 18 
LASI- 

OSO 5t 
1975 Aug 8 

aec: 

19 , 5 Apr. II 

f(l0.7cm)(IO 1 11/ ’) 

75.5 

70 1 

70 I 

70 


/■(Lo) (crg% cm a $ *) 

4 02 

5 46 * 207, 

4 05 1 :ot; 

4.25 

none 

L0) (crg>cm* 2 $“ *) 

none 

0.078 ± 207. 

none 

none 

0.050 


• Rotiman 1977. 
t Vidal* Madia r 1977. 
t llmtctcggirr 1970. 
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Fig. 5 — Full Son average La and IS profiles w avelengths ii 
values of the last column of Table 1. 
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io the left The conversion into A units can be made using the 


the first time, to point from an orbiting observatory 
with an absolute accuracy of nearly l \ The LPSP and 
LASP instruments were both operated from LASP 
(Boulder, Colorado) with “resident” LPSP scientists 
and g"est investigators involved in the daily opera- 
tions. iarget selection alternated daily between LPSP 
and LASP until 1976 April and then week!). Fora more 
complete description of the operations command 
generation and quick-look facilities, see Jouchoux and 
Hansen (1978). 

a) Pointing System Problems 

The pointing system of the “sail” secticr of OSO 8 
uses either one of two Sun sensors (SEAS), designed 
by Hughes Aircraft Company, mounted on each instru- 
ment. These devices coalign the SEAS pointing axis 
with the optical axts of the associated telescope. This 
corrects for drifts of the optical axis with respect to 
the mechanical structure of the instruments. 

Because of an electronic problem, the SEAS on the 
LPSP instrument failed after 56 days in orbit ar. all 
subsequent operations made use of the other SEAS. 
Consequently thermal and other drifts between the 
LPSP axis and the LASP axis had to be known. To 
determine these, we measure relative positions 
( 1 V. ± Y) of the solar limb (in the X • Y frame of ’he 
satellite) during one orbit using images from the 


internal raster mode. We corrected for these drifts by 
programming the secondary mirror when it was neces- 
sary to stay within 1" of the target. 

Variations in the SEAS scale factor and zero point 
(Sun-center line of sight) proved to be more trouble- 
some. A weekly determination of the absolute four 
positions of the solar limb in the X - Y frame of the 
satellite was necessary. Using such data, from space* 
craft and internal rasters, an extrapolated scale and 
zero point could be found for the particular day of 
observation. This proved to be successful, and we were 
able to position targets at the very center of our field 
of view, often without the need for corrections. Finally, 
the repeatability of the pointing system at the limb was 
found to be within I * or 2\ but with occasional jumps 
of 5". 

6) Target Acquisition 

Nearly 80^, of the orbits under LPSP control were 
dedicated to studies of selected targets, often as small 
as a few seconds of arc. This mode of observation 
from an unmanned observatory involves a fairly 
complex nroccdure which requires considerable care 
and dispatch from the observe . Wc illustrate tins in 
Figure 6 and describe the acquisition of the core of a 
sunspot. 
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As soon as a spot is visible on either a Ha picture 
(taken daily by NOAA 'NBS in Boulder) or the Ca u K. 
picture transmitted by telephone from Sacramento 
Peak Observatory, we determine its Stonyhurst 
coordinates for the time of the photograph. These 
coordinates are then transformed and oriented into 
the satellite spin frame of reference for the piojeeted 
time of observation. Finally the scale distortion of the 
pointing system is corrected for. The instrument 
commands and associated pointing commands are 
generated and sent to NASA on the day before the 
date of observation. During the day of observation an 
internal image of the spot, 64' x 64% is executed in 
the far wings of Ca n or Mg n during a real-time pass 
of the satellite over a ground station. These real-time 
data are received by telephone via GSFC on the PDP 1 1 
computer at Boulder, where they are subsequently 
decoded and displayed as an image. Any corrections 
to the position of the spot are determined from the 
image and are sent by telephone to NAS A who uplink 
the corrections to the spacecraft. The time delay 
between the real-time pass and execution of pointing 
corrections can be as short as the time interval between 
two successive ground station passes, namely, U hours. 
This apparently straightforward operation is made 
more difficult because of the need to extrapolate the 
scale distortion parameters and to correct for the 
drifts of the line of sight already mentioned. 

c) La Modulation 

The La signal was discovered to be occasionally 
modulated with an amplitude which may reach \0% 
of the signal at precisely the rotation period of the 
spacecraft wheel except during the first 2 minutes 
after sunrise w hen sometimes a period distinct ly shorter 
than the wheel period was found. 

All attempts made to detect a similar phenomenon 
in the other channels have failed, suggesting that it is 
not caused by a pointing problem. Indeed we have 
sometimes observed oscillations in the pointing axis 
with periods equal to that of the wheel rotation and an 
amplitude of 0 75, but these affect all channels at the 
same time. 

The phenomenon has to be taken into account when 
analyzing time series and profiles of the La line. As 
described in Paper I, the duration of every individual 
measurement is the product of 0.16 s gate lime by a 
power of 2, and the most commonly used values are 
10.24 s (64 grating steps) and 20.48 s (128 grating 
steps). The period of rotation of the wheel of the 
spacecraft varies from 10.7 to 9.5 s and for spectral 
scans with a time base equal to or larger than 10.24 s, 
the modulation induces a “beat 0 of period ranging 
from infinity to 131 s. 

Figure 7 plots the raw data for the first moment of 
the wavelength of the line versus time. A strong H00 s 
period is evident. If the individual data points are 
corrected for the modulation with a period equal to 
that of the wheel and the first moment is recomputed, 
then the circles in Figure 7 show that the K(K»s 
oscillation is suppressed. 



Fig. 7 — Variation with time of the position of the first 
moment of la. The solid line variation shows a strong oscilla- 
tion of ^800s period Circles represent the same data after 
correcting the La profiles for the instrumental modulation: 
the 800 s oscillations have vanished. Abscissa is time on 
seconds); ordinate is grating step number. 


This proves without ambiguity that the phenomenon 
is purely instrumental and that the possibility of its 
solar origin, should be definitely disregarded. 

V. PRELIMINARY RESULTS 

Spacecraft (playback) data are sent by NASA to 
Boulder by telephone line and recorded there on 
magnetic tapes. Mass production tapes of this data are 
prepared by NASA and mailed to the Centre National 
d*Etudes Spatiaies (CNES) (Toulouse, France), which 
is in charge of the processing and distribution of the 
final data tapes to the LPSP investigators and guest 
investigators. The results presented below have been 
obtained mostly from playback data at Boulder. 

a) Observational Program 

Table 4 presents a summary of the various types of 
observations that were programmed during the first 
8363 orbits. No distinction is made in the table with 
regard to pointing control. Orbits with problems in 
either the instrument, the spacecraft, or the command 
system represent only 2% of the total. 

The column headings in the table describe tne modes 
of operation of the instrument while the rows indicate 
the scientific question or the solar feature under study. 

Orbits labeled “ Line Profiles” usually correspond to 
spectral scans ranging from 64 to 512 grating steps 
centered at either the six chromospheric lines or the 
O vi line at 103.2 nm. The “Others” columns indicate 
studies in either the wings of the Ca ti and Mg n lines 
or the O v and Si m lines at 121. 8 nm, 120.6 nm, 
respectively. Other lines such as O i, N i, etc., were 
observed in the L0 channel by making use of the differ- 
ent grating orders. During these orbits the satellite 
was in the pointed mode. 
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Under “Spectroheliograms" we classify orbits for 
which most of the time is spent rastering, using either 
the internal or the satellite raster mode. Morphology 
and long-term variability studies of solar features were 
of interest here. 

During orbits labeled “Velocity Fields and Oscilla- 
tions*’ we spent most of the time in the pointed mode 
studying rather small areas of the disk, not more than 
64* x 1". I let c, special spectral scans vvete made. I or 
example* short -period waves (7' < 40 s) were searched 
for by scanning rapidly through the line profiles using 
wavelength position separated by 16 grating steps. 

Under “Transients" are grouped orbits observed 
with the fast, low-resolution satellite rasters. The short 
time constant of such modes allows one to observe 
rapidly propagating shocks and any other rapid 
phenomena. 

"Large A Scans" represent the full scanning capa- 
bility of the grating and were performed cither to study 
lines in the low orders of the L3 channel or to standard- 
ize the Mg ti and Ca it line cores with respect to their 
far wings. 

For most orbits, sunset and sunrise experiments were 
performed in order to measure atmospheric extinction 
at various wavelengths of interest (see $ V.e below). 

“Chromospheric Network Studies" include not 
only morphology and time evolution but also line 
profiles for center-to-limb and cell-network compari- 
sons. “Limb Studie r " include orbits dedicated to the 
shape of the limb, spicules, and the vertical extension 
of the solar atmosphere in various lines. 

“Photometric Calibration" has already been des- 
cribed in III; for these orbits the satellite is pointed 
at disk center in a quiet region. 

The last live rows of Tabic 4 inclu . scientific as 
well as instrument or satellite calibrate:, orbits. Also 
included are orbits devoted to eclipse observations, to 
studies of the geocoronal hydrogen line, and to the 
search for lines such as Fe xmi 121.64 nm and He n 
102.5 nm. 

6) Quiet Sun and Chromospheric Studies 

Figure 8 (Plate 21) represents an example of our 
study of the network and the quiet Sun. Such internal 
raster images along with associated profiles will permit 
an intercomparison of cell and network properties. 

Figure 9 compares profiles obtained at the center of 
a cel! and in a network fragment. One can notice the 
strong variation in intensity, particularly in M^ u k 
and La, between the two regions. This is due partly to 
a decreasing line background contribution as well as 
differential temperature sensitivity. The asymmetry of 
La is always much less pronounced than in the case of 
Mg n h and A, presumably because of smaller velocity 
gradients at the La height of formation. 

Tigurc 10 shows average quiet Sun profiles of La 
and Lfi taken with a 6" x 2' slit at disk center and near 
the limb (/« - 0.14). Notice the two lines of O i at 
130.48 and 130.59 nm in the wings of Ltf (observed in 
the 1 1th order). By the insertion into the light beam 
ol a Mg! j Idler one can block out all photons below 



1 15 nm and observe only the two O i lines. Differenc- 
ing permits one to reconstruct the L£ profile. The 
result of this procedure is shown on Figure 11. The 
remaining slight anomaly at grating step ISO, might be 
due to the Hi line of He n (102. 5“* nm). Center-to- 
bmb measurements at the corresponding wavelength 
show appreciable limb brightening Instigations are 
um!cr way to confirm this observation. 

The variation ot the distance between peaks of the 
L« and L£ lines from center to limb is apparent 
together with the variation of the ratio of the peak to 
core intensities (see Table 5). The values in Table 5 at 
M ^ 0.14 are close to those calculated bv Vernazza 
(1972). 

One noticeable feature is the reversed intensity 
asymmetry between the La and I i shortward and 
longward peaks at the disk center. In La the short- 
ward peak is generally higher than the longward peak 
while the reverse holds in L£ At the limb, both La 
and Li profiles become symmetrical. It is possible 
that this effect may not be intrinsic but is either an 
instrumental effect or the effect of unresolved lines at 
Li. The matter is under study. 

To study the morpho!og\ f network fragments and 
their evolution in time, monochromatic images have 
been obtained simultaneously in the six lines alter- 
nately with broad-band images in O vi (103.2 nm) for a 
number of observing sequences (Fig. 12). Preliminary 
analyMs of a 20 hr sequence shows that significant 
evolutionary changes can occur over a 12 hr period. 
The larger size A' 2v fragments are easily identified in the 
La Li ( — H Lyman) and O vi images where they do 
not appear as extended as indicated by the ATM data 
(Reeves 1976). This is supported by an analysis of the 
La brightness distribution. Applying the method of 
Skumanich. Smythe, and Frazier (1975) to both the 
present OSO 8 sequence and ATM data, one finds a 
fractional La network area of 37% and 41%. respec- 
tively. The ratio of mean network to mean cell bright- 
ness proved to be 1.9 and 2.1, respectively. For 
comparison the OSO 8 Ca n distribution yielded a 
fractional area of 27% and brightness ratio of 1.3. 


c) Quiet Chromospheric Oscillation and Transients 

We have already mentioned in § IIlc (cf. Fig. 3) our 
successful detection of the 300 s oscillation of photo- 
spheric lines in the wings of Ca n H and K. 


TABLE 5 


Center-to-Limr Comparisons of the Distance between the 
Blue and Red Peaks and of the Ratio between the Average 
Peak Intensity and the Core Intensity of the La and 1 $ 
Lines 

SFPAR ATION OF 

Pi aks (nm) / p „ h A:or. 


I.INI |i = I O ft O 14 01 = I 0 0 . =■ 0.14 

La 0 04 3 0 05 1.3 + 50% 2 2 

l (J (H>27 noil 1.4 + 50% 1.5 
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Fig. 9 - Companion between simultaneous observations of Ca n K, Mg it k s and La profiles in the network and at the center 
of a supcrgranulatton tell. The slit site is I" x 10'. Units arc counts per counting gate. 




The Mg it k and Ca ii K lines (Fig. 13) show oscilla- 
lions of **200 s period. We parametrize the profile of 
Mg n k as a difference of two Gaussians: 

/(A) - E, exp (~^)*[ I - A, exp (—-*)*] • 

where /T, and A, rcprcsenl the intensities of the emis- 
sion and absorption components of the profile, respec- 
tively, and /* a and A 2 the average wavelength position 
(or ** first moment of the wavelength ") of the emission 
and absorption components. We should point out 


that no physical meaning is to be attached to this 
parametri/ation. The time behavior of these various 
parameters is shown on Figure 14. 

The oscillation of the parameter £ 2 or average emis- 
sion position covers a range of ±0.001% nm or 
±2 kins 1 . For A 2 or average ubsorption position, 
this value is doubled. This clear difference might be 
interpreted as the amplification, with increasing height 
in the atmosphere, of a wave, if one assumes, a*, usual, 
that the core of Mg n A is formed at a higher altitude 
than the wings. 1'he average ** bluer" position of £* 
(emission) compared to A, (absorption) reflects the 
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Fig. 14. — Time variations of the various parameters of the 
analytical function used to represent MgttA. Observations 
are made at disk center. 


fact that the average Mg n k profile in the Sun is 
asymmetric, with k 2 v stronger than knr. This is 
illustrated also on Figure ! 5, where we see that the line 
shape varies from a strongly asymmetrical profile 
(* 2 i > h 2 r) to a nearly perfect symmetrical profile 
(k 2 v - frjT). 


In Figure 16 we see the variation over 40 $ of the 
Mgii and Ca ti K profiles observed simultaneously 

The Mg ii lines exhibit the same periods as the Ca it 
lines. We have a broad set of obsemtions which con* 
tain wave trains tasting in general no longer than a few 
cycles with periods ranging from 250 s down to 1 30 s. 
A search for shorter periods was undertaken, in 
particular by guest observers, but no clear evidence 
has yet emerged from these investigations at this early 
stage of data analysis. 

The good results obtained in Mg n encouraged us to 
search for a possible oscillation of La. The large 
contribution function of the line, which tends to 
smooth out the effects of any wave on the profile, 
together with the low photon count in this channel 
made this observation a particularly difficult one. We 
first tried to detect intensity fluctuations by integrating 
the number of photons over ±0.025 nm from line 
certer. We did not find any obvious evidence of varia* 
tions, other than random, a result in accord with the 
previous attempts made from studies of Skylab results 
(Vernazza et at. 1975). 

To overcome the low photon statistics problem, we 
tried to correlate the shape of La with that of Mg li k 
profiles. We definitely see evidence for a correlation, 
the bluer k profiles corresponding to redshifted La 
profiles. The amplitude of the shift is of the order of 2 
grating steps (3km$~ l ). More work is under way, 
but we may state at this stage that the oscillations seen 
in Mg ii k have an influence higher in the chromo- 
sphere, at the altitudes where La is formed (Artzner 
etai 1978). 

Several orbits were devoted to the study of oscilla- 
tions in the O I, Si lit. and O vi lines, but have not yet 
been analyzed. Transient and short time phenomena 
have been observed. Numerous tachograms dedicated 
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Fio. 1 5.— Mg it k profiles observed at the maxima and minima of the oscillations of the line center or gravity, showing obvious 
occurrence of a>ymmctrical blue peaked profiles at maxima and symmetric-^ or slightly asymmetrical red peakett profiles Jt minima 
Wavelength* nutcase to the left 
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Fio. 16. — Observation of a temporal event, simultaneously 
in Ca n K and Mg w k as it propagates through the line profile. 
Wavelengths increase to the left. Vertical scale is in relative 
units. 


to the investigation of transients and other flows in the 
chromosphere have been programmed and are under 
study. 

d) Study of Sunspots and Active Regions 

Sequences of high -resolution images and spectra 
taken r apart show the evolution of spot morphology 
and profiles in space and time. Umbral flashes have 
been observed as well as oscillations. Nearly simul- 
taneous observations of lines from the photosphere to 
the transition region should allow us to study the 
vertical structure of spots. 

A preliminary analysis of spacecraft rasters made 
with the r x iO' entrance slit, and LPSP internal 
rasters made with the T x T, T x 3' entrance slit, 
show f * it regions of strongly enhanced O vi 10X2 nm 
emiss, , very often tend to be distributed other than 
directly above sunspot umbrae. This is the case both 
for some unipolar single spots as well as multipolar 
spot groups. However, we do find circumstances when 
single sputs show the transition region "plume" 


directly over the spot as reported by Foukal et at. 
(1974). 

Figure 17 shows this phenomenon, where isophotes 
are represented in photospheric light. La, and O vi. 

We have also found examples of steplike changes 
in the distribution of the enhanced O vt emission 
(Fig. 1S>. If we avsume that the enhanced 0 vi emission 
is associated with a specific magnetic field connectivity, 
then our results imply “step** changes in field connec- 
tivity. The general O vi emission in the active region 
was found to be approximately 10 times brighter, with 
strongly enhanced features approximately 100 times 
brighter, than the average quiet Sun. Simultaneous 
observations in La, L0, Mg u, and Ca n also show 
the general active-region enhanced emission as we)! as 
strongly enhanced features. The La and L£ features 
are identical, but show systematic horizontal displace- 
ment with respect to the Ca n and Mg u emission 
features. Little or no correlation is found with the O vi 
structures. With regard to the profiles of the resonance 
lines, the self-reversal is w*ak or absent above regions 
as shown on Figure 19. Presumably the line-forming 
region has reduced opacity. 

e) Studies of Prominences 

Temporal evolution of, and velocity field in, promi- 
nences were studied with consecutive monochromatic 
spacecraft rasters. Figure 20 shows, at the top, Ca it 
images at wavelengths ranging from -0.012 nm to 
+ 0X)3 nm from line center. The main loop is clearly 
visible only at line center; this is confirmed by spectra 
constructed from the different monochromatic images 
of the prominence which show a single emission peak 
with a FWHM of 0.020 nm which is typical of quies- 
cent prominences (Engvold and Livingston 1971). This 
is consistent with a mean turbulent velocity of 
9 km s' 1 and a temperature 7* e * 8500 K. At the 
bottom of Figure ‘20 are two simultaneous La and L0 
images. 

Several observations were made to study the thermal 
structure and evolution of active and eruptive promi- 
nences. Figure 21 shows an internal raster performed 
with a slit of r x P above active region McMath 
No. 14127. 

The first three simultaneous images in La, L0, and 
Ca it K a v show a loop system with a rather faint 
contrast in K 2 v as compared to the plage at the limb, 
but a higher contrast in La (and L0). The next three 
images are separated by 1 1 minutes (raster repetition 
rate) and arc made in O vi 103.2 nm. They trace out 
the high-temperature evolution of the region (O vi 
being formed at approximately 350,000 K ). The first 
O vi image shows a faint high loop, the second some 
residual brightness around the foot of the loop, and 
the third an important enhancement at this same point. 
Then within 20 min the loop disappeared, as can be 
seen on the next three images in La, L/3, and Ca n K. 

Many images and spectra have been obiainc * that 
will allow the study of temperature, density, and 
velocity variations during loop evolution (Vial et ai 
1978).' 
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Fig. !$.— illustration of step changes in held connectivity of spot evolution. The two upper pictures are sr 'tliste rasters made m 
the wings of Mg iU. The two tower pictures are taken simultaneous!} in O vt. The evolution of the sumpc*. group is continuous 
when observed m white light and discontinuous m O vt 
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Figure 22 represents observations of an active 
region at the limb on 1975 July 7. The slit was parallel 
to the limb and probably intersected it slightly, as 
indicated by the presence of scattered light in the 
wings of the Ca it profiles. On the left portion of 
Figure 22, full and dotted lines correspond to profiles 
taken at positions separated by onK l\ The emission 
maxima in Ca u and Mgn are displaced short ward with 
an amplitude corresponding io about 20 and 1 5 km s’ 1 , 
respectively. These maxima may in fact correspond to 
K 2 v and kyv f while IC 2 r and k 2 r are barely visible and 
show up only as asymmetries in the line profile. La is 
not displaced and shows the geocoronal reversal with 
no aprarent self-reversal. The distance between its 
emission peaks is only 0.022 nm, leading to an optical 
depth of only 500 at line center. On the right panel of 
Figure 22, solid and dotted lines show temporal 
variations for a 22 min interval. Orbital Doppler 
effects have not been corrected for on the figure; but 
if they are corrected for. a shift to the longward is 
still found for Ca n and Mg n with an amplitude of 
5 km s’ 1 . Moreover, Ca ii and Mg ii lines have 
essentially the same intensity while La has increased 
by a factor of 2. 

/) C romosphcre-Corona Transition Lines 

The two lines Si ill 120.65 nm (3 s l S-3p l P°) and 
O vi 103.19 nm (2 s 2 P-2p 2 P°) are formed in the 
chromosphere-corona transition region at 40.000 K 
and 350,000 K, respectively (Jordan 1969). They are 
observed with the LPSP instrument with a spectral 
resolution of 0.002 and 0.006 nm, respectively. 

The shape of transition-region line profiles may 
indicate whether there is any propagation of either 
acoustic or magnetohvdrodynamic waves (McWhirter 
1977). The presence of such waves may he sympto- 
matic of coronal heating mechanisms. 

Figure 23 shows an average (single orbit) Si 
profile at the center of the disk (quiet Sun) observed 
with a resolution of r x 40*. The FWHM is 0.015 nm 
and, if we assume that the line is optically thin, the 
rms (line-of-sight) nonthermal velocity is 22 kms' 1 . 
This value is 4 km s’ 1 higher than the values obtained 
by Nicolas et ai (1976) for the Si lit lines at 128.9 and 
189.2 nm. 

Quiet and active Sun profiles of the O vi line at the 
center of the disk are given on Figures 24*7 and 24 b. 
This line is optically thin, and the KWH M is 0.021 nm 
nearly identical for both quiet and active Sun; the 
rms line-of-sight nonthermal velocity is 30 km s’ 1 . A 
departure from a purely Gaussian profile can be 
noticed in both profiles. The line appears asymmetric 
and may indicate the cITcct of a velocity structure in 
the region of formation of the line. 

Figure 24r shows a quiet limb O vi profile, averaged 
over several positions above the limb ( + 2* + 6"). The 
FWHM is now 0.026 nm, equivalent to a rms line- 
of-sight nonthermal velocity of 3Kkms‘ 1 . This is 
larger by II kms' 1 than the value previously quoted 
by Moe and Nicolas (1977). 

The chromosphere and transition region heittht 
distnhution is shown on f igure 25 as derived fioin 


TABLE 6 

Characteristics of Species in the Earth's atmosphere 
Observed with the LPSP Instrument 


Range of 
Altitudes m 
the Earth 

OSO $ Atmosphere 

Channel Component Wavelength (nm) ikno 


Mg H 03 277.7-282.3 60-80 

Mg It OH 281.6 >75 

La 02 120.6-122.3 90-110 

La H 121.6 >500 

Lfi 02 101.7-103.2 150-200 

L0 H 102 5 >500 


internal raster scans simultaneously in the wing of 
Mg ii A, Mg ii /r 3 . O v 121.8 nm, and N l 1 19.95 nm. 

As is apparent from the figure, the A 3 chromosphere 
appears to have an additional contribution, which may 
likely be due to spicules that appear to peak at about 
the same height as the O vi component. Separate 
measurements of O vi (and for the far wing of Mg n h) 
show a similar behavior as O v. These results would 
argue for an inhomogeneous transition between 
chromosphere and corona (cf. Doschek, Feldman, and 
Tousey 1975). 

g) Aeronomy Investigations 

At orbital sunsets and sunrises, the solar UV light 
is absorbed by successively denser layers of the 
Earth's atmosphere. Vertical distribution of number 
densities of several components ma\ be studied by this 
technique as indicated in Table 6. The light in the Ca 
channels is no: attenuated and provides a pointing 
reference. 

The measurement of the width and depth of the 
hydrogen geocoronal absorption was undertaken 
during several orbit days. This is a new and very 
promising observational technique to measure simul- 
taneously the exospheric temperature at each point 
of the orbit and the atomic hydrogen density at the 
exobase, which may solve the question of what 
mechanism(s) control(s) the hydrogen distribution at 
the exobase. The preliminary results of aeronomy 
investigation have been published in Vidal-Madjar 
et ai (1976). 

VI. CONCLUSION 

We have described here the actual performance of 
one of the most complex solar physics instruments 
launched into space and the main results obtained with 
the first high-resolution multichannel UV and visible 
spectrometer placed in orbit by OSOti. For the first 
time, an absolute pointing accuracv of nearly T 
could be achieved in orbit with real time opera! ions. It 
undoubtedly represents the largest and most complex 
experiment of the French space program in solar 
physics. Although the bulk of the data has not yet 
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Fig. 19 — Profiles of Ca n K. Mg n k. Le, and LS observed simultaneously over an active region All lines are free from the centra! 
self-reversal usually observed over quiet Sun areas The two peaks in the L/i profile are the two O i lines. For conversion into A 
units, see Table I. Wavelengths increase to the left 


been examined in detail, preliminary analyses show that 
the performance of the instrument was nominal and 
at times be>ond nominal expectations. The results 
presented here are only isolated samples of what has 
been obtained in the first 18 months. The instrument 
continues to perform nominally and has begun its 
third year of operation. This will allow us to obtain 
more data on act ; ve regions and particularly on flares 
which were very rare during the first 18 months. 
Indeed, only one flare, that of 1977 April 19, has been 
observed so far (Jouchoux ct al. 1977). 

The operation of the instrument has been very 
exhausting, and we have benefited from the assistance 
of many people. Our experience with regard to the 
remote management of an entirely automated complex 
instrument is, we feel, of great value for similar 
experiments in the future. 


The accomplishment of this experiment would not 
have been possible without the support of CNES and 
particularly of Dr. A. Lebeau, former Director of 
Programs and Plans, and Professor M. Levy, former 
President. Wc would like to thank collectively the 
NASA and CNES engineers who have contributed to 
this experiment. 

The operations of the instrument from Boulder 
would not have been possible without the kind 


hospitality of LASP, in particular of its Director, 
Professor C. Barth. We also thank the LASP OSO 8 
staff for their contributions. We are indebted to NASA 
and LASP for access to space on the American 
calibration rockets. The excellent spirit of cooperation 
and the dedicated service of the OSO 8 Control Center 
at the Goddard Space Flight Center was certainly a 
key to the success of real time acquisition and of the 
daily programming of observations in general. The 
observations of sunspots, active regions, and flaring 
regions could not have been done without the generous 
assistance of NOAA, Big Bear Solar Observatory, 
Meudon Observatory, Sacramento Peak Observatory, 
and Lockheed Research Laboratory. We wish to 
express our warmest acknowledgements to these 
numerous and often anonymous people who played 
such an important although thankless role in the 
daily work required by the continuous observation of 
the Sun during several years. Highly appreciated were 
the contribution of Drs. P. Bruston and M. Malinovsky 
of LPSP in the preparation and checkout of observing 
programs. Invaluable and continuous support in the 
daily operations was provided by M. Brusion (Mrs.) 
and N. Dionnau of LPSP. We also thank J. Borsen- 
berger ( Inst it ut d'Astrophysiquc dc Paris) and B. 
Phissamay ( LPSP) for their important contribution. 
Last but not the least, all ofthc Guest Investigators 
who have assisted our team in the operation and 
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Fig. 2?.— Profile of the Si lit 120 6 nm line observed at Sun center with a resolution of I* 
left. Hor conversion into A units, see Table I. Counts represent an average over one orbit 
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QUIET ACTIVE LIMB 

Fig. 24.— io) Quiet and (6) active Sun profiles of O vi at Sun center. Each profile is an average of 20 individual profiles. The full 
lines represent the observations; the dashed line, a best Gaussian fit. The resolution is 1* x 40', Wavelengths increase to the left. 
For conversion into A units, see Table I . Profile tr) refers to the limb; it is also an average, but the individual profiles were observed 
at another orbit, and the units cannot be compared with those of profiles (a) and (6). 
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Fic. 25.— Limb "darkening" profiles, in Mgu A3, the far wing of Mg it k, Ni 119.95 ntu, and O v 121.83 nm. Distance is 
measured in arcsec (arbitrary origin). 
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ABSTRACT 

/ 

Observations indicate that .he temporal variations of wavelength of the reversal of the solar 
H La and Mg k lines are correlated. 

Subject headings: Sun: atmospheric motions — Sun: chromosphere — Sun: spectra 


I. INTRODUCTION 

The solar photospheric and chromospheric lines 
exhibit ubiquitous quasi-period ic oscillations of inten- 
sity and wavelength. The presence or absence of 
periodic fluctuations higher in the solar atmosphere has 
been investigated by measuring the intensity variations 
of spectrai lines formed between 10 4 K and 1.5 X 10* K 
(Wmazza ft at. 1973). the micr wave emission of the 
Sun around 10* K (Avery 1976), the intensity and 
position of the C ii 1356 A line formed around 20.000 K 
(Chipman 1977), the intensity and position of the C iv 
1554) A line formed around iOO.OOU K (Bruner 1977), 
and the line profile of the Fe xiv X5 *03 corona! une 
(Tsuhaki 1977). 

We report here the first observations of the solar 
HLa 1216 A lire formed near 25.000 K with temporal, 
spectral, and spatial resolution adequate for the study 
of solar velocity fields. 

u. observations 

The LPSF instrument on hoard OSO 6 is described 
by Artzner et at. (1977). The observations reported here, 
mostly at disk center, had * l’X 3", \" X 10", or 
1" X 20" spatial resolution, a spectrai scanning incre- 
ment of, respectively, 2.4, 2.3 and 1.6 km s* 1 *or the 
La, Mg k -* 1 Ca K lines, and a lime resolution of from 
10 to 40 s. A sequence of observations consists of 100 to 
180 successive spectra scanning the I-a line over ±0.53 
A, the Mg k line over ±1.4 A, and the Ca K line over 

± 1.1 A. 

In order to discriminate between actual solar oscilla- 
tions and pointing variations, we have measured the 
pointing stability by three different methods (Bonnet 
et a!. 1978). As a result, we estimate that the pointing 
drifts randomly at a rate smaller than 1" per 3 minutes 
of time. 

As a chetk of the instrumental stability over 1 hour 
and of the procedure to com|>cnsatc for the Dopplcr- 
Fizca j component of the velocity of the space r aft, wc 
have ..veragvd Lo s|»ectra duri. g three consecutive 20 
minute intervals of the 60 minute daylight |>ortinn of an 
orbit, llu constant jmsition of tile geuroronul line 
demonstrates th t the wavelength drifts by less than 

* Lockheed Research Lalioraiurics, I'alu Alto. 


4 mA over one orbit. The observed spectral resolution 
does not vary, but the shape of the geocoronal absorp- 
tion is wider during the third part of the orbit day than 
during the first part. This measured geophysical effect, 
most likely due to the heating of the terrestrial atmo- 
sphere between local sunrise and sunset, will be in- 
vestigated elsewhere and gives us confidence in the 
ability of the instrument to detect minor changes in the 
shape and position of spectral lines during 1 hour. 

The low measured dark current (1.1 counts s~*) 
enables us to make use of data taken through orbit 
380, when the sensitivity was 35 counts s~* at the blue 
La peak with 0.01 A and 1" X 10" resolution. 

Note that the simultaneously measured counting 
rates may vary from 10 (La wing^ to 100,000 (K2v 
peak of Ca K), because of the combined increase of 
solar dux and instrumental sensitivitv from 1216 to 
4000 A. 

in. data analysis 

Anv plot of th. data exhibits quasi-periodic oscilla- 
tions of the Ca K, Ca H, Mg k 2795, and Mg h 2803 A 
lines. These are reported in Artzner et ai (1978). As 
indicated by Bonnet et al. (1978), the La signal for a 
constant solar input is modulated in phase with the 
rotation of the spacecraft wheel, with a period of the 
order of tO s. Therefore a mathematical procedure (as 
used by Artzner et a /. 1978) applied to the raw, noisy, 
uncorrccied La data will fail to re- .‘a! the true time 
behavior of this lin*. because it will essentially reflect 
the “10 second modulation” of the signal. As the mean 
intensity of the La line, measured over 10.24 or 20.48 s, 
is much less sensitive to this modulation than the 
velocity measurements, we have integrated the La 
profiles over ±0.41 A and we have computed the aver- 
age power spectrum (for 16 sequences) of the time 
variation of this quantity (Fig. 1). The result agrees 
with the observations of Vcrnazza et at. (1975). For 
each sequence wc have measured the standard deviation 
of • he values of the integrated intensity (Table 1). For 
the sequences with low average numbers cf counts per 
profile, the intensity variation is found to be primarily 
due to statistical noise, but for the sequences with 
higher count rates, the standard deviation is consider- 
ably larger than photon statistics alone would imply. 
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Pic. 1 —Temporal variations of the La intensity; bandwidth 
±0 41 A. Average of 16 power spectra. 

TABLE l % 


La 1216 A Line (bandpass ±0.41 A) 


Average 
Numlier 
Duration of Counts 

(minutes' per ProtVe 

Standard 

Deviation 

around 

Average 

40 

200 

23 

49 

220 

14 

55 

2 .W 

16 

42 

240 

13 

35 

270 

16 

2$ 

290 

U 

3“ . 

2>m 

20 

37 

300 

21 

39 

350 

16 

53 

340 

20 

40 

350 

^7 

42. ... 

370 

h 

56 . . 

400 

24 

47. . . 

450 

62 

53 

520 

35 

56 

520 

2V 


As of now, we have analyzed 21 orbits of the La data. 
The analysis has been conducted in such a way as to be 
insensitive to the 10 s modulation : we classify into three 
classes the Mg k spectra according to the relative 
strength of the k 2 peaks and then compare average La 
spectra corresponding to each of the three classes. Class 
I has a high ratio of (Mg k 2r intensity/Mg intensity), 
while classes II and III have successively lower values 
of this ratio. This procedure reduces statistical noise 
and should reduce the effect of the “10 second" modula- 
tion, because the solar oscillations observed in Mg k 
and the spinning of the spacecraft should have a 
random relationship. As a result, after correcting for 
the Doppler component of the velocity of the space- 
craft, average class I, II, and III profiles are computed, 
not only for the Mg k line but alfo for the simultaneously 
observed H La, Ca K, Ca H, and Mg h lines. 

1. For the Mg k channel. % *crage class I profile 
actually has a strong (blue k • ■ , red peak) ratio, as a 
verification of the procedure. 

2. A very similar effect is seen for the Mg k t Ca H, 
am! Ca K channels. 

3. For the I .a line, the average class J and class 111 
profiles are identical in the wing, but shifted in the 
central reversal. The integrated intensity of the average 


“class 1“ La profile is equal to the integrated intensity 
of the average “class III” La profile, whereas the 
integrated intensity of the average “class I" Me k 
profile is stronger than for the average “class III" 
profile. 

We have fitted to the La and Me k line profiles a 
six-parameter formula. /, Xc, HY, .1, Ag, H u, 

1M - , eI p[-(^)’] ^ 

x|>-'«p[-C-ir)’]S 

The parameters Ac and IFe refer to the wings and the 
parameters Xu and IFe to the position and width of the 
central reversal. 

The data ,oints within ±35 mA of the geocoronal 
line were excluded from this fit. The 2 ero of the wave- 
length scale is, for the La observations, fixed by the 
geocoronal absorption (no correction for the solar rota- 
tion is needed, as the observations reported in Tables 2 
and 3 are at disk center). The zero of the wavelength 
scale for the Mg k line, at the present stage of data 
reduction, may be affected by a systematic error of 
±0.01 A. For the sake of comparison, the wavelengths 
have been conx-erted in meters per second. The numl ers 
quoted do not imply that such velocities are presen; in 
the solar atmosphere; the derivation of actual solar 
velocities is be>ond the scope of this Letter. 

Table 2 shows the results of these fits for the three 

TABLE 2 

Position and Intensity of Spectral Features 



of Average Spectra* 



High 

(Mg*.. Mg*.,' 
Ratio 
Class I 

Average 

Ratio 

CW II 

Low 

(Mg*j. Mg*lr> 

Ratio 
Chf* III 

La 

Xo-H telluric 
(central 
reversal) . . . 
Xf-He telluric 

(wings) 

It 

. +3800 ms“‘ 

. +1100 ms“* 
8 9 

+2300 m %* 

+ 1100 ms-* 
8.9 

+950 m r* 

+780 rr s* 1 
90 

Mg k 

Xo central 

Xe emission 
Arbitrary units 
it emission . . 

. +4150 ms* 1 
+330 m s’* 

315 

+2950 m r » 
+ 730 ms*' 

283 

+1720 m*" 1 
-t- 1350 in s' 1 

269 


Famt Mn i 
photospheric 
line ter<- 

average +75 m s* 1 -25 m s' 1 +50 m s~ l 

M geocoronal 
ab*orpt‘ on 

(aero average) —60 mC 1 M m s“* 


• C*iwifictl according t«» the K|,/Kj, ratio; three equal claw* 
for each time sequence. 
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classes of profiles defined above. From the same set of 
data, by averaging of faint and bright sequences, we 
have computed typical La and Mg k faint and bright 
profiles, and applied the same fit with a six-parameter 
formula (Table 3). From the comfwrison of Tables 2 
and 3, it ajipuars that, for the I-a ami Mg k lines, if the 
iwsititm of the centra! reversal is related to i lie atmo- 
sjdieric velocity fields, the instantaneous, time-resolved 
velocity fluctuations are greater than the large-scale 
(chromospheric network) permanent velocity features. 

We have also computed three average profiles by 

TABLE 3 


Position and Intensity or Spectral Features 
of Averace Spectra* 



Average of 5 
Brightest Sequences 
out of 16 

Average of 5 
Faintest Sequences 
out of 16 

La 

Xfl-H telluric central 

reversal 

Ar-H telluric wings . . 

!t 

4-2600 m s-> 
4-900 m S-« 
12.5 

4-2100 m s" 1 
-r330 m s* 1 
5,4 


Xa central reversal . . 
Xr emissior. . 
it 

-*►3000 m s“ l 
-r 1300 m s" 1 
391 

4*3060 m s* 1 
~450 m sec 
241 


* Brightest sequence versus faintest sequence. 


sorting out at random the observed profiles into three 
classes to derive an estimate of the statistical precision 
of our measurements. We find the La, photospheric 
line, and geocoronal line velocities are accurate to about 
± 100 m s~\ while the Mg k velocities are accurate to 
±25 m s The La intensity values an accurate to 
±0.1 in our units, and the Mg k intensities to r 2. 

IV. CONCLUSION 

At this stage of OSO 8 data reduction, we cannot yet 
present the power spectrum of the spatially resolved 
solar La line velocity fluctuations; nevertheless, we 
have demonstrate that the La centra! reversal does 
exhibit wavelength fluctuations positivelv correlate 
with the oscillations of the chromospheric Ca K and 
Mg k lines. 

Temporal variations in the derive velocity of the 
central reversal feature of 2500 m s' 1 in Mg k are 
accompanie by variations in the same sense of approxi- 
mately 2800 m s'* in the central reversal feature of La. 

This Letter is the continuation of efforts of teams on 
the OSO 8 project a; Verrieres-le-Buisson, Gre^nbtlt, 
Los Angeles, and Boulder. Special acknov ledgme: k goes 
to A. Jouchous, who o! .:ned special observing se- 
quences for this program. The computations were 
carried out on the CDC 7600s of XCAR and CNES. 
The CNES funded the fabrication and operations of the 
instrument under contracts 70-220. 71-202. 72-202, 
73-202, 74-202. 75-202, 76-202, and 77-202. 
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A Dynamical Representation of the Solar Qiromosphere , J.W. 
IZIBACKEP, Lockheed and P. GOUTHBROZE, Laboratoire De Physique Stellaire 
et Planetaire * Me present investigations of the non-linear hydrodynamic 
state of a model solar atmosphere and the time dependent, ionized calcium 
and magnesium resonance line profiles emitted by the chromosphere. We 
calculated the one dimensional response to an excitation 1.6 Mm below the 
visible surface using a finite difference code, and the resulting motions 
exhibit the well known 300 second oscillation of the photosphere and 200 
second oscillation of the chromosphere. Both oscillations have the chracter 
of standing waves. We shall discuss the formation of these oscillations and 
their dependence upon parameters of the calculation. Time sequences of 
resonance line profiles display asymmetry and intensity variations that result 
from both oscillations. This work was supported by NASA contracts NAS8-32356 
and NASW-3053, the Lockheed Independent Research Program and the Centre 
National d’ Etudes Spatialies. 
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ABSTRACT 

i wo of the most representative chromospheric lines. 

Mg II k and Ca II K are used to study the formation of optical- 
ly thick lines In a time-dependent, one-dimensional model of 
the solar atmosphere, which was discussed by Leibacher, Gouttebroze 
and Stein . (1979). Time sequences of these line profiles 
are calculated for two kinds of atmospheric motions : propagation 
of a pulse through the atmosphere and free oscillations. The 
mechanisms of formation (especially the displacements of the 
emitting layers) are studied for different parts of the profiles. 
Finally, the deformations of the profiles are analyzed using 
methods also suitable for observations, and the resulting 
parameters are compared to physical variables in order to 
evaluate the diagnostic methods. 



INTRODUCTION 


Among the known motions of the solar atmosphere, the so- 
called 300 second oscillations, which are now known to originate 
from the sub-photos pheric regions of the sun, have been the 
object of numerous studies. On the contrary, while some early 
observations of the Ca II K resonance line (Jensen and Orrall, 
1963) showed periodic variations in the core, formed in the 
chromosphere, with shorter periods (170 seconds), we remain far 
from an understanding of these chromospheric oscillations. 
Nevertheless, the observations of the Ca II K line continued 
to provide knowledge about chromospheric motions (see for 
instance Wilson and Evans, 1971, Liu, 1974, Beckers and 
Art zner, 1974, Cram, 1974). More recently, chromospheric motions 
have been seen in other lines : especially U.V. emission lines 
with the 0S0-8-Colarado experiment (see, e.g. Chipman, 1977, 

1978, White and Athay, 1979a, b, Athay and White, 1979) Mg II 
and L a lines with the 0S0-8-LPSP experiment (Bonnet et al . 

1978, Artzner et al^. 1978). Among the theories developed to 
interpret these observations, the more detailed are probably 
those of Heasley (1975) and Cram (1976) who studied the deform- 
ations of the Ca II lines during the propagation of a pulse 
through the atmospeh. . Recently, Durrant, Grossmann-Doerth 
and Kneer (1976), discussed the temporal evolution of Ca II K 
line profiles and concluded In favour of a mechanism of 
chromospheric oscillations similar to that proposed by Leibacher 
(1971). This encouraged us to undertake a more detailed study 
of chromospheric hydrodynamics (Leibacher, Gouttebroze and stein. 



1979* hereafter Paper I) and to compute the emergent intensities 
for various spectral features relevant to chromospheric studies. 
These features are conti nua (Lyman continuum or centimeter 
microwaves), optically thin emission lines, and optically 
thick lines. The other aim of this study was to provide the 
framework for an interpretation of the optically thick lines 
(Mg II, Ca II and L a) observed with the 0S0-8 LPSP experiment, 
which is planned in a further paper of this series. In the 
present paper, we Investigate the Mg II k and Ca II K resonance 
lines and discuss variations of these lines produced by either a 
propagating pulse or free oscillations. 

II - THE INITIAL SOLAR-TYPE ATMOSPHERE 
As a first step in the computations, we must define a 
solar- type atmospheric model in hydrostatic equilibrium, 
providing a medium for the propagation of the waves. This 
atmosphere must cover a wide range of depths in order to include 
both the region where the five-minute oscillations originate in 
the convective zone and photosphere, and the whole chromosphere, 
which is the object of the present investigation, from the 
temperature miriimjm up to the base of the corona. The definition 
of such a model atmosphere Is subject to two general constaints : 
those arising from the frequencies and amplitudes of the actually 
observed oscillations, which will be referred to as "hydro- 
dynamical constaints", and those provided by the intensities 
observed in the various lines and continue of the solar spectrum, 
which may be summarized as "radiative transfer constraints". 



a - HydAodyviamical conA&uUnti 
First of all, the variation of the temperature with height 
in the convective zone is prescribed by the outward flux of 
energy, following the mixing-length theory. Once this variation 
is known the position of the bottom of the atmosphere (assumed 
rigid for computational purposes), is determined by the *i ve- 
in inute oscillations. Five minutes represent four tires . the 
time of propagation of the waves between the bottom and the top 
of the convective zone (the top corresponding to the photos pheric 
layer whose acoustic cut-off period is shorter than five minutes). 

A similar requirement exists for the chromosphere. As 
explained in Paper I, the eigen-frequencies of chromospheric 
oscillations are sensitive to the temperature structure of the 
chromosphere. The observed periods of oscillation are not so well 
defined for the chromosphere as for the photosphere or the 
convection zone, and range from about 150 to 250 seconds : For 
Instance, Jensen and Orral 1(1963) found a period of about 
170 seconds for the K 3 absorption feature of Ca II, Liu (1974) 
a rather broad distribution about 180 seconds, and Beckers and 
Artzner (1974) a repetition rate of 200 - 50 seconds. These 
variations may be attributed to the inhomogenity of the solar 
chromosphere. 

b - Radiative. Pian6^eA c omtMuwU 
The solar spectrum Is known with some accuracy over a 
large range of frequencies, which allows the determination of the 
average temperature structure of the photosphere and chromosphere, 



through the use of semi 'empirical models. These methods and 
some recent models have been reviewed by Avrett (1977). Some 
discrepancies remain between different semi •empirical model 
atmospheres based on different kinds of data, which allows 
some degrees of freedom In the choice of atmospheric para- 
meters. For instance, models based on continuum data, such as 
those of Gingerich et al_. (1971) or Vernazza, Avrett and 
Loeser (1976) exhibit temperature minima between 4100 and 
4200 K. On the other hand, atmospheric models derived from 
the shape of the Mg II or Ca II resonance lines, taking into 
account partical frequency redistribution effects, have a some- 
what higher temperature minimum of about 4400 or 4500 K, as in 
that of Ayres and Linsky (1976), 

In the higher chromosphere, the increased importance of 
inhomogeneities gives riseto larger variations between the 
various temperature models. Above 9000 K (the upper chromosphere) 
a constraint is provided by the non-LTE calculations of hydrogen 
ionization required for the determination of electron and hydrogen 
densities in hydrostatic equilibrium. These computations include 
the determination of intensities in the Lyman a and B lines, 
which may be compared to observations (owing to inhomogeneities 
in these regions, an agreement within a factor of 2 or 3 for 
the Integrated intensities may be considered as satisfactory). 

In order to reproduce the observed width of chromospheric 
lines (and especially the separation between the red and blue 
peaks of the Mg II or Ca II lines), our atmospheric model must 
include a non-thermal isotropic velocity field which increases 



with altitude. These motions produce a "mlcroturbulent" pressure 
whleh must be Included In the hydro-static or -dynamic equations. 

In order to get a self-consistent representation of the 
atmosphere. It should be noticed that the Inclusion of this 
effect may alter significantly the depth scales In the 
chromosphere, and consequently the eigen-frequencies of oscillation. 

c - ComputcutionoUt mytkotiU 

Our atmosphere Is defined by the electron temperature (Te) 
and mlcroturbulent velocities (v T ) as functions of mass Integrated 
along a vertical line (m). From this basic set of parameters, 
we obtain the corresponding altitudes, densities, pressures* and 
el ectron-to-hydrogen density ratios by solving iteratively the 
equations of hydrostatic equilibrium, hydrogen and metal ioniz- 
ation. A difference-equation method is used to solve the equations 
of hydrostatic equilibrium, taking into account the gas, 
electronic and mlcroturbulent pressures. To determine the 
Ionization of hydrogen, we consider a 3-level and continuum 
atom, and solve self -consistently the equations of statistical 
equilibrium and radiative transfer in the L o, L e, H 0 and 
Lyman continuum transitions. The photoionization rates in the 
Balmer and Paschen contlnua are fixed by setting the rad'* Mon 
temperatures to 5000 and 4800 K, respectively. In the parts 
of the atmosphere which are optically thin in these transitions. 
Deeper In the atmosphere, they are assumed to be equal to the 
local electronic temperature. The effects of partial frequency 
redistribution are taken into account in the computation of 



the Lyman lines. The otter elements contributing to the 
electron density are treated in the LTE approximation 
(their contribution Is really important only In the temperature 
minimum region, where hydrogen is scarcely ionized). 

d - The adopted model 

As explained above* the present hydrodynamics! calculations 
require a model atmosphere extending from deep within the 
convection zone to the low c or or a , which constitutes a range of 
depths far more extensive than is usually the case for radiative 
transfer computations. For this reason, the atmosphere used here 
is a synthesis of several models frem different origins, with 
local modifications. First, our convective zone is taken from 
Baker and Temesvary (1966). Our photosphere and chromosphere 
are similar to those of Vemazza, Avrett and Loeser (1973 and 
1976) with the following modifications : first, the temperature 
minimum region is raised to 4500 K, following Ayres and Linsky 
(1976), in order to get a better fit to the observed Ca II and 
Mg II resonance line w*ngs. The temperatures are slightly lowered in 
in the upper chromosphere (by less than 200 K) and the so-called 
Lyman e-plateau (by about 2000 K) for two reasons : first, the 
large temperature excursion/ , roduced by waves in these regions 
give strongly non-linear variations of the source functions, so 
that the averaged Intensity is higher than the intensity 
corresponding to the average-temperature model ; secondly, 
some dissipation of mechanical energy occurs in th?se regions 
(provided that the waves have sufficient amplitudes) which 



produces an Increase of the average temperatures with time. 

A similar lower L 6-plateau was .-'roposed by Lites, Shine 
and Chipman (1978). The microturbulent velocities adopted here 
are also lower (by about 20 percent) in the chromosphere than 
those of Vemazza, Avrett and Loeser (1973), since an extra 
broadening of the average profiles may be expected from the 
atmospheric motions. A summary of the initial and compute'' 
atmospheric parameter t is given in Table 1. 

Ill - HYDRODYNAMICAL COMPUTATIONS 
As described in Paper I, we have calculated the one- 
dimensional, non-linear dynamical state of this model atmosphere 
for various excitations and boundary conditions, as well as for 
different treatments o f the energy equation and ionization balance. 

We have included the effect of an isotropic, unresolved velocity field 
of "micro-turbulence" which is required for the radiative 
transfer, in the hydrodynamic calculations. Vertical variation 
of this velocity field, which remains constant in time, con- 
tributes to the mechanical support of the atmosDhere. Thus at a 
given (gas) pressure in the atmosphere, the density is lower than 
It would have been in the absence of the “micro-turbulence". 

To discuss the formation of optically thick chromispherlc lines, 
we have selected three time sequences illustrating : a) the 
propagation of a pulse through an initially motionless atmosphere, 
b) well established chromospheric and photospheric oscillations. 



RIGINAL PAGB B 
OF POOR QUALHY 

whose amplitudes correspond to those observed, and c) a similar 
set of oscillations, whose amplitude was reduced substantially 
to investigate the effect of non-linearities upon the hydro* 
dynamic and line formation results. The pulse introduced at 
the bottom of the atmosphere for experiment a) was one half 
cycle of a sinusoid in velocity, lasting 100 seconds and reach- 
ing a maximum velocity of 0.01 km/sec. That is to say the bottom 
boundary of the atmosphere was raised by 318 metres. Experiments 
b) and c) were excited by five cycles of a 300 second motion of 
the bottom boundary, whose amplitude was respectively 6.7 
and 1.05 m/r, corresponding to the total vertical motions of 
640 metres and 100 metres. 

Experiment a) might represent the impulse generated by a 
granule and is of additional interest here to compare our 
results with previous work (viz. Cram, 1976). The momentum transfer 
and heating of the chromosphere result in a very substantial 
difference between the mean state of the dynamical atmosphere 
and the initial atmosphere, so that this excitation is not 
appropriate for the study of the long term fluctuations of the 
line formation process. Experiments b) and c) result in a much 
more efficient transfer of energy to the resonant oscillations, 
so that for comparable resulting amplitudes in the oscillations 
negligible heating occurs and the mean atmosphere remains very 
similar to the initial atmosphere. Thus we have some hope that 
the resulting line profiles bear some resemblance to the 
observations. 





a simplified energy equation to expedite the calculations. 

Radiative transfer (and loss) are neglected in the internal 
energy of the gas, and ionization variations are considered 
to be snail fluctuations about the mean state. Locally, the 
state variables are related adiabatically, and the adiabatic 
exponent y (the ratio of specific heats) reflects the state of 
collisionally controlled ionization in the initial atmosphere. 

When the gas is completely ionized or completely neutral, y 
is 5/3. However, when an abundant species is partially ionized, 
changes in the internal energy of the gas primarily produce changes 
in the Ionization -because of the large ratio of the ionization 
energy to the thermal energy- and thus the temperature fluctuations 
are substantially reduced, y is reduced to 1.1 in the 
hydrogen convection zone. Above the temperature minimum, where the 
ionization is maintained by radiative, rather than collisional, 
processes which are slow compared to the dynamical time scales 
(Kneer and Nakagawa , 1976), we hold the ionization constant. 

The dynamical aspects of these experiments have been 
discussed in Paper I. He present in Figures la, b and c the 
variatior of the velocity ami temperature for the intervals 
during which we have investigated the line profiles. Two 
essential points should be recalled. First, during the propag- 
ation of energy, velocity and pressure (and hence temperature) 
vary in phase, as is seen in experiment A. However, they vary 



In quadrature for both trapped and evanescent oscillations. This 
difference of relative phase has a significant effect on the 
spectral diagnostics. Second, since for established oscillations 
the chromospheric velocity field Is primarily non -propagating, 
the chromosphere moves more or less in phase, without important 
vertical velocity gradients, so that even though the particle 
velocities are substantial fractions of the sound speed 
mechanical dissipation Is not significant. In addition, the 
absence of large vertical gradients reduces the importance of 
averaging along the line of sight Inherent in spectral line 
formation, so that there is good reason to suppose a priori 
that line shifts will reflect the chromospheric velocities. 

IV - COMPUTATION OF LINE PROFILES 

Let us now study the temporal variations of the Mg II k 
and Ca II K profiles corresponding to the atmospheric motions 
described above. For that purpose, we compute temporal sequences 
of profiles corresponding to the three hydrodynamical experiments. 
In every case, the moving atmospheres are sampled every 10 seconds. 
For experiment a), the computations are restricted to the interval 
from 200 to 400 seconds, when the primary pulse travels through 
the photosphere and chromosphere. For the "standing oscillations" 
experiments b) and c), we concentrate our attention on the inter- 
val from 1600 to 2600 seconds, where a sequence of typical 
chromospheric oscillations occurs. 

Because of the short time constants associated with the de- 
excitation of atomic levels, especially for the resonance lines. 



compared to the dynamical time scales, we use a quasi -static 
approximation to compute the line profiles, to the other 
hand, the photoionization rates depend mainly on the continue, 
and principally on the hydrogen Lyman continuum, which are 
expected to vary slowly, as seen in the precedi no section. 

He thus take these intensities as constant, and equal to those 
of the initial atmosphere . 

The computation of line source functions and profiles 
requires a simultaneous solution of the radiative transfer 
equations for the lines and statistical equilibrium equations 
for the atomic level populations. The radiative transfer 
equations are expressed in lagrangian coordinates, as in our 
previous computations (Gouttebroze, 1977, hereafter G 77). Methods 
using this "comoving frame" formulation have been described by 
Noerdlinger ami Rybicki (1974) and Mihalas et al (1976). The 
large numbers of profiles corresponding to long time sequences, 
e.g. 100 profiles both for Mg II and Ca II in experiment b), 
lead us to make some drastic simplications in order to keep the 
computer time within reasonable limits. He make three principal 
simplications : first, we reduce the frequency mesh, for 
intensity integration to 25 points instead of the 39 used in 
G 77 (for moving atmospheres, we must use symmetric frequency 
meshes, so that these numbers correspond for a static atmosphere 
to 13 and 20 points, respectively). This reduction is the main 
source of uncertainty, but it allows a significant reduction 
of the computer time used in matrix inversions. He use the 
Eddington approximation for the Integrations over direction. 



Secondly, we use average net radiative rates for transitions 
other than the one computed . Me first calculate the complete 
set cf line: corresponding to the average atmosphere, self- 
consistently. He then compute the time sequence of profiles 
for the transition under study, assuming that the net radiative 
rates in other transitions remain equal to the average- 
atmosphere rates. The computing time for an instantaneous 
profiles is thus equal to that required by a 2-level atom 
calculation. This change has negligible effects in the case 
of Mg II, the k line formation process being close to that of 
a 2-level atom transition. The effect on the Ca II lines is 
somewhat more important (although not essential), this 
simplification being equivalent to neglect of the Influence 
on the K line of variations in the Ca II infrared triplet. 
Finally, we replace the exact redistribution (for the 
scattered fraction which is coherent in the atom’s frame) 
by the Rjj N function of Argyros and Mugglestone (1971), which 
requires only the computation of an exponential and a Voigt 
function. This procedure is particularly efficient when a fast 
Voigt function routine is used. Our numerical experiments on 
the Mg II k line have shown that this change accounts for an 
error of a few percent in the emergent intensities, much less 
than from the first approximation. 

For Mg II, we use a 3-level and continuum model atom 
described in G 77. The Ca II model atom has 5 levels and a 
continuum, and includes the H and K resonance lines, and the 
infrarea triplet. The atomic parameters are taken from Shine 



and Linsky (1974), except for the line broadening parameters 
for which we use the empirical values of Ayres (1977). 

Before looking at the results of these computations, it 
may be useful to bear in mind some simple considerations about 
the effects of the variations of the physical parameters 
(velocity, temperature and density) on the lines. The velocity 
acts by shifting the absorption (and emission) profile by the 
Doppler effect. The relation between the frequency of a spectro- 
scopic feature and the velocity at a given lagrangian al'. tude 
is thus straightforward when this feature is formed in a limited 
part of the atmosphere (which is the case for optically thin 
lines or, as shown below, for the central reversals of the 
Ca II and Mg II lines). The temperature and the density vary 
generally in phase and both contribute to increase the intensity. 
For a two-level atom (a situation not too far from that of our 
resonance lines), the intensities depend upon two parameters : 
first, the Planck function B which, in the ultraviolet and at 
chromospheric temperatures, is given by 

n \ 2h\l , y . / hv . 
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(the symbols have their usual meaning) ; 
secondly, the ratio e of collisional to radiative de-excitation 
rates : 


c 


% 




which depends mainly on the electron density (and thus the 
density p of the gas) provided that C 21 does not vary rapidly 
with temperature. As limiting cases, the emission in an optically 
thin line depends linearly on cB, the photon creation term, 
while the intensity goes to B when one gets closer to L.T.E. 

An interesting intermediate case is that of a semi-infinite 
isothermal atmosphere (whose temperature and density vary in 
phase everywhere) : the emergent intensity at line center then 
varies as Bt'T. 

To evaluate the relative contributions of density and 
temperature fluctuations to the intensity fluctuation, one 
may consider an elementary isentropic pressure increase dp, for 
which density and the temperature increase such that : 

dNe _ d_£ = _1 dp and dT = y - 1 d£ 
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The intensity variation is : 
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where n is 0, 1/2 or 1 depending on the problem under consider- 
ation as discussed above. 
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In chromospheric conditions, the effects of the temperature 
variations on intensities are thus generally far more important 
than those arising from density changes. 

It should be emphasized that velocity and temperature 
changes do not act independently. For instance, an increase of 
temperature in an internal slab produces a more important increase 
of intensity when it is coupled with a Doppler shift, because 
the photons created escape more easily from the atmosphere. 

Such an effect occurs during the passage of a pulse, as described 
in the next section. 

V - VARIATION OF LINE PROFILES ASSOCIATED WITH THE PROPAGATION 

OF A SINGLE PULSE 

The temporal sequences of profiles obtained between 200 
and 390 seconds are shown in Figure 2a and b for the Mg II k 
and the Ca II K line, respectively. The behavior of the two 
lines is roughly the same : the profiles being initially 
symmetrical, the blue peak begins to increase with respect to 
the red pe k (around 250 s for Ca II, and 280 s for Mg II), 
reaches a maximum (around 300s for Ca II, 330 s for Mg II), 



after which the asymmetry reverses (around 330 s for Ca II, 

360 s for Mg II). In the Ca II sequence, a new "blue asymmetry" 

(i.e. with a ratio I(K2V) / I(K2R) larger tnan 1) appears 
about 380 s, due to the passage of the secondary pulse. These 
changes may be roughly interpreted as follows (with the help 
of the "contribution functions" of Figure 5) : At rest, the 
k2 peaks are formed principally in the low chromosphere (around 
1000 km for Mg II k). The pulse travelling upwards through the 
atmosphere produces an increase of temperature and density 
(and thus photon creation) coupled with an upward velocity. 

When the pulse has reached the formation region of k2 (at rest), 
it begins to produce an emission of photons shifted to the 
blue. Only the most blue-shifted photons then escaDe through 
the motionless overlying atmosphere, with its centered absorption 
profile, which produces an increase of the blue peak intensity. 

On the contrary, the red peak continues to be formed principally 
in the low chromosphere with small velocities and remains almost 
steady, which produces a blue asymmetry. Before the pulse reaches 
the k3 formation region, the central intensity begins to increase due to 
the incoherent scattering within the line core of photons created by tne 
passage of the pulse in underlying regions. When the pulse 
reaches the "k3 slab", the wavelength of the reversal is shifted 
to the blue with respect to the mean position of the peaks, which 
results in a red asymmetry. These results are consistent with 
those that Cram (1976) obtained for the Ca II K line. The delay 
of about 30 seconds between the events in the Ca II and Mg II 



lines corresponds to the time for the pulse to travel, from the 
formation altitude of a part of the Ca II profile, to the 
formation altitude of the corresponding part of the Mg II 
profile. 

In order to analyse these variations of profiles more 
closely, we tried to determine more accurately the positions 
and intensities of the peaks and central absorption of the 
lines by fitting parabolae to these three features. The 
variations of the three intensities for the Mg II k line are 
shown in Figure 3. The intensity increase resulting from the 
passage of the pulse, the delay between changes in k2V and 
k3, and the variation of asymmetries appear clearly. 

Having simultaneously the values of the physical para- 
meters (particularly the velocities, temperatures and electron 
density) and the positions and intensities in the lines, we 
investigated the correlations between these two kinds of para- 
meters, in order to establish the diagnostic possibilities of 
these "optical ly-thick" transitions. The correlation between 
the shift of k3 and the velocity of the slab where the optical 
depth at line center is equal to 1 in the average atmosphere 
is remarkably good for Mg II (Figure 4a), and not so good but 
nevertheless reliable, for Ca II (Figure 4b). The central 
reversal of these optically thick lines is thus a possible 
diagnostic tool to determine the vertical velocities in some 
parts of the chromosphere. On the contrary, we failed in our 
attempt to find a similar relation between the position of the 
k2 peaks and the velocity at a definite lagrangian altitude. 



The reason appears more clearly when one looks at the contribut- 
ion functions. 

The contribution functions used here represent on- 
tribution at geometrical depth z to the vertical *■ erging 
intensity, per unit of z : 

c(z,v) = S(z,v) k ( z , v) exp I- t ( z , v) | 

(S : source function ; v : frequency ; < : absorption coefficient, 
t : optical depth). These functions for the central reversal and 
the peaks are shown in Figures 5a and b, respectively, for the 
Mg II k ’ine. The contribution functions for k3 are sharply peaked 
at a definite and rather constant lagrangian altitude, i.e. 
the main part of the emission comes from a thin slab of the 
chromosphere, which is always about the same. For Ca II K3, 
these functions (not represented here, but similar to those of 
Figure 12e related to oscillations) cover a wider, though limited 
(about 500 km), range of depth, so that the Doppler shift of the 
reversal is the result of an average of the velocities across 
this region. For the peaks, on the contrary, the contribution 
functions are spread over regions as thick as 1000 to 1500 km 
(for Mg II as well as for Ca II), covering the main part of 
the chromosphere, where the density varies by several orders 
cf magnitude. This prevents us from using the displacements 
of the peaks as a practical tool to determine chromospheric 
velocities. Nevertheless, Figure 5b shows an interesting peculiar- 
ity of the blue peak : its contribution function exhibits a sharp 
peak which follows, during nearly 100 seconds, the passage of 



the pulse. This might have applications in the diagnosis of 
such pulses from spectroscopic data. 

The phase difference between velocity and temperature 
fluctuations represents an important characteristic that 
distinguishes between different types of waves. As shown 
above, the velocity in the k3 slab is closely related to the 
Doppler shift of k3. Unfortunately, the determination of the 
temperatures from the intensities is not so straightforward : 
due to the scattering effect mentioned above the intensity in 
k3 depends on the temperatures and densities in a large 
region of the atmosphere with a thickness corresponding 
approximately to the thermal izati on length. Thus, the 
temperature-intensity correlation is accurate only when the 
temperatures arr approximately in phase throughout the chro- 
mosphere, as will be the case in the next section. In t 1 - 
case of a pulse propagation, the temperature and velo'i w 
variations are approximately in phase, but the k3 intensity 
begins to increase about 40 or 50 seconds before than the pulse 
has reached the corresponding slab. This effect may be seen in 
Figure 6, where we plot the reduced k3 intensity : 
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and the corresponding reduced wavelength and temperature. 



VI - VARIATIONS OF LINE PROFILES ASSOCIATED WITH NON-PROPAGATING 

OSCILLATIONS 

Let us now examine the variations of the Ca II and Mg II 
resonance line profiles emerging from the freely oscillatino 
atmospheres defined by the experiments a) and c) of § III. The 
profiles were calculated between 1600 and 2590 seconds, and thus 
begin 100 seconds after the end of the excitation phase. Figure 
7a shows the whole sequence of Mg II profiles corresponding to 
the large amplitude oscillations. The intensities are varying 
with a quasi-periodicity of about 200 seconds, as expected; the 
maximum corresponding to the minimum o f the altitude (the 
p *essure and the temperature of the gas then reach their 
maxima). At this minimum altitude point, the velocity changes 
quickly, the gas moving downwards before beginning to move 
upwards. For such general, one-dimensional motions, the 
decrease of density wit altitude produces an increase of the 
velocity amplitude. If the gas moves downwards, the formation 
region of k3 is shifted to the red with respect to the under- 
lying emitting layers, which results in a blue asymmetry. 
Conversely, a similar upward motion produces a red asymmetry. 
Changes of a blue to a red asymmetry associated with maxima 
of intensity may be seen on Figure 7a around sequence numbers 
169, 191, 209, 229 and 248. We obtain similar results with 
other experiments. Figures 7b, c and d show the last 20 
sequences, corresponding approximately to a cycle, for large 
amplitude oscillations (Mg II and Ca II), and smell oscillations 



(Mg II). 7b and c show that the oscillations are practically 
in phase in Mg II and Ca II, the change of asymmetry occurring 
between sequence numbers 248 and 249 in each case. This con- 
trasts with the result obtained for a pulse propagation. 

Another feature visible on Figure 7c is the presence of profiles 
with a single peak (sequences 243 to 245) which are associated 
with large velocity gradients and relatively low chromospheric 
emission. Profiles of this kind occur frequently on the Sun. 

One may notice that both singly-peaked and doubly-peaked 
profiles are explained by the same theory of line formation 
without difficulty. (Nevertheless, we must recognize that 
singly-peaked profiles are less frequent in our computations 
than in the observations ; this might certainly be modified 
through adjustments in the initial model atmosphere). 

Applying the same method of least square fit to the peaks 
and reversal as in the preceding section, we obtain the variations 
of intensities shown on Figure 8a and b for the large and small 
amplitude oscillations, respectively. In the case of large ampli- 
tude, the variations of intensity are sharply peaked, with 
positive excursions reaching an order of magnitude, while the 
variations remain almost sinusoidal in the other case, (One 
should notice that the variation of k3 is in phase with that 
of the average of the two peaks). 

As seen above, an important property of non-propagating 
oscillations is that velocities and temperatures vary in quadra- 
ture. As the temperatures are approximately in phase throughout 
the chromosphere, and thus the intensities, we obtain a similar 



phase difference between the intensity and the Doppler shift 
of k3. This is shown in Figure 9 which represents these two 
quantities in reduced form : one may see that the instants of 
maximum acceleration correspond approximately to the extrema 
of uhe intensity. 

In data analyses, it is a current practice to use the Fourier 
transform for the determination of phase shifts. The results of 
such an analysis for the variations of shift and intensity in k3, 
are summarized in Figure 10 : the power spectra of the shift and 
the intensity have a peak at 5 mHz, as expected, and several 
harmonics (the harmonic at 10 mHz is particularly strong in the 
case of the intensity). The phase difference between the components 
of the Doppler shift and intensity is close to 9Q C for the 
frequencies corresponding to the fundamental and the first 
harmonic, in accordance with the direct comparison of the oara- 
meters in Figure 9. 

In Figure 11, we compare the Mg II k profile emitted by 
a "mean atmosphere" -obtained by averaging the physical para- 
meters (temperature, density, velocity) over the 100 sequences- 
with the average of the 100 profiles of the sequence. The first 
aim of this comparison is to investigate the effects of non- 
linear averaging of the intensities which is expected to be 
important for the ultraviolet spectrum. This effect appears 
clearly for the average intensity which is about 25% higher 
in the peaks than that emitted by the average atmosphere. 

Another aim is to try to explain discrepancies between the 
observed profiles and those computed from static atmospheres. 



As a general feature, observed Mg II k profiles have broader 
peaks and less deep central reversals than the calculated 
ones. From that point of view, the mot ion -averaged profile 
of figure 11 is far more satisfying than the static profile. 

Another puzzling feature is the average blue asymmetry of 

the observed profiles. Although our motion-averaged 

profile has a blue peak slightly higher than the red 

one, the difference is small and does not compare well 

with the observations. In our one-dimensional scheme, the 

addition of a mean downflow would be required to explain the 

observed average asymmetry. Nevertheless, many effects, 

arising in particular from multi -dimensional hydrodynamics 

and radiative transfer, remain to be investigated before 

concluding in favour of the reality of such a downflow. The 

formation process of the peaks and reversal appears when one 

examines the contribution functions, defined as in § V 

( r igure 12). The atmosphere undergoes a series of short, 

high temperature, contraction phases and long, low temperature 

exapansion phases. The central reversals are formed, approximately, 

at a constant lagrangian altitude (Figures 12d and e) with 

sharp intensity maxima in the neighborhood of the minimum 

altitude times. As was the case for the pulse propagation 

studied above, the regions of formation of the k2 peaks move 

over a range of lagrangian altitudes of about 1000 km. During 

the "cool" expansion phase, the maximum contribution for both 

peaks is located in the low chromosphere (around 1000 km for 

Mg II K) (Figure 12b). As one gets closer to the maximum contraction 

phase, the region of formation of the blue peak moves upwards with 



increasing intensities (Figure 12a). At t'is point, the intensity 
in the blue peak decreases rapidly and a red peak appears in 
the high chromosphere, close to the formation region of the 
k3 reversal (Figure 12c). Then, the region of formation of 
this red peak moves downwards with decreasing intensities, 
wh' h is the beginning of a new expansion phase. As in the pulse 
propagation experiment (§V), the central reversal is formed in 
a narrow range of quasi constant lagrangian altitudes and thus 
produces a good correlation between the Doppler shift of these 
reversals and the velocity in the corresponding layers. This is 
valid for Ca II as well as for Mg II, which gives a diagnostic 
tool to determine the velocities at two different altitudes 
and thus discriminate propagating from non-propagating waves. 

As seen in § V, the k3 intensity diagnostic !* not so easy to 
use : although some general correlation exists between temperature 
and intensity, the latter is the resi.it of a non-local transfer 
process and thus depends on both temperatures and densities 
throughout a wide range of altitudes. The fact that the oscillatory 
motions are in phase throughout this range gives rise to a good 
correlation, which however is of very limited diagnostic value. 

VII - CONCLUSION 

We have constructed a time-dependent, one -dimensional model 


of the solar chromosphere satisfying radiative transfer (line 
profiles) and hydrodynami cal (oscillation periods in the range 
150 to 250 seconds) requirements. The aim of computing Ca II 
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TABLE 1 


INITIAL ATMOSPHERIC PARAMETERS 


m 

T 

V T 

(g/cm 2 ) 

(K) 

( km/s ) 

7.90 (-6) 

1 000 000 

8.0 

7.95 (-6) 

1 000 000 

7.9 

8.05 (-6) 

316 000 

7.7 

8.08 (-6) 

158 000 

7.6 

8.10 (-6) 

79 400 

7.6 

8.12 (-6) 

39 800 

7.5 

8.18 (-6) 

20 000 

7.4 

8.40 (-6) 

19 000 

7.0 

8.80 (-6) 

17 800 

6.2 

8.85 (-6) 

12 600 

6.1 

8.90 (-6) 

9 000 

6.0 

1.00 (-5) 

7 410 

5.9 

1.26 (-5) 

6 790 

5.7 

1.58 (-5) 

6 640 

5.5 

2.51 (-5) 

6 520 

5.1 

3.98 (-5) 

6 430 

4.7 

1.00 (-4) 

6 310 

4.0 

2.51 (-4) 

6 190 

3.3 

6.31 (-4) 

6 030 

2.5 

1.58 (-3) 

5 780 

1.8 

3.98 (-3) 

5 370 

1.0 

1.00 (-2) 

4 920 

1.0 

2.51 (-2) 

4 570 

1.0 

5.01 (-2) 

4 510 

1.0 


h N H N e /N H 

( km ) (cnf 2 ) 

2 745 6.83 (+8) 1.20 

2 438 6.88 (+8) 1.20 

2 076 2.19 (+9) 1.20 

2 035 4.37 (+9) 1.20 

2 021 8.60 (+9) 1.20 

2 014 1.68 (+10) 1.20 

2 003 3.42 (+10) 1.05 

1 977 3.73 (+10) 1.04 

1 935 4.28 (+10) 1.02 

1 931 6.18 (+10) 9.12 (-1) 

1 928 9.04 (+10) 7.37 (-1) 

1 884 1.26 (+11) 6.41 (-1) 

1 812 1.83 (+11) 5.14 (-1) 

1 749 2.50 (+11) 4.03 (-1) 

1 635 4.60 (+11) 2.23 (-1) 

1 536 8.15 (+11) 1.19 (-1) 

1 359 2.33 (+12) 3.76 (-2) 

1 200 6.39 (+12) 1.45 (-2) 

1 054 1.74 (+13) 6.92 (-3) 

920 4.73 (+13) 2.93 (-3) 

797 1.32 (+14) 8.65 (-4) 

686 3.63 (+14) 2.57 (- 4 ) 

582 9.78 (+14) 1.48 (- 4 ) 

508 1.98 (+15) 1.32 (- 4 ) 



1.00 (- 1 ) 

4 500 

1.0 

435 

3.96 (+ 15 ) 

1.20 

(- 4 ) 

2.00 (- 1 ) 

4 570 

1.0 

361 

7.80 (+ 15 ) 

1.14 

(- 4 ) 

3.98 (- 1 ) 

4 740 

1.0 

285 

1.50 (+ 16 ) 

1.17 

(-*) 

1.00 

5 040 

1.0 

179 

3,54 (+ 16 ) 

1.26 

(- 4 ) 

2.51 

5 570 

1.0 

64 

8.05 (+ 16 ) 

1.70 

(- 4 ) 

3.98 

6 170 

1.0 

0 

1.15 (+ 17 ) 

3.63 

(- 4 ) 

5.01 

6 920 

1.0 

35 

1.30 (+ 17 ) 

1.24 

(- 3 ) 

6.31 

8 320 

1.0 

77 

1.35 (+ 17 ) 

8.91 

(- 3 ) 

7.94 

9 550 

1.0 

- 125 

1.45 (+ 17 ) 

3.17 

(- 2 ) 

1.00 (+ 1 ) 

10 300 

1.0 

- 181 

1.66 (+ 17 ) 

5.65 

(- 2 ) 

1.58 (+ 1 ) 

11 500 

1.0 

- 306 

2.24 (+ 17 ) 

1.14 

(- 1 ) 

2.51 (+ 1 ) 

12 700 

1.0 

- 455 

3.03 (+ 17 ) 

1.92 

(-D 

3.98 (+ 1 ) 

13 700 

1.0 

- 626 

4.24 (+ 17 ) 

2.59 

(-D 

6.31 (+ 1 ) 

14 700 

1.0 

- 819 

5.98 (+ 17 ) 

3.25 

(-1) 

1.00 (+ 2 ) 

15 800 

1,0 

- 1 036 

8.40 (+ 17 ) 

3.97 

(-1) 

1.58 (+ 2 ) 

16 900 

1.0 

- 1 278 

1.19 (+ 18 ) 

4.59 

(-D 

2.51 (+ 2 ) 

18 100 

1.0 

- 1 550 

1.70 (+ 18 ) 

5.19 

(-D 

3.98 (+ 2 ) 

19 400 

1.0 

- 1 851 

2.44 (+ 18 ) 

5.75 

(-D 



FIGURE CAPTIONS 


Figure 1 

Normalized velocities and temperatures as functions of 
time (seconds) for 15 selected lagrangian altitudes (kilo- 
meters) indicated at the left of the figure. The extrema of 
the velocities are written in the right side. 

a) propagation of a pulse ; 

b) large amplitude oscillations ; 

c) small amplitude oscillations. 

Figure 2 

Time sequence of line profiles, computed every 10 seconds, 
during the passage of a pulse through the atmosphere. The 
sequence n° 21 corresponds to t = 200 seconds. 

a) Mg II k line ; 

b) Ca II K line. 

(in abscissa : relative wavelength in angstroms). 

Figure 3 

Time variations of intensities in the Mg II k line 

-1 -2 -1 -1 

(time in seconds, intensity in erg s cm sr Hz ). 

Square symbols : central reversal (k3). Triangles : 
blue peak (k2V). Crosses : red peak (k2R). 

Figure 4 

Comparison of time variations of the position of the 
central reversal, converted into velocity (squares) with 
the velocity of the slab corresponding to an optical 
depth unity at line center (velocities in km/s). 

a) Mg II k ; 

b) Ca II K. 



Figure 5 

Mg II contribution functions versus eulerian altitude 
(in km) for 20 time sequences (pulse propagation). 

a) k3 (the dots correspond to constant lagrangian altitudes) 

b) k2 (solid line : blue peak, dotted line : red peak). 

The contribution functions are shown to \% of their peak 
value. 

Figure 6 

Time variations of reduced Mg II k3, intensity (squares), 
Doppler shift (triangles) and reduced temperature 
(crosses) in the k3 formation slab (pulse propagation). 

Figure 7 

Sequences of line profiles during chromospheric oscillations, 
separated by 10 seconds. Sequence n° 161 corresponds to 
t * 1600 seconds. 

a) and b) large amplitude oscillations and Mg II k line ; 

c) large amplitude oscillations and Ca II K line. 

Figure 8 

Time variation of Mg II k line intensities in the k3 
reversal (squares), k2V (triangles) and k2R (crosses) 
peaks . 

a) large amplitude oscillations ; 

b) small amplitude oscci nations. 

Figure 9 

Reduced intensity (squares) and wavelength (triangles) 
versus time for the Mg II k3 reversal (large amplitude 
oscillations) . 



Figure 10 


Fourier transform of the time variations of Doppler shifts 
and intensity of Mg II k3 (large amplitude oscillations), 
Squares : power spectrum for Doppler shifts (arbitrary 
units) ; triangles : power spectrum of intensity ; 
crosses : phase difference between Doppler shift and 
intensity components (degrees). 

Figure 11 

Comparison of time-averaged profi the Mg II k line 

(squares) with the profile corresponding to the average 

-1 -1 -1 

atmosphere (triangles) (intensities in erg s cm sr Hz , 
wavelength in anstroms) (large amplitude oscillations). 

Figure 12 

Time sequences of contribution functions versus eulerian 
altitude in kilometers (large amplitude oscillations). 

The dots indicate constant lagrangian altitudes. 

a) Mg II k2V ; 

b) same as a) normalized at each time to show more 
clearly the low chromospheric contribution which dominates 
the "quiet" profile ; 

c) Mg II k2R ; 

d) Mg II k3 ; 

e) Ca II K3. 
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ORIGINAL RAGE uS 
OF I’OOR QUALITY 


c 'rgy to the corona away from active regions and its dissipation 
has become increasingly attractive. While it is tempting to consider 
coronal models that are scaled down aralogs of flate models where e r. e r c \ 
is released gradually, ve should recall that it has been the t equi r r.t 
of rapid ewergy release that has driven the search for flare mechanist.-, 
and that processes which are too slow to describe flares, and c^r. 

discarded, may be relevant for the quasi-steady nesting of the ctr.r.a. 

Although we are not discussing generation ir.ee nan isms here, w<. sn 
point out that the observed fields are roughly in equipart i tion w;t. :..e 
kinetic energy density at the top of the convection cone \JVcKO.<,r t-r.u 
Weiss, 1S76). The seed field for amplification i : tne convec t:v t -^tions 
presumably results from the dynamo interaction between rotation .?nu 
convection in deeper layers. Thus, the magnetic field associated 
ing will vary with the spin down time of a star, which in general will 
be unrelated to the evolutionary time scale determined by nuclec ryr.n.c.- 

In the convection zone : gas pressure /’macne tic pressure > 1 

and convective velocities move the magnetic line- of :\-ra wit*. li:t 
resistance. Magnetic fields provide (a. leas: > t.urec broad meaner it 
for transmitting energv from the convention z o r. e t o tne corona* 
waves propagating along lines of force, 2) current.- inched in : ■; 
corona by wa**es, twists and motions of the magnetic mir.t - . and 

3) emerging : 1 ux carrying stored energv upwrrus. Iiss:ta:*cr. c: 
energy occurs primarily as tne uissipa: :or. of inched current.-. Alt 
where the magnetic energy j& in t for;., of Al;v c n waces, visco^t -a 
and wave mode coupling also ocoot . 

Alfven Waves 

The field may serve (moi or less passive!*; as ^ medium r trans- 
mitting the shaking or twisting of field lines (Wer.tzel, 197c; l. 
and Kaburaki, 1974; Hollweg, i960). Twisting cr jiggling of a m^unet^ 
line of torce can be produced directly by convective motions or r- 
thermal over-stability in the convectively unstable layer. 7i.c 
ance will then propagate away as a wave, or if slow enough, as a q-asi- 
static change in the i.cld. Several magnet o-acou* t i : -zravi ty w. s ^ 
modes exist. One can understand the essence of magnetic waves by 
sidering only Alfven waves. Their restoring force is magnetic tension, 
so they act like waves on a string — the magnetic field line — wit. roti.n 
transverse to the field. The temperature, aensi:>, pressure and r.ig;ut;c 
field strength (§ 0 + {'&) are constant. Energy nropagate^ along :..c field 
at the Alfven speed 



Because of its density dependence, the Alfven speed increases hr 
a factor of 10^ between the photosphere and the corona, so that - 

stantial reflection of Alfven waves will occur. Since the magnet, 
field in the coron is inhomogeneous, the internal /lfven waves will 
be slightly modified and propagate as surface waves wi*h the average 
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ABSTRACT 


Hie outer atmospheres of stars must 

be heated by some non- therm 


energy flux to produce chromospheres and coronae. We discuss processes 
which convert the non-thermal energy flux of organized, macroscopic 
motions into random, microscopic (thermal) motions. Recent advances in 
our description of the chromosphere velocity field suggest that the 
acoustic waves observed there transmit very little energy, and hence 
are probably incapable of heating the upper chromosphere and corona. 

The apparent failure of this long held mechanism and the growing appre- 
ciation of the importance of strong magnetic fields in the chromosphere 
and corona have led to hypotheses of heating by the dissipation of cur- 
rants (both oscillatory and quasi-steady). This follows discoveries 
in laboratory and ionospheric plasmas and work on solar flares, that 
instabilities can concentrate currents into thin high current density 
filaments where they dissipate rapidly. 


INTRODUCTION 

The major advances in understanding the Sun and other stars 
the second half of this century have resulted from our appreciaticr. of 
violations of the paradigm: stellar structure (both interior and atmos- 

pheric) is determined by thermal energy transport (i.e., radiation, 
convection and to a much lesser extent conduction). It now appears 
that the outer atmospheres of virtually all stars are substantially 
hotter than radiative equilibrium would have predicted and that stellar 
winds are important energy losses for the oute- atmosphere. The con- 
ventional wisdom has become: a thermodynamic engine converts internal, 

thermal energy (maintained by thermal energy transport) to external, 
macroscopic energy which then propagates, or is transported, with 
little interaction with the material through which it propagates, up- 
wards into the corona. Generally, production of propagating energy is 
not very efficient, €10~ . However, where the bulk of the total out- 
ward energy flux escapes from the star, i.e., within the outer several 



photon mean free paths * the macroscopic energy now represents a sub- 
stantial proportion of the energy density, and dissipation of this 
energy can substantially change the temperature structure. In addition, 
the mechanical energy density directly contributes to the momentum bal- 
ance, and hence structure, of the overlying atmosphere. Dissipation of 
the macroscopic energy results from microscopic processes (viscosity, 
conductivity and resistivity) which reduce gradients and increase rancor 
internal motions; i.e., increase entropy. By dissipation, we understand 
the ultimate transformation of macroscopic, organized energy to micro- 
scopic disorganized ones, and we don't discuss the transformation of ont 
propagation mode to another, which "dissipates" the mode, but not the 
non- thermal energy. 

We shall discuss the wave modes (acoustic, internal gravity and 
Alfven) which have been considered, and the role they appear to play in 
the Sun, as well as the less fully studied hypothesis of heating by 
current dissipation. Ve do not discuss the generation of this propaga- 
tion energv here, and the reader is referred to Stein and Leibacher 
(1979, 1980). 

Finally, we should point out that the identification of the various 
candidates for propagating the required non-thernal energy flux pro- 
ceeded from observations of the time dependent, horizontally inhomogen- 
eous structure of the Sun. The Sun is where we learn the physics to 
understand the stars. On the other hand, the variation of non-thernal 
heating with stellar type provides a stringent test for all generation, 
propagation and dissipation theories. We should note several recent 
reviews of heating mechanisms: Cram (1977), Cram and Ulmschneidcr 

(1978), Rosner, Tucker and Vaiana (1978) — particularly the Appendix — 
and Wentzel (1978) . 

ACOUSTIC (PRESSURE) WAVES 

Propagating pressure fluctuations were the first non-radiative flux 
suggested to account for the corona's high temperature. These steepen 
to form discontinuities due to the dependence of their propagation speed 
on temperature: the upwardly moving, hotter part of the wave propagates 

faster than the downwardly moving cooler part. The steepening is accel- 
erated in a stellar atmosphere, since the wave amplitude increases as 
the density decreases, to keep the energy flux constant in the absence 
of dissipation and refraction. When the steepening progresses to the 
point that the wave changes appreciably over a collision mean-f ree-path, 
diffusion (viscosity and conductivity) acts against further steepening, 
converting the organized motion of the wave into thermal motions, and 
a shock wave is said to exist. In stellar coronae, the high thermal 
conductivity of the hot plasma is generally more important than viscosity 
As the pressure jump across the wave increases, the dissipation (or 
entropy increase) grows very rapidly: initially as the third power of 

the pressure jump. It is important to note that it is the sharp gradient 
or discontinuity, produced by the wave which produces the dissipation; 



^ pooh quaut? 


the existence of large velocities (comparable with the sound speed) is 
not sufficient. In particular, long period vertical motions will raise 
and lower the atmosphere more or less in phase, with the velocity vary- 
ing exponentially — not sinusoidally with height — and even at large amp- 
litudes "phase" is not propagated vertically. Radiative exchange (losses 
from hot, compressed gas and gains by the cool, rarefied gas) also will 
dissipate the flux of an acoustic wave. 

Short period waves (10 - 100 seconds), created by turbulent notions 
in the convection zone, should be capable of providing the heating re- 
quired near and immediately above the temperature minimum (Ulmschneider , 
Schmitz, Kalkofen and Bohn 1978, and references therein). However, 
they rapidly dissipate their energy and do not appear capable of any 
heating higher up. These waves lose a very substantial part of their 
energy by radiative damping in the photosphere before forming shock 
waves. Thus, the flux available for heating is a small number which 
is the difference of two large numbers (the generation and radiative 
damping), a situation which amplifies numerical as well as physical un- 
certainties. 

Although the dissipation of these waves is fairly well understood, 
the total power emitted and its temporal spectrum remains uncertain. 
Extrapolations to other stars have been suggested (Ren 2 ini, Cacciari, 
Ulmschneider and Schmitz. 1978), but the predicted, very rapid increase 
in the chromospheric emission with luminosity for sc.ar type stars does 
not appear to occur (Linskv and Ayres, 1979). Short period waves 
remain "unobservable" because their wavelengths are comparable to the 
widths of the regions emitting the spectral lines used as diagnostics, 
however they do contribute to the non-thermal width of spectral lines. 

Longer period acoustic waves are seen as spatially resolved dis- 
placements of photospheric and chromospheric lines . Theory and obser- 
vation have achieved rare agreement in the description of the "five 
minute oscillations" as acoustic waves trapped within certain ranges 
of temperature (Leibacher, 1977; Zirker, 1979). These standing waves 
transmit virtually no energy: neither in the cavity where upward and 

downward propagating waves, of equal flux. Interfere constructively, 
nor outside the reflecting boundaries of the cavity where the mode may 
continue to exist by "tunneling" as an evanescent wave. Thus, although 
M.e velocity amplitudes may be substantial, since pressure and velocity 
vary nearly in quadrature, little energy would be propagated. In addi- 
tion, since there is no vertical phase propagation, the large vertical 
gradients required for diffusive dissipation do not occur. 

The radiative losses of the upper chromosphere and corona appear 
to exceed the most optimistic estimates of the acoustic flux (Athay 
and White, 1976; Bruner, 1979): assuming that the total non-thernal 

width of chromospheric lines is due to acoustic waves propagating up- 
wards at the maximum possible energy propagation (group) velocity. 

For other stars, the situation is likely to be similar: short period 

waves lose a majority of their flux by radiation and dissipate the 



remainder in a relatively abort distance. Longer period waves develop 
substantial amplitudes only if they are trapped, and propagate little 
energy. 


INTERNAL GRAVITY (BUOYANCY) WAVES 

Buoyancy restoring forces, in a convectively stable atmospnere, 
will drive predominantly horizontal motions known as internal gravity 
waves (IGW). These waves are abundant in the earth's atmosphere, and 
will be generated by slow motions, such as penetrative convection. and 
velocity shear, such as that associated with the variation of the super- 
granulation velocity with altitude. Because of their short vertical 
wavelengths, they are not prominent in the observed macroscopic motions 
of the Sun, but they presumably make a significant contribution to the 
unresolved, "microscopic" velocity field. Dissipation through ordinary- 
viscosity am* conductivity will occur, and as they propagate upward? 
and their amplitude grows, the waves will "break" and produce shorter 
wavelength IGW turbulence along their fronts — like whitecaps. Large 
amplitude waves can drive mean flows, and conversely TJW's can be 
absorbed by a steady flow whose velocity equals the wave’s phase velocity. 
The wave energy goes into the mean flow, and it is not dissipatr- in the 
sense of returning to thermal energy. Radiative damping tends to make 
a wave isothermal, which destroys the buoyancy driving force, and Ler.ce 
is particularly effective in dissipating IGW's. 

The energy propagation speed for IGW’s is ± the sound speed, so 
the arguments of Athay and White (1978) for the energy flux carried if 
the excess line broadening were all taken as propagating mechanical 
energy, also apply for IGW's. However, because their velocity is pri- 
marily horizontal, only a small part of their amplitude contributes to 
line broadening at disk center. 


MAGNETIC FIELDS (CURRENTS) 

Recent observations have vastly increased our appreciation of the 
importance of magnetic fields on the Sun, and hence other stars. Strong 
fields have been found away from active regions: the super-granulation 

network consists of strong (greater titan 1000 gauss), concentrated (less 
than a magameter horizontal extent) fields, and flux emerges at all 
latitudes. X-ray observations indicate that the corona consists almost 
entirely of filamentary and loop structures which have been identified 
as magnetic flux tubes. The role of magnetic fields in heating solar 
and stellar flares has long been known, as has their importance in 
maintaining the enhanced emission of active regions. While the impor- 
tant physical processes are becoming clearer, we are still far from an 
accepted model for these phenomena (Spicer and Brown, 1980 ). Nonethe- 
less, as the difference between a spherically symmetric, one dimensional 
chromosphere or corona and reality becomes more obvious, the idea that 
the magnetic field plays an essential role in the transport of mechanical 
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c* irgy to the corona away from active regions and its dissipation there, 
has become increasingly attractive. While it is tempting to consider 
coronal models that are scaled down analogs of flare models where energy 
is released gradually, we should recall that it has been the requirement 
of rapid energy release that has driven the search for flare mechanisms, 
and that processes which are too slow to describe flares, and have been 
discarded, may be relevant for the quasi-steady heating of the corona. 

Although we are not discussing generation mechanisms here, we should 
point out that the observed fields are roughly in equipartition with the 
kinetic energy density at the top of the convection zone (Peckovt-r and 
Weiss, 1978). The seed field for amplification by the convective notions 
presumably results from the dynamo interaction between rotation and 
convection in deeper layers. Thus, the magnetic field associated heat- 
ing will vary with the spindown time of a star, which in general will 
be unrelated to the evolutionary time scale determined by nucleosynthesis 

In the convection zone i (* gas pressure/magnetic pressure) is large 
and convective velocities move the magnetic lines of force with little- 
resistance. Magnetic fields provide (at least) three broad mechanisms 
for transmitting energy from the convention zone to the corona: 1) 

waves propagating along lines of force, 2) currents induced in ti.e 
corona by waves, twists and motions of the magnetic foot points, and 
3) emerging flux carrying stored energy upwards. Dissipation of mu cr. otic 
energy occurs primarily as the dissipation of induced currents. Althoucn 
where the magnetic energy is in the form of Alfven waves, viscous damping 
and wave mode coupling also occur. 

Alfven Waves 

The field may serve (mot or less passively ) as « medium for trans- 
mitting the shaking or twisting of field lines (Wentzel, 1976; Uchida 
and Kaburaki, 1974; Hollweg, 1980). Twisting or jiggling of a magnetic 
line of force can be produced directly by convective motions or by 
thermal over-stability in the convectivelv unstable layer. Tne disturb- 
ance will then propagate away as a wave, or if slow enough, as a quasi- 
static change in the field. Several magneto-acoustic-gravitv wave 
modes exist. One can understand the essence of magnetic waves by con- 
sidering only Alfven waves. Their restoring force is magnetic tension, 
so they act like waves on a string — the magnetic field line-wit.- motion 
transverse to the field. The temperature, density, pressure and magnetic 
field strength (® 0 + d®) are constant. Energy Dropagate^ along Lite field 
at the Alfven speed 

a • B / V 6np . 

Because of its density dependence, the Alfven speed increases by 
a factor of 10^ between. the photosphere and the corona, so that sub- 
stantial reflection of Alfven waves will occur. Since the magnetic 
field in the corona is inhomogeneous, the internal /lfven waves will 
be slightly modified and propagate as surface waves wi *.h the average 
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Alfven speed. Where the local Alfven speed equals this average Alfven 
speed a resonance occurs , at which the transverse velocity and electric 
field become large. 

The Alfven speed is independent of Alfven wave amplitude , since 
the density and magnetic field remain constant. Hence. Alfven waves 
do not steepen and shock. However, they can dissipate by microscopic 
processes (viscosity and resistivity). The damping lengths for both 
viscous and Joule heating increase as the cube of th' magnetic field 
and the square of the wave period. 

Viscous damping is more important than Joule heating in the corona 
and vice versa in the photosphere and chromosphere. In the corona the 
problem is turned around. The damping lengths are so long that dissi- 
pation is very small and attention has focussed on coupling to ether 
modes which dissipate more rapidly. Surface Alfven waves have very 
large amplitudes at their resonant point, so most of their enerev will 
be dissipated by Joule and viscous dissipation in the th<r. resonant 
layer. The damping rate is determined by the race at which energy can 
propagate into the resonant layer. Mode coupling occurs when the 
Alfven wave amplitude becomes large, and a compressible, sound wave is 
formed which can dissipate efficiencly (Derby, 1976). Coupling alsr 
occurs in the presence of inhomogeneities, by producing ar. uncertain. v 
in the wave number, of the order of the reciprocal cf the inhomegeneity 
scale size. For two waves of the same frequency and wavenumbers differ- 
ing by Ik, the coupling will be of order 

l/Lj&k! (Melrose, 1974). 

In particular, Alfven waves propagating nearly along the field car. couple 
efficiently to fast mode waves. 

Current Dissipation 

Slow motions of the magnetic foot points in the convection zone 
. ly force flux tubes against one another in the corona, which will induce 
currents between them perpendicular to the field direction. If the mag- 
netic fields are oppositely directed, a neutral sheet will be produced. 
Buffeting of magnetic flux tubes in the photosphere by granule and 
supergranule motions will produce twists and Alfven waves which propa- 
gate up into the corona and induce currents there. Emerging magnetic- 
fields may also carry energy stored as twists Q x B / 0) up into the 
corona, where B is small. Any unbalanced magnetic force (J x p) ••ill 
cause the tube to move and relax (Parker, 1974). Field aligned urren's 
can remain in the relaxed coronal flux tube (Parker, 1974; Hollwcg, 19t>0). 

Current or magnetic field dissipation is a diffusive process due 
to single or collective particle collisions. The heating rate is 
Q « r»J 2 , where J is the current density and r, is the resistivity. The 
characteristic resistive diffusion time scale is -t R » 4nL 2 /nc 2 . 
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For typical coronal parametnre tr v 10* yrs; coo long to be significant, 
this Joule dissipation tine can be reduced either by reducing the width 
L of the region through which the currents flow or by increasing the 
resistivity hy increasing the effective collision rate. To get signi- 
ficant current dissipation by any mechanism, the dissipation must occur 
in such small volumes that the transfer of the resulting beat to the 
rest of the corona is a serious problem. 

Tearing Instability 

The size of the current region can be reduced by the tearing 
instability in a sheared magnetic field. Parallel currents attract one 
another end tend to clump. The dumping of current produces a fluid 
flow that forces the sheared magnetic field into X-tvpe neutral points 
Fi lamented currants are produced with small enough length scales so that 
the classical Coulomb collision resistive diffusion time becomes small 
and the magnetic field can tear and reconnect, and the currents can 
dissipate (Drake and Lee. 1977; Bateman. 1976, chapter 10). 

This mechanism has been invoked for the violent energy release in 
flares. (See Spicer and Brown. 1980.) However, short wavelength tear- 
ing modes release less energy than the long wavelength instability thought 
to occur in flares, and so cay produce a more tranquil heating appropri- 
ate for coronal flux tubes. The shorter wavelength modes distort the 
field lines more, which produces a greater restoring force, and also 
have a smaller volume of magnetic energy they can release. The tearing 
instability has a lower threshold than the current driven instabilities 
which lead to anomalous resistivity (discussed below) . 

Thermal Instability 

Thermal instabilities can also reduce the si ’e of the current region 
by producing current filamentation, because Coulomb resistivity decreases 
with increasing temperature. Hence, Joule heating will increase the 
temperature, whicn reduces the resistivity, which increases the current 
density, which leads to enhanced heating. Current filaments parallel 
to B result. This instability can act as a trigger for other instabili- 
ties. 

Anamolous Resistivity 

When the current density, J = n e e becomes large enough 

that the drift velocity approaches the electron thermal velocity, then 
substantial numbers of electrons will tend to run away and generate 
several different types of electrostatic waves. These waves bunch the 
ions, so that the electrons collide with the electric field of a large 
collective charge rather than that of a single ion. This scattering of 
electrons by the waves increases the effective collision rate, the rate 
of momentum transfer and hence the resistivity. The enhanced resistivity 
due to electron scattering by plasma waves is called "anomalous resist- 
ivity", and since it occurs in conjunction with current filamentation 
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will shorten tbe resistive diffusion ties trewBioMly (Fapmdopoious, 
1977; Iseatti folub, Ceppl end Vaiana , 1978). 

Short period acoustic waves (period < acoustic cut-off periods still 
appear to he the most satisfactory mechanism for heating the low ci.ror.o- 
Sphere, these waves shock before reaching the corona, so increasing 
their flux only increases the temperature of the lev chromosphere . 

Longer period waves are trapped by temperature gradients, and although 
their amplitudes can increase by constructive interference, they transr.it 
little energy. 

Internal gravity waves are probably present in stellar atmospheres, 
ted while they may be important dynamically (line broadening, accelera- 
tion of steady flows), they probably do not contribute directly to heat- 
ing. 


Magnetic fields have recently received a great deal of interest, 
largely as a result of the discovery of very strong fields and the 
inhomogeneous structure of the corona. Progress has proceeded by in- 
voking mechanisms previously considered for solar active regions and 
flares: Alfven waves dissipating motions excited by convective buffet- 

ing of flux tubes, anomalous resistivity dissipating field aliened cur- 
rents induced by twisting of the field in and below the photosphere , 
and reconnection of field lines in regions of strong field gradients 
produced by the flux tube twisting and relative motions of different 
flux tubes. 

For the new, magnetic field associated heating mechanisms, ev*n 
more than the older, wave heating mechanism, small spatial structures 
are of fundamental importance and we must rely on solar observations 
and theory to lead the way to applications to other stars. 
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1. 1KTR0DUCT10K 

Ladies and Gentlemen, we now reveal to you the secrets of how to create chaos 
out of order. The existence of a chromosphere or corona requires the existence of 
motions. A chromosphere or corona requires some non-radiative heat input. There 
has to be some kind of motion, either oscillatory or quasl~static, to transport the 
energy up to the chromosphere or corona. This ordered motion may be observed as 
chaos: microturbulence, macroturbulence, line asymse tries or shifts. Of course, it 

is necessary to actually compute the effects of motions on line profiles in order to 
see what will really happen. 

We review the properties, generation and dissipation mechanisms of three kinds 
of waves: acoustic, gravity and Alfven waves. These are not the only kinds that can 

exist, but they will give you some idea of most of the range of wave properties, at 
least for the low frequency waves f^r which plasma effects are unimportant. They 
are pure cases. These different wave modes are distinguished by their different 
restoring force— pressure for acoustic waves, buoyancy for gravity waves, and magnetic 
tension for the Alfven waves. Their properties are summarised in Table I. From an 
observational viewpoint, the most important properties are the relation between tem- 
perature and density variations (which change the intensity) and fluid velocity (which 
shifts the line). Acoustic waves are compressive: In propagating waves the temper- 

ature and density vary in phase with the velocity, which is parallel tc the energy 
flux. However, for standing or evanescent waves the temperature and density are 90 e 
out of phase with the velocity. Gravity waves are slightly compressive: The temper- 
ature and density vary oppositely to each other and 90° out of phase with the velocity, 
which is parallel to the energy flux. Alfven waves are not compressive: The temper- 
ature, density, pressure and the total magnetic field strength remain constant and 
the motion is transverse to the energy flux. We will come back to tho$e properties 
in more detail later. 



There are basically cue kinds of generation mechanisms: One is direct coupling 

from the convective motions to the wave motions, either inside the convection sene or 
penetrating into the photosphere; and the other is thermal overstability. There are 
only a few basic dissipation mechanisms also: Radiation can destroy the restoring 

force and damp a wave. The other major dissipation processes occur by collisions 
which diffuse momentum and energy, and produce viscosity, thermal conductivity, and 
resistivity. Sometimes, instabilities can clump the particles so that collisions 
occur with a large collective "particle" rather than an individual one. This in- 
creases the effective collision rate and enhances the diffusion of momentum and 
energy. These diffusive transport processes dissipate acoustic waves in shocks, 
gravity waves in shear layers, and Alfven waves by viscous or Joule heating. 


II* ACOUSTIC WAVES 


A. PROPERTIES 

This material is well-known, so let us quickly run through the basics. As we 
said, the restoring force for acoustic waves is the pressure. The energy flux is 

2 

r « pu 8 pu v , 

g 

and the group velocity is 

V g * Bd-N^A- 2 )**, 

where s is the sound speed. The acoustic cutoff frequency is 

N « s/2H = yg/2s. 
ac 

Acoustic energy propagation can occur only for u > (2r/200s for the sun). In 

the absence of dissipation or refraction the flux must be constant, and the sound 
speed is roughly constant throughout the photosphere and chromosphere and increases 
by a factor of 10 going up to the corona, hence the velocity amplitude will scale 
roughly as , 


Acoustic waves can propagate or be evanescent or standing. The essential difference 
is that propagating waves transport energy, but evanescent or standing waves don't. 
In propagating waves with vertical wavelength snail compared to the scale height, 
the pressure, temperature and density vary in phase with the velocity: 


6o _ u 6T _ , , N u u 

p s’ T 1 s’ P 8 

In standing or evanescent waves although the pressure, temperature and density fluc- 
tuations are of the same order as in propagating waves, they vary 90 degrees out of 
phase with the velocity ■: 


p 


u 4T 
1 s' T 


<»-» ;• * ■ - * «t> i • 



Since the energy flux is the average aver a period of the pressure times the veloc- 
ity, the average flux will he sere. Also the vertical phase velocity ef evanes- 
cent waves will be infinite, so that the motions will he in phase all the way up and 
down through the atmosphere. Figure 1 Shows a portion of the diagnostic diagram 
that Jacques Beckers showed yesterday, to remind you thst the acoustic waves occur 
in the high frequency region. Later we will come to gravity waves, which occur in 
the low frequency region. (For acre details, see Lighthill, 1978.) 
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B. GENERATION 

1. Direct Generation by Turbulent Notions 

Acoustic waves ma> i>e generated directly oy the turbulent convective motions. 

This is what is usually called the Eighth! 11 mechanism. The radiated power is roughly 

the energy, density in the turbulent motions divided by the time scale for the turbu- 

2 

lent motions, times some efficiency factor. The turbulent energy density is pu , 
the time scale is the eddy to turn over time, which is the length scale of the eddies 
divided by their velocity, t * i/u, and the efficiency factor is the wave number of 
the wave times the size of the eddy to some power, (ki) n 
Thus, the radiated power is 


P 



(U) 


2n+J 




tte exponent ft depends on the kind of admission t Xf there is monopole emmission* 
which corresponds to a mass source* then r. • 0} there is no mass source in the con~ 
vection tone* If there is dipole emission* which corresponds to a momentum source* 
which corresponds to an asternal force* then n • 1; in a uniform medium there would 
be no external force and no dipole emission* but in stars there is an external gravi 
tational field so there is some dipole emission. .Finally* for quadruple emission* 
which corresponds to the action of the Reynolds stresses* n ** 2, and that is the dor. 
inant process (Stein* 1967). For acoustic waves: 

k * w/ s and u * * u/f* 


ki * u/s * M* 

the Mach number of the turbulent motions. What you get is the familiar result that 
the radiated power is proportional to the eighth power of the turbulent velocity: 

»•*?-(?) 

The turbulent velocity that one chooses is very sensitive to the model that one take 
for the turbulence* and therefore the emission is very uncertain. But if one r.akv£ 
some crude estimates for the sun, 

p - 10' 6 , s - 10 6 , u/s - 1/4, 


then 


7 2 

F- PI '10 ergs/cm s. 


One can also use mixing length theory to see how the flux will depend on Stella 
properties. From mixing length theory 


and 


where 


0 - (f ^ t, 
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dT/dz - (dT/dz) 
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is the sups rad iabatic temperature gradient. Hence 

F - pc p ATu - pc p (g/T)** B 3/2 l 2 . 
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From hydrostatic equilibrium 
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this means that the flux decreases very rapidly as you come down the main sequence, 
end increases rapidly as you go up to the giants and the supergiants. There are sore 

problems with this* Linsky and co-workers find that the Mgll flux is a good measure 

of the chromospheric emission and they claim that the ratio of the Mgll flux to the 
total flux of the star is independent of g, which is contrary to what the Lighthill 
mechanism predicts (Basri and Linsky, 1979). Also If you look at the cool main 
sequence stars you find that the predicted flux is much less than the scaled chromo- 
spheric losses* The predicted wave flux may be increased by including effects of 
molecular hydrogen cn the specific heats and the adiabatic gradient* Just how much 
is not known. Ulmschneider and Bohn are working on that now. But as of the moment 

there is still over an order of magnitude discrepancy in those results (Schmitz and 

Ulmschneider, 1979). 

The turbulent motions may also directly excite the 11 five-minute 99 oscillation 
In a steady state the amplitude of a given mode will be determined by the balance 
between turbulent generation by the Lighthill mechanism and dissipation by turbulent 
viscosity (Goldreich and Keeley, 1977): 

3 

PU A 5 2 

-p- (k>r = vie e k , 

where 

V “ U A X * 


Hence, the energy density of an oscillation mode k will be 


3 2 3 

c k * p> * U A k * 

vhete A is the size of the eddy whose turnover time equals the oscillation period, 
u^/A * For a Kolmogorov turbulence spectrum, where 

and H is the scale height, which is assumed to be the size of the largest turbulent 
eddies which contain most of the turbulent energy, the oscillation mode energy density 

** t ( H \un f .\+n 

h -N l") («-) • 

2* Thermal Over&tability 

Overstability is an oscillating, thermal instability. There arc several kind* 
of thermal instabilities. The one that works for acoustic waves is the k~ mechanist, 
or the Eddington Valve. If you have an opacity which increases with temperature, then 
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«feQa the gas is compressed, it gets hotter, the opacity goes up. It blocks the flow 
of radiation, so heat accumulates which raises the gas pressure, which means there is 
ttiore pressure expanding the gas than would be obtained from just compressing the gas, 
and so it will have a stronger expansion than its compression ano the amplitude 
will increase* When you actually calculate the growth rates as Ando and Osaki (1975) 
did, you find that they are very slow* The time scale for a mode to grow is about a 
thousand periods, and that Is so long that the turbulent viscosity has a chance to 
destroy the overstability* On the other hand, as we have seen, the turbulent motion 
may also direetly excite the inodes* And since ve certainly see then In the sun, we 
know something Is exciting them* It ought to be pointed out that the calculations 
show that the fundamental mode is stable. It is, however, seen on the sun, although 
at a somewhat smaller amplitude than the higher inodes. If the calculations are right, 
at least the fundamental mode must be excited by some other mechanism besides thermal 
overstability* So thermal overstability may or may not work for the five minute oscil- 
lation. Some people have proposed a mechanism of Doppler shifted line opacities as 
a generation mechanism for sound waves in stellar winds. 


C. DISSIPATION 


What about the dissipation of acoustic waves? 


1. Radiation 

In the first place photons can transfer energy from the honor to the cooler 
regions of a wave which will reduce the restoring force and damp the wave. Calcula- 
tions show that about 90% of the wave energy of the acoustic waves is removed in the 
photosphere. Radiative damping also alters the phase of the temperature and density 
relative to the velocity (Noyes and Leighton, 1963). 


2* Shocks 

The other damping mechanism, of course, is shocks. / the wave propagates Its 
front steepens, and when the thickness of the wave front becomes comparable to the 
mean free path of the particles, one gets a shock. It should be remarked that a shock 
has to do with the steepness of the gradient, not with the size of the velocity. You 
can have a shock where the velocity amplitude is small compared to the sound spt-od. 

The distance a wave must travel for a crest to overtake a trough and a shock develop 
is 

«-2»ln <l*£l ifj). 

Short period acoustic waves will dissipate near the temperature minimum, but longer 
period waves with periods around the acoustic cutoff period (200 sec) will dissipate 
higher up* And the five minute oscillation which is evanescent doesn* steepen at 
all until it gets high enough to become nonlinear. The dissipation length is 


L * 


{ 


X/ (M-l) 
X 


weak shocks 
strong shocks 
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III. GRAVITY WAVES 


A. PROPERTIES 

In gravity waves the restoring force is buoyancy* which is similar to convection. 
The different thing about gravity waves is that, while for acoustic waves tr.erc is a 
natural speed, the sound speed, for gravity waves there is a natural frequency, tht 
buoyancy or Brunt-Vaisala frequency at which a blob will oscillate if displaced: 



where £ is the superadiabatic temperature gradient. In an isothermal atmosphere 


N •* (y - U^g/s. 


Gravity waves only propagate at frequencies less than this natural buoyancy frequer. ; 
and the buoyancy frequency is only real and nonzero in conveetively stable region 5. 
You cannot have gravity waves propagating in a convectively unstable region. They 
cannot propagate or bo generated inside the convection zone, only by motions in the 
stable photosphere. Gravity waves propagate energy in a particular direction, whic*. 
depends on frequency. The cosine of the angle between the flux and the vertical is 

cos 6 ® w/N. 

. jy the sun, for values that are appropriate for the granulation. 


T . - - 10-20 min and N. . *0.03* 2*/3xin, 

granulation Tmin * 

the direction of gravity wave energy propagation is 

cos', * 1/5, 6= 75°. 

Figure 2 shows what you would see if you looked at a schlieran photograf^h of gravity 
waves produced by an oscillating source. The solid lines are the wave crests and the 
dashed lines are the wave troughs. The lines intersect at the source. Energy 
Is propagating radially outward at the angle theta and the velocity of the fluid is 
also radial, parallel to the energy flux, but the phase propagates perpendicular to 
the energy, 

V //Fib* 

You sec different waves moving across the fan, while the fan extends out further and 
further with time as the energy gets out further and further. The group velocity of 
gravity waves is 


sin 0. 
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For solar granulation N/k * 1 km/s. Figure 3 shows the group velocity for gravity 
(and acoustic) w^ves as a function of frequency and direction. Notice that one would 
observe low frequency, long wavelength waves' first. The phase relations between tem- 
perature, density, pressure and velocity for low frequency gravity waves are similar 
to evanescent acoustic waves: 


«£. « , H «I « _i «E . -ff- H << H. 

p 8*1 S*p ®^H S 8 

Pressure fluctuations in phase with ve3ority but very small 
density are 90° out of phase with velocity (Pittway and Hines, 
Chapter 4). 


, whil<? temperature and 
1965; Lighthill, 1978, 


B. GENERATION 


Gravity waves are ubiquitous in the atmosphere of the earth; they are produced 
by any slow motion. Therefore, they should also be present in the Sun. There are 
two common generation mechanisms. 

1. Penetrative Convective Motions 

One of the main ways gravity waves will likely be produced is by the penetrative 
convective motions, >ie can think of the penetrative convection as blobs pushing 
on the boundary of the stably stratified layer. The amplitude of the wave produced 
will be comparable to the amplitude of the penetrative motion, 


u * u . . 

‘‘Wave 1 1 -penetration 1 

This has been verified in laboratory experiments (Townsend, 1966). However, since 
only frequencies that arc less than tie buoyancy (Brurn-Vaisala) frequency can propa- 
gate, only that part of the penetrative convective power that satisfies * ® k. V < N 
will contribute -o the production of gravity waves. This mechanism is similar to the 
Lighthill mechanism, but with an efficiency rr r one. 

If you make a rough estimate for the soi^r granulation, taking a velocity of 
1 km/sec and the appropriate length scales, then 

F - pu^V * 3x10 ^ x 10*** x \ 10^ * 10® erg/cm^s. 
gz d 

This flux will, however, be greatly reduced by the strong radiative dissipation c f 
gravity waves, which we discuss below. 

2. Shear 

Gravity waves :an also be produced by the shear that will 'rise from thv super- 
granule motions Supergranule fi~*« have a cellular structure. Conservation of mass 
requires that a gradient of the vertical momentum flux produces a horizontal momentum 
flux. Braking is large in the photosp*vre id produces a large horizontal flow there. 
Even though the horizontal momentum flux is small in the chromosphere, the chromo- 
spheric horizontal velocity is large, b- cause of the small density. The horizontal 



supergranule flow is observed to decrease from * 0.8 km/ & In the low photosphere to 
_ •. 0.4 km/ s in the low chromosphere and then increase to - 3 ka/s in the mid-chromo- 
sphere (November, et al. 1979). Vftiere the else of shear becomes comparable to the 
buoyancy frequency. 


d U H <*> 

da 


>_ min £ 


k,n 


cool 


] 


(where T eoo ^ is the radiative cooling time) the shear layer becomes unstable and 
radiates gravity waves. Most of the energy is radiated near 


k « U/fT l' H , 


and the growth times are of order 

Y * 10~ 1 dl? H /dz 

(Lindze~, 1974). This mechanism will operate in the low photosphere where the cool in 
time is short and in the high chromosphere where the shear is large. 


C. DISSIPATION 

Now do gravity waves dissipate? 


1. Wave Breaking 

Gravity waves steepen, but instead of forcing a shoe**, front they form a thin 
shear layer, where the fluid velocity changes direction over a very short distance. 
When that shear becomes comparable to the buoyancy frequency, dV/dz * N, turbulent 
will develop along that wave front. Small scale motions are produced which dissipate 
the wave motion and damp the wave. To find the condition on the wave amplitude for 
breaking to occur, we need to calculate du/dz . Let u * (u, 0, w) where u is the :.c*. 
2 ontal and w the vertical component of the velocity. For gravity waves the Boussir.os 
approximation hoiJs, so 

V . u *0, 


which implies that 
The wave equation for gravity waves is 


dw . . 

— B - i k„ u. 
dz H 


80 


7 ? ♦(? - , ) k " 2 “■ 0 • 


Hence, the condition for gravity wave breaking is 


»- * {7 - *)¥ > 


1 
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2* Radiative Draping 

As we mentioned, there is very severe radiative damping of the gravity waves. 
Radiation tends to rake a wave isothermal, which destroys the buoyancy restoring 
force. The radiative cooling rate is 


y R 



16 *pT 


pc, 


-(* 




1* j ^ optically thin 


1 

3 



optically thick. 


where ic is the inverse of the photon mean free path (Spiegel, 1957). The opticallx 
thin damping time cc/16- d\ increases rapidly with height and exceeds 10 nin above 





700 ka above • 1. Heaping is severe when the radiative cooling time is shot ;er 

than tk wave period. The long period, abort wavelength waves are the most highly 
daoped. Figaro 4, free Barbara Mi hales* thesis (1979), shows the fraction of flux 
reoainiag at each height for several wavelengths and periods. In order of magnitude, 
the flux decreases by 10^ between the bottom of the photosphere and the mid -chromo- 
sphere , for paraneters appropriate to granule produced gravity waves. On the other 
hand, gravity waves generated by penetrative motions a few hundred kilometers above 
the bottom of the photosphere suffer reduced radiative dissipation and can transmit 
22 of their flux to the chromosphere. Since the radiation fc'iich produces the dtefip.: 
of gravity waves also allows for easier penetration of the convection, the wave flux 
reaching the chromosphere from penetrative convective motions high up in the photo- 
sphere is only slightly less than that from the bottom of the photosphere with its 
severe radiative damping. 


3. Critical Layers 

Gravity waves have one other property which is very different from sound waves: 
they can interact very strongly with the mean horizontal fluid flow. In particular, 
they can give the fluid a horizontal acceleration. A layer where the horizontal 
fluid flow velocity is equal to the horizontal phase speed of the wave, is a 'Vri t i 
layer.” Here in the fluid frame, the dopplcr shifted frequency goes to zero, 
w-k.u -> 0, so the waves propagate along the critical layer, cos-Z * (_-k.u)/N - 0 , 
and their energy is absorbed. Such a critical layer can arise from the horizontal 
super granulation flow or may be produced bv absorption of horizontal gravity wave 
momentum. For the sun, the horizontal gravity wave phase speed is about 1.5 kn/s, 
which will match the horizontal supergranulation flow somewhere in the mid-chromo- 
sphere. The transmission through a critical layer is 


trans. 

r 

r inc. 


exp 


{-H ri - if} • 


where the Richardson number, 

Ri - K 2 /(dl/<5z) 2 , 


is the order of 50. Hence, the absorption at critical layers is very large (Looker 
and Bretherton, 1967; Acheson, 1976). Such critical layers will occur in the chromo- 
sphere and will produce localized dissipation of gravity waves there. 

The important thing to emphasize about gravity waves is that because they pro- 
pagate energy mostly horizontally and because of radiative damping, one does not 
expect too much heating from then. But they are a source of chaos and may contribute 
substantially to any microturbulcnce, because they are certain to be there. Their 
horizontal wavelengths are comparable to granules, < lMm, and their vertical wave- 
lengths are only 1/4 as large, which is comparable to the scale height or smaller. 



IV. ALFVEN WAVES 


A. PROPERTIES 


Roving on apace, we now come to Alfven waves* Here the restoring force is 
sagnetic tension. The group velocity is the Alfven velocity. 


Unlike the sound speed, the Alfven speed increases substantially between the photo- 
sphere and the corona, by a factor of 10*. The flux is 


F 


- *u 2 


a 





and is parallel to the magnetic field* In the absence ox dissipation or refraction, 
the flux remains constant, so the velocity a*nplltude scales as 

„ . ,-U* ,-»« 

The velocity is perpendicular to the magnetic field, i.e. transverse to the direction 
of energy propagation. Hence Alfven waves act like a vibrating string- There is r.c 
compression, 

£T « fp ® £f * 0, 


so there are no opacity changes. Only Doppler shifts affect the line profiles. 

The total magnetic field strength, + 6B|, is constant, so the wave is polar- 
ized in part of the arc of a circle, along which the magnetic vector swings back a:\c 
forth (see e.g. , Bazer and Fleischnan, 1959; Barnes and Hollweg, 1974). 

Most people up until recently have considered uniform magnetic fields. But ve 
know that in the corona the magnetic field is very inhomogeneous. Luckily it turns 
out that waves in inhomogeneous fields are rather similar to the internal waves in & 
uniform field. If you consider a thin flux tube surrounded by plasma in a weaker 
field, several different modes occur (Roberts and Vebfc, 1978; Kilson, 1979; Ventztl. 
1979a). First, there is an azisymetic mode which is just the torsien.i] Alfven vavt , 
propagating along the flux tube at the Alfven jpeed inside the tube, c ■ a... Seccrc, 
there is another axisyxnmet ic mode which is like the slow mode, a sound wave propage 
ing along the flux tube with 

2 2 2 it 2 2 * 
c * a 1 s 4 /(a i + s i ). 

It has a short wavelength and its amplitude is concentrated at tne surface of the tube. 
For these waves to excite waves outside the tube, th*: phase velocity inside the tub* 
would have to be greater than the sound or the Alfven speed outside. In general that 
is only possible for these modes when the inside density is less than the outside 
density. Third, there are other modes which arc not axisymmetric , and act like vi- 
brating string modes. Their phase velocities are essentially an average of the Alfven 
speed inside and outside the tube. 



c 2 • (B^ 2 ♦ B e 2 )/ 4» (P t ♦ o e >. 

these modes can have a resonance where the phase velocity of the tube mode equals 
the local ven speed. At that resonant point the amplitude of the perpendicular 
tto the magnetic field) velocity and electric field becomes large. Other perturbed 
quantities are unaffected by the resonance. Thus, even in an inhomogeneous corona 
the waves are similar, although not identical, to those in a homogeneous corona. The 
main difference is that for these tube or surface wave modes there exists a resonance 
at places in the flux tube where c * a. 

B. GENERATION 

What about Alfven wave generation? 

1* Convective Motions 

Alfven waves can be generated by the convective mot * ons , similar to the Lighthill 
mechanism for the sound waves. But* because the magnetic field channels the motions, 
monopole rather than quadrupole emission occurs. The radiated power is 

P • ^ (kt). 

For Alven waves, 

. u u 

k “ I ’ - “ 7 ’ 

SO 

k£ - a " 

Thus 

(Kulsrud, 1965: Kato 1968). Rough estimates for the sun p. 10~®, u * 10*, a - 10® 
predict an Alfven wave flux from strong field regions of 

F - Pt - 10® erg/cm 2 s. 

However, anything which jiggles a magnetic field line will also generage Alfven 
waves, so granules will generate Alfven waves and supergranules will generate Alfven 
waves. Granules have larger velocities and so are more important. They produce a 

flux - 

F * pu a. 

For typical granule velocities of 1 km/s, 

F* 3x10 2 x 10*® x 10® * 3x10® ergs/cm 2 s. 

(See, however, Hollweg, 1979.) 

What we are really interested in is the average flux, so we have to Include the 
fact that the flux tubes will spread with height and that the Alfven speed increases 
with height. The waves produced by granular motions have fairly long wavelength so 
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they see the change in Ai.*veo speed roughly as a discontinuity* In this case the 
transmission coefficient is 


,.1Vl 


(■i sf 


* 4 — 
*2 


10 


-3 


The average flux is 


photo 


A 


corona 


photo 

a * 

corona 


where A is the flux tube area* Since the magnetic flux is constant along a flux 

- 1/2 

tube* BA * constant , so Aa »p , Thus, the average flux will be 

F » (^corona^photosphere) 1 ^ 2 



5 2 

* 3x10 erg/cn s, 

which is fairly substantial, enough to heat the corona. 

1 have not made a distinction between Alfven waves and magne toby cor dynamic wave*. 
In a strong field the Alfven and the fast mode are similar and the slew node is an 
acoustic wave propagating along the flux tube. 

2* Thermal Owe ratability 

Alfven waves may also be generated by thermal overstability. This is not the * 
mechanism, but the Cowling- Spiegel mechanism (Cowling, 1957; Moore and Spiegel, 1966 ). 
Here buoyancy acts as a driving force which tends to destabilise the system and make 
it depart from equilibrium* Magnetic tension acts as a restoring force which tends 
to bring it back to equilibrium and radiative transfer decreases the destabilizing 
effect of the buoyancy force. Since there will be less destabilizing effect on the 
way back than there was on the way out from equilibrium, the magnetic tension will 
return the system to equilibrium faster than it departed, and the wave amplitude 
will grow. Because* the buoyancy must be destabilizing this mechanism works only in 
convectively tins table regions. One can calculate the grew-*' rate by equating the 
rate of working of this buoyancy force with the kinetic energy of the waves (Parker, 
1979). The buoyant force is 

h j, AT 

c Apg * Pg “J. 

The temperature fluctuation is the temperature difference between the adiabatically 
displaced fluid in the wave and mean temperature at its displaced level, reduced by 
diffusive radiation cooling: 

AT « AXB ^ — , 
coo] 

When* A is the wove amplitude, > is the wave length, f, is the superadiabatic temper u- 
ture gradient, t is the period, and the diffusive radiation cooling tire is 



X 4 


*V 

c 


!gas . PC - A pc 


cool cfc - E 

1 mfp & rad 


A * 3 • 

at *coT J 


Which la assumed to be much greater than the period* The growth time is the time 


1 2 

it takes the buoyant work, vF^, to supply the wave energy j ov , Where v * Aa * AA/t 


Thus 


y ’ 1 V F b - i p» 2 . 


-|2 . !S2£1 

1 ~ 2 2 
Tw T , t , , T 

cool eddy cool 


Bow does this growth rate depend on stellar properties? For an opacity 


■ • W* ■ »*•’ »». 


hydrostatic equilibrium gives 


r**Z « g'V 6 . 


SO 


g/T. 


The wave frequency is 
and the cooling rate i> 

From mixing length theory 9 


2 J/ 2 n 2 / i 2 t, 2 -*6 < _+7,-2 

w*a/A - B/cA * B g T a , 


r 1 , 

cool 


" ..-‘V 1 X- 2 , 


PC Tk*>. 

p 


-2 g£ 

t e -fit- 

eddy T 


te) 


2/3 


- g 1 ' 6 T 5 ' 25 , 

where t ° H 8 T/g. Hence the growth rate varies with stellar gravity* surface temper- 
ature 9 and magnetic field as 

T - g 2, V 25 „- 2 . 


C* DISSIPATION 

1* Viscous and Joule Heating 

How do Alfven waves dissipate? Alfven waves don*t steepen and form shocks, 
because the Alfven speed is Independent of wave amplitude, since the magnetic field 
strength is constant and they are not compressive. They can still dissipate by par- 
ticle collisions which produce Joule or viscous hocting, However* the damping lengths 



The Joule and 


are large, although they may be comparable to coronal loop dimensions, 
viscous heating rates are 


Tj » n J 2 • nc 2 k 2 i«B> 2 /(4*) 2 
and 

» v k 2 u 2 . 


The damping rate is 
where 


and the damping length is 


Y - Q/2E, 

E « pu 2 * 6B 2 /4i! , 
L * a/y. 


For Joule heating 

i * ci. 3, 2 2 
Lj a a /ric u 


where the resistivity is 


For viscous heating 


. 10 ™ »V c 


" ‘ ■ IT 1 ■ ‘ • “' V3/2 - 


Ly a 2pa 3 /vu 2 


• !<>” „->«T - 5 'W„. 


where the viscosity is 


v* \ " Ve 2/v coll * 10 " 15t 5/2 ' 


Viscous damping is more important than Joule heating in the corona and visaversa in 
the photosphere and chromosphere. If the fields are weal:, then the Joule dissipation 
will be significant in the photosphere. Hence Alfven waves can only get through fron 
the photosphere up into the corona in strong field regions. That used to be a serious 
problem before we knew that fields come in little patches of high field strength. 

Now it is not. For typical coronal loop parameters (n“10^cm~^,T-2xlO^°K, B -100G, 

L ~10^cm) the viscous damping length is 

ly * 10VCK, 

which it comparable to the loop length for short period waves. 



Because the corona is inhomogeneous, the Alfven waves are 
and have a resonance where the tube phase speed is equal to the local Alfven speed* 

In this resonant region the wave amplitude is large. Hence largo currents and large 
Joule heating will occur in the narrow resonant layer. The rate of heating is con- 
trolled by the rate at which energy can flow into that resonant region, and has been 
calculated by lonson (1978) and Uentzel (1979b): 

y ■ irwk^Ar (Aa 2 /a 2 ) (2 + P i /r e ♦ P e /P t )~ l * 

They thought this was the rate at which waves were radiating energy away. It isn't. 

It is the rate of energy flow into the resonant region (Hollveg, 1979). There is a 
problem of how the energy released in this very small resonant region volume is 
transferred to the rest of the large coronal volume where it is needed. Nobody has 
figured out how that is done. This is a problem for just about every type of Alfven 
wave dissipation, except for the Joule and the viscous dissipation. The reason is 
that in order to rapidly dissipate the Alfven wave energy the currents must be clumped, 
and so the dissipation occurs in a small regio:. 

2. Mode Coupling 

In the presence of inhomogenie ties the Alfven wave will couple to other wave 
modes. Coupling will be large between wave modes whose wave vectors* difference is 
comparable to the inverse of the inhomogeniety scale length. The coupling ratio 
between two modes is roughly 

law. r 1 . 


really tuba 




where Ak is the difference in k between the two modes and L is the length scale of 
the inhomogenie ty (Melrose, 1977). In a strong field the major coupling occurs be- 
tween Alfven and fast mode waves that are propagating in the direction of the magnetic 
field, because both of them are propagating nearly at the Alfven speed, so they will 
stay together. The difference in wave vector is 


Ak B L ( c ) ” 


u> Ac 
c c ’ 


where c is the 
direction 


phase speed of the mode. For propagation close to the magnetic field 


2 . 2 _2 
a + s 6 , 


2 

c 




a 2 (1 - e 2 ). 



sssft as 


Ac _ C +* C A 


« e 2 /2. 


(“ij ® 2 ) -(>- > 2 ' 2 ) 


and 


Ak/k * e‘/2. 


Thus the coupling ratio is 

« 2 (jo.) - 1 e" 2 . 

When 6 < (w/H^) * where 2 ^ is the ion- cyclotron frequency* ion cyclotron effects 
increase Ak. So the maximum coupling occurs at this critical angle and is 

I AkL I" 1 * (kL)" 1 C^/-) 


(Melrose, 1977). There will not be much coupling if the waves are not propagating 
along the field direction, nor will there be much coupling between the Alfven and 
slow modes. The fast mode, if it gets any energy, will form shocks and dissipate * 
so that this is a round about way in which the Alfven waves can dissipate their ene 


3. Alfven Wave Decay 

The Alfven waves can also decay. If a set of waves satisfy the resonance condi- 


tion, 


bl 

O 


“ 2 , 


k 

*»o 


h 


* *2 , 


which is essentially energy and momentum conservation, then one wave can decay into 
two. In this case a forward moving Alfven wave can decay into a backward moving 
Alfven wave and a slow mode pressure wave, at the rate 

(Kaburaki and Uchida, 1971). For a broad spectrum of incident waves, only those that 
stay in resonance for a decay time, i.e. ha'*e Ak/k * Aw/u, < >/u , can decay. So 
the decay rate is 


(Sagdeev and Galeev, 1969). There seems to be some disagreement between the calcu- 
lated decay rates and the fact that one sees the Alfven waves in the solar wind at 
the earth. 



4* Current Dissipation 

Finally* I want to talk a little about current dissipation. Currents are pro* 
duced not only by. Alfven waves* but by an; twisting notion of the magnetic flux tubes* 
for instance a quasi-static twisting motion. Current or magnetic field dissipation 
is a diffusive process due to single or collective particle collisions. The charac- 
teristic resistive diffusion time scale is 




4irl 2 /nc 2 . 


4 

For typical coronal parameters ^10 yrs* too long to be significant. This Joule 
dissipation time can be reduced either by reducing the width L of the region through 
which the currents flow or by increasing the .esitivity by increasing the effective 
collision rate. If some instability or resonance filaments the current so the current 
density is high in a small region* then there can be significant dissipation of the 
currents. If that occurs and if the current density J * n^ev becomes large 

enough so that the drift velocity approaches the electron thermal velocity then sub- 
stantial numbers of electrons will tend to run away and generate several different 
types of electrostatic waves. These waves bunch the ions, so that the electrons 
collide with the electric field of a large collective charge rather than that of a 
single ion. This scattering of electrons by the waves increases the effective colli- 
sion rate* the rate of momentum transfer and hence the resistivity. The enhanced 
resistivity due to electron scattering by plasma waves is called "anomalous resis- 
tivity*" and since it occurs in conjunction with current fllamentation will shorten 
the resistive diffusion time tremendously (Papadopolous* 1977; Rosner et al* 1978; 
Hollweg* 1979). Also if the current density becomes large it will develop large 
shears or gradients in the magnetic field* which will lead to tearing mode insta- 
bilities. Parallel currents attract one another and tend to clump. The clumping 
of current produces a fluid flow that force* the sheared magnetic field into X-type 
neutral points. Filamented currents are produced with small enough length scales 
so that the classical* Coulomb collision, resistive diffusion time becomes small and 
the magnetic field can tear and reconnect and the currents can dissipate (Drake and 
Lee* 1977). This mechanism has been Invoked for the violent energy release In flares 
(see Spicer & Brown* 1980). However, shorter wavelength tearing modes distort the 
field lines more* which produces a greater restoring force, and also have a smaller 
volume of magnetic energy they can release* so they may produce a more tranquil quasi- 
static heating appropriate for coronal flux tubes. The tearing Instability has a 
lower threshold than the current driven instabilities vtUch lead to anomalous resis- 
tivity. To get significant current dissipation by any oechanism* the dissipation 
must occur in such small volumes that the transfer of tbe resulting heat to the rest 
of the corona is a serious problem. 



V* CONCLUSION 


In conclusion, there are many kinds of different vave motions in the sun; acous- 
tic, gravity, Alfven waves, and other kinds of magneto-acoustic-gravity wave modes; 
maybe even higher frequency waves like whistlers. All of these ordered motions may 
contribute to the chaos observed in stellar atmospheres. In order to develop diag- 
nostics, somebody has to take self-consistent calculations of these waves with the 
right velocities, temperatures, densities and pressures and calculate the effects 
on the line profiles of each wave mode. 

To summarize the major roles of the waves we have discussed today: Acoustic 

waves can heat the low chromosphere but not the corona. The evidence is partly ob- 
servational: the observed non thermal line widths due to waves in the upper chromo- 

sphere are too small, and also theoretical: Increasing the driving amplitude at the 
bottom of the atmosphere only increases the dissipation of the acoustic waves in the 
chromosphere, but doesn’t increase the flux through the transition region to the 
corona. 

For gravity waves the motion is mainly horizontal. Their amplitudes will be 
comparable to the amplitudes of the penetrative convection (granulation)- They r.av 
contribute to the observed microturbulent velocities, but they are unlikely to be 
important in the heating. 

Alfven waves will dissipate primarily by highly clumped currents in very small 
regions, so there is the problem of how to get that energy from the small volume 
where the dissipation occurs to the larger volume of the corona. The nice thing about 
Alfven waves is that they are observed in the solar wind, and they seem to be impor- 
tant in providing an energy and momentum input to the wind. Someplace between the 
photosphere and the earth, where the wind is observed, those Alfven waves must be 
produced. 





TABLE I. 
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GRAVITY 
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PROPERTIES 
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Pressure 

Buoyancy 

Magnetic Tension 


u // F // k 

«/ / F 1 k 
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1 

Propagating: 

cos9 * uj/N 

V g » a 


6T ~ §2. ~ u/b 

T p p 

it * JL. 4 «u/b 

P k H S 8 

6T * 6p * 6r 8 0 


Evanescent 


| B o + $Bl * const 


6T . 6p . ip „ . ij 
T p P s 

6p . 6T . . • 

“ ~=T i U/8 
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GENERATION 

Convective Motions ' 

Penetrative 

Convection 

Convective 
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Thermal 

Over stability 


Thermal 

Overstability 

DISSIPATION 

Shocks 

Shear 

Viscous and 
Joule Heating 
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collisions 

radiation 

Critical 

Layers 
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Instabilities 
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Radiation 
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ATTACHMENT 


SOME EFFECTS OF SIRONG ACOUSTIC WAVES ON STRONG SPECTRAL LINES 

Pierre Gouttebroze (Laboratoire de Physique Stellaire et Planet ai re, 
Verri£res-le-Buisson , France) 

John Leibacher (Lockheed Palo Alto Research Laboratory, Palo Alto, CA) 


We have studied the formation of optically thick lines in time 
dependent , non-linear hydrodynamic model of the solar chromosphere. 
Models of the 200 second, chromospheric oscillation indicate that the 
emission peaks of self-reversed profiles such as those of MgH and 
Call are firmed at very different depths depending cm the phase of the 
oscillation, while the central absorbtion feature is emitted at a very 
nearly constant mass depth. The figure shows the Mgll k line emitted 
by the mean atmosphere, including "microturbulence" (triangles) and 
the mean profile (squares) . In addition to the substantial intensity 
increase interior to the emission peaks, one should note that the peaks 
are broadened only towards line center; i.e. the intensity fluctuations 
are symnetric outside of the peaks and strengthen more than they weaken 
within the peaks. A more detailed version has been submitted to the 
Astrophvsical Journal. 
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